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UNIVERSITY OF SOUTHAMPTON 


ABSTRACT 

FACULTY OF ENGINEERING 
DEPARTMENT OF AERONAUTICS AND ASTRONAUTICS 

SOME ASPECTS OF WIND TUNNEL MAGNETIC SUSPENSION SYSTEMS WITH SPECIAL 
APPLICATION AT LARGE PHYSICAL SCALES 

by Colin Paul Britcher 

Despit;e over a quarter of a century of development, wind tunnel 
magnetic suspension and balance systems (MSBSs) have so far failed to 
find application at the large physical scales necessary for the majority 
of whole model aerodynamic testing. Recent developments, such as the 
cryogenic operation of wind tunnels and advances in superconducting 
electromagnet technology have greatly reduced the apparent cost of a 
large, or perhaps more specifically, high Reynolds Number MSBS (LMSBS) . 
Many difficulties remain, however, and three are addressed in this 
thesis. 

A powerful method of magnetic roll torque generation is essential for 
any LMSBS. Two variants of the new Spanwise Magnet scheme are studied 
herein. Spanwise Permanent Magnets are shown to be a practical method 
and are experimentally demonstrated using the Southampton University 
MSBS, though precise evaluation of maximum torque capabilities has not 
been possible. 

Extensive computations of the performance of the Spanwise Iron 
Magnet scheme indicate potentially powerful capability, limited 
principally by current electromagnet technology. Some experimental 
verification of the computed performance at low applied field levels 
is presented. 

Aerodynamic testing at extreme attitudes is shown to be practical 
in relatively conventionally configured MSBSs. Preliminary operation 
of the Southampton University MSBS over a wide range of angles of 
attack is demonstrated. 

The impact of a requirement for highly reliable operation on the 
overall architecture of LMSBSs is studied. It is shown that the 
system's cost and complexity need not be unduly increased, provided 
certain of its unique characteristics are exploited. 
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INTRODUCTION 


1.1 Magnetic Suspension and Balance Systems (MSBS) with emphasis 
on Large Systems (LMSBS) 

The first MSBS for wind tunnel application was constructed in 
Prance by ONERA and first reported by Tournier and Laurenceau in 
1957 (1) . This system successfully demonstrated 5 component 
control of a range of simple models in work such as wake and drag 
studies up to intermediate supersonic Mach numbers. The potential 
advantages of wind tunnel testing without mechanical model supports, 
including the elimination of support interference and more straight- 
forwcird dynamic testing, fired the imagination of researchers in 
numerous institutions and quite soon many systems conceptually 
similar to the early ONERA system were under construction. 

A total of eleven MSBSs of various designs could be considered as 
constituting this first wave of effort and are tabulated briefly 
belov/. All bar one of the institutions involved presented papers 
at the First International Symposium on MSBSs in 1966 (2) . 


Table 1.1 

MSBSs pre-1970 

(listing in . 

alphabetical order) 

Institution 

Degrees of 
freedom 
controlled 

Approx . test 
section size 

/ 

a/ cm 

Application 

AEDC/NASA Langley 

5 

30 

Wake studies /R&D 

Univ„ Michigan 

1 

4.5 

Low Re. No. sphere 
drag 

MIT (A) 

5 

10 

Static/dynamic 

MIT (B) 

5/6 

9-15 

Static/dynamic/ 

R&D/Magnus 

NASA Langley 

1 

11 

R&D 

ONERA (A) 

5 

5.5 - 8.5 

Drag/base pressure 

ONERA (B) 

5 

26 

Base pressure/heat 
transfer 

Princeton Univ. 

3 

13 

Wake studies 

R.A.E. Farnborough 

5 

18 

Sting effect/Magnus 

Univ. Southampton 

6 

13 - 18 

Static/dynamic/ 

R&D/Magnus 

Univ. Virginia (A) 

3 

9 

Cone & sphere drag 


(data from Ref. 3 and various) 
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Two further MSBSs were coiranissioned in the early 1970' s, 
including the first system utilising superconducting electro- 
magnets. A Second International Symposium was held in 1971 (4) . 

Table 1.2 MSBSs post-1970 

Oxford Univ. 3 8 Low density sphere drag 

Univ. Virginia (B) 3 13.5 Dynamic stability /R&D 

Contemporary with these developments there was an increasing 
recognition in the aerodynamic testing community that the majority 
of existing wind tunnels were of too low a Reynolds number capability 
for the test purposes to which they were put. Further, it was clear 
that the bulk of production testing over the coming years would be 
at Mach numbers no higher than low supersonic ( <2.5), thus 
including the regimes where Reynolds number and support interference 
tend to have their greatest influence. Interest in MSBSs waned, 
due principally, it is thought, to the realisation that scaling of 
existing technology MSBSs to the large wind tunnels then known to be 
necessary, would be impractical for one reason or another, for 
instance E/M power requirements, with associated high capital and 
running costs, or lack of adequate roll torque capability. Indeed 
the modern day resurrection of the cryogenic wind tunnel concept, 
now universally accepted as the best approach to high Reynolds number 
testing in the crucial transonic regime, came about as a part of 
efforts to reduce the size of wind tunnels necessary for a high 
Reynolds number MSBS. 

Relatively recently however, various factors have acted to increase 
the apparent viability of a LMSBS. Concern over the accuracy and 
validity of support interference corrections is now widespread, 
particularly and rather paradoxically in the case of the new breed 
of high Reynolds number cryogenic wind tunnels, since these tunnels 
are generally operated at high pressures. 

MSBSs of the type considered to date are inherently capable of 
dynamic testing, which otherwise presents severe difficulties with 
mechanical model supports, and several small systems have been 
successfully used for this purpose. 

A related capability is that fixed model attitudes may be 
selected or changed rapidly and with little restriction, potentially 
manifesting a highly productive system. 
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Further potential advantages of a LMSBS are thought to lie in 
the areas of high angle of attack (extreme attitude) testing, where 
support interference can become rather intractable, and as an 
alternative solution to the difficulties involved with operation of 
complex conventional strain gauge force baleinces in the low 
temperature environment of the cryogenic wind tunnel. 

The above factors, coupled with the recent rapid developments 
of large scale superconducting E/Ms for power generation and nuclear 
physics research indicated that revived effort aimed at developing 
a LMSBS would be worthwhile. This effort was and is principally 
directed by NASA Langley Research Center with the intention of 
developing a LMSBS at that institution. 

1.2 Status of MSBS development circa 1978 

Many of the aforementioned small MSBSs were in no way intended to 
be representative of large systems , rather being specialised designs 
tackling specific and specialised problems . Only two systems are 

known to have been constructed with scaling to larger systems 
specifically in mind, namely MIT(B) and UVa(B). Nevertheless much 
existing technology was thought to be applicable to a future LMSBS 
and some of the relevant design features will be reviewed here. 

The reader must carefully note that the comments below refer only 
to the state of developments at the time that this project commenced 
(late 1978) . Considerable developments are reported in later 
chapters herein and it must be assumed that corresponding or 
complementary developments have taken place elsewhere. 

1 .2,1 E/M c onfigurati ons 

Most current MSBS E/M configurations can be considered as falling 
into one of three broad categories as follov/s (See Fig. 1.1) : 

I) Geometrically orthogonal simple E/Ms 

e.g ONERA (A) & (B) , Univ. Southampton 

II) Convoluted E/M e.g MIT(B) 

III) Orthogonal force , isotropic model e.g UVa (A) fi (B) 

Class I above is characterized by separated E/Ms often of simple 
solenoidal form arranged in a relatively straightforward configura- 
tion around the test section, generally with the total number of 
E/Ms equal to or only marginally greater than the number of degree-s 
of freedom controlled. 
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class II features E/Ms, other than the axial field E/Ms, usually 
of non-solenoidal form, arranged in relatively complex ways around 
the test section with the intention of achieving relatively high 
uniformity and symmetry of all fields as well as high magnetic 
efficiency. The number of E/Ms is not strongly influenced by 
the number of degrees of freedom controlled. 

Class III is inherently limited to three degree of freedom 
control (zero magnetic torques), being originally developed to 
perform dynamic testing. This class cannot be considered viable 
for a LMSBS where principally static testing is to be performed 
since there is no possibility of selecting arbitrary model attitudes. 

Numerous subclasses of Class I exist, mainly classified by the 
E/M disposition in a cross section through the test section, such 
as "L”, "V", etc. It seems certain that any future LMSBS will 
exhibit considerable E/M symmetry for various reasons, including 
E/M efficiency and system reliability, so only subclasses such as 
the classical "+" and "X" need be considered. 

Considerable efforts have been expended over many years in 
attempts to derive "optimum” E/M configurations but early design 
studies of superconducting E/Ms for a LMSBS have indicated that 
E/M configuration and geometry is likely to be forced predominantly 
by technological limits of superconductors and by requirements for 
roll control (5) . All configurations studied herein fall into 
Class I above, but further discussion of E/M configurations may be 
found in Ref. 6. 

1.2.2 Power supplies 

The class of supply almost universally used to date has been the 
phase controlled multiphase rectifier, with thyristor switching 
elements in all but the most venerable units. This type is easily 
scaled to high powers. One institution is known to have experimented 
with rotating machinery for primary power supplies (MIT) although it 
is difficult to achieve high output slewing rates with conventional 
generator designs. 

The high reactive powers that would be experienced in dynamic 
operation of LMSBSs make some form of energy storage in the supply 
most attractive, in order to raise the system's overall energy 
efficiency by regenerating power from the load E/Ms. The only MSBS 
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operated with superconducting E/Ms to date (UVa(B)), utilised power 
supplies that were specially constructed but can now be regarded as 
of relcitively conventional design, a notable feature being 
capacitive load energy storage and regeneration (7) . 

Bipolar operation will be essential in at least some of the 
suspension E/Ms in any realistic LMSBS but relatively few small 
MSBSs have utilised bipolar supplies, mostly using gravitational 
or D.C, current biasing to achieve bidirectional control field 
capability. 

No power supply of the required characteristics for LMSBS 
application is known to have been demonstrated at design power 
levels remotely approaching the LMSBS requirement. 

1.2.3 Position sensors 

Optical methods of model position sensing have been employed in 
the overwhelming majority of cases. In turn, the majority of these 
methods are detail variants of the simple analogue shadow movement 
detection systems described in Appendix 1. These systems are 
inapplicable to LMSBSs because of the lack of immunity to light path 
degradation.. ONERA developed a target scanning/tracking system 
using vidicon tubes (8) , which is thought to represent a broad class 
of sensor suitable for LMSBS use, albeit with the magnetic field 
sensitive vidicons replaced with more modern solid state CCD or CID 
imaging arrays. Oxford University also used a scanning method (9) 
more closely approaching conventional T.V. techniques although the 
system was only used to detect the centroids of small spheres. 

The MIT Electromagnetic Position Sensor (EPS) (10) has often been 
viewed as a first choice for a LMSBS position sensor since the non- 
optical approach accrues, in principal, considerable flexibility and 
versatility. At the time of writing that view must be considered 
still valid since the system has amassed hundreds of hours of operation 
at small scale, although the problem of electromagnetic noise pickup 
is severe enough to demand very careful design for a large system. 
Numerous advanced optical methods have been proposed and some sub- 
jected to preliminary evaluation. None has so far been demonstrated 
on a comprehensive and fully representative basis at any scale. 

1.2.4 Control systems 

These have almost universally consisted of simple analogue feed- 
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back networks. Digital controllers have recently been demonstrated 
with an idealized single degree of freedom system (MIT) but no 
comprehensive application has been attempted. There does not 
appear to be serious difficulty with the control hardware for a 
LMSBS, the technology of reliable high speed digital flight control 
computers developed for civil, military and space applications 
appearing adequate. 

Existing control algorithms generally do not encompass the aero- 
dynamic characteristics of the model, nor the magnetic or inertial 
cross-couplings between. model degrees of freedom which exist in all 
MSBSs, rather controlling each degree of freedom separately with 
aerodynamic loads being small with respect to model mass. This 
approach is unlikely to yield the best performance from a LMSBS 
since, firstly, the aerodynamic loads may be large with respect to 
model mass and, secondly, the aerodynamic and magnetic cross-couplings 
may be severe, particularly at extreme model attitudes. More 
sophisticated control algorithms, perhaps including self-adaptation 
or self-optimisation with changing test conditions, thus seem 
necessary for effective operation of a LMSBS. No such algorithms yet 
exist. 

1.2.5 Miscellaneous 

No MSBS yet exhibits any viable hardware redundancy or utilises 
any failure tolerant features and, principally for this reason, all 
systems tend to be relatively unreliable in operation. 

Several methods of roll control exist but all are unsatisfactorily 
weak in torque capability. Thus, whereas the generation of adequate 
forces and torques in degrees of freedom other than roll appears 
practical at large scales, the generation of adequate roll torque 
does not. 

1.3 Contents of this thesis 

The underlying emphasis is towards the concept of a LMSBS though 
the majority of contents are unspecific in their application. At 
the commencement of the current phase of investigations, there 
appeared to be several specific and fundamental obstacles to be 
surmounted, or shown to be surmountable, before detail design of 
a LMSBS could commence. These include (not in any particular order 
of importance) : 


Generation of adequate roll torque. 


1 ) 

2) Provision of adequate overall system 

reliability and integrity. 

3) Design of very large capacity power amplifiers 

(loosely referred to in the context of 

MSBSs as power supplies) . 

4) Design of large multipurpose A.C. capable super- 

conducting E/M arrays. 

5) Development of large scale, precise and versatile 

position sensors. 

6) Development of advanced control algorithms including 

full magnetic decoupling and some self-adaptive 

features. 

Of course only a few of these topics could be addressed within 
the scope of a study of this kind. 3) and 4) above are quite 
beyond the capability and experience of this University. It had 
been felt initially that 6) above would represent a worthwhile and 
significant early advance but this view was quickly modified. 

Other problems, such as l)-5) quoted, were recognised as representing 
potential blocks to the development of a LMSBS, whereas advanced 
control algorithms need not be made available until quite late in 
the construction of such a system. Since it seems certain that, 
for reasons of reliability and versatility, the major part of the 
contarol system of a LMSBS would be implemented with some form of 
digital computer, merely realistic estimates of required computer 
power would be required in the early stages of design. 

Some attention was paid to 5) above but it was eventually 
concluded that the technological advances of certain sensing devices, 
such as solid state photo-detector arrays, currently occurring for 
application in other areas, may alleviate this problem in the 
relatively near term. 

The bulk of this thesis is therefore directed towards problems 
1) and 2) above, also the matter of extreme attitude testing in MSBSs 
which, although potentially very advantageous for a LMSBS does not 
represent a crucial block to its development. 
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2 . 


THE SOUTHAMPTON UNIVERSITY MAGNETIC SUSPENSION AND BALANCE 
SYSTEM (SUMSBS) 


2.1 Historical aspect 

Studies of magnetic suspension systems began at Southampton 
in 1959, construction of a system beginning in 1962, with the 
principal objective of performing dynamic testing in conventional 
wind tunnels without mechanical model supports. Work was under 
the direction of Dr. M. Judd, latterly and to date Dr. M.J, Goodyer. 
It was quickly shown that dynamic testing was indeed feasible 
although the high frequencies of oscillation required at small scale 
led to the development of the two-mass "tuned" model (11 ) . 

Extensive Magnus force testing was performed on ballistic-type 
models with reverse Magnus force, absent in corresponding sting 
supported tests, being discovered (12) . Several winged models 
(AGARD-G etc.) were tested in 6-component suspension and methods 
of rapid data acquisition were developed (13) . Studies of scaling 
of MSBSs (14) revealed difficulties that were also being identified 
by other groups of researchers, particularly that the power con- 
sumption of conventional copper conductor E/Ms in a MSBS scale as 
a high power of characteristic dimension. It became clear that 
systems designed for wind tunnels of realistic scale and Reynolds 
number capability, ambient temperature operation being mandatory at 
that time, would be quite out of the question, consuming megawatts 
of steady power. 

SUMSBS fell into temporary disuse but was re-commissioned in 1978 
and adapted to demonstrate operation with the University's 0.1m 
low speed cryogenic wind tunnel, the first known operation of a 
MSBS with such a tunnel (15, 16 ) . 

2.2 Modifications to SUMSBS to achieve current status 
2.2.1 E/M configuration 

As previously mentioned it seems certain that any future LMSBS 
will exhibit considerable E/M symmetry. It thus appeared desirable 
to reconstruct SUMSBS, previously of the unusual " I " configura- 
tion (11) to a more representative arrangement. Subsequently, it 
was realised that extensive E/M symmetry would be advantageous to 
suspension of models at extreme attitudes (Section 7). Further, 
spanwise magnet rolling moment generation systems (Sections 4, 5, 

6) require in general at least a quadrupole field disposed in the 
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cross-sectional plane of the test section. Thus, in order to 
demonstrate such roll control systems, suspension at high angles 
of attack and representative LMSBS configurations, the previous 
E/M layout of SUMSBS was augmented and modified by the addition of 
two new "vertical" E/Ms, identical to the six "vertical" and 
"lateral" E/Ms existing, and provision of twin symmetrical (fore-and- 
aft) "drag" E/Ms, replacing the single E/M previously used (Fig. 2.1; 
Ref. 11; Appendix 1) . 

The system appears to fall into the "+" subclass of Class I as 
defineid earlier but in order to accommodate spanwise magnet roll 
control, 6-component suspension takes place with the model's axes 
inclined at 45° to the horizontal, in the roll sense (Fig 2.2), 
wherein the system falls into the "X" subclass. 

2.2.2 Power supplies 

In common with the majority of contemporary MSBSs, SUMSBS utilised 
monopolar three-phase thyristor E/M power supplies. These were 
satisfactory insofar as they were relatively simple and had operated 
almost faultlessly for over a decade and a half. The low frequency 
switching inherent to the mode of operation does introduce con- 
sider<ible ripple in the output current and although this was not 
unduly troublesome in normal suspension there is an undesirable 
frequency response cutoff associated with the firing frequency 
(150Hz with U.K. mains supply). The monopolar operation presents 
obvious restrictions and the single quadrant capability of the 
supplies necessitated the installation of ballast resistors in series 
with each E/M to shorten the time constant of current decay. This 
also eases control system design, reducing the lag between E/M 
current and applied voltage, but would be quite unacceptable in a 
LMSBS due to the enormous power wastage that would result. 

Since some bipolar supplies were necessary for SUMSBS effectively 
to demonstrate suspension at high angles of attack (Section 7) it 
was decided to take the opportunity to solve many problems simul- 
taneously and incorporate a more sophisticated type of power supply, 
that is Pulse Width Modulated (PWM) transistor switching supplies, 
discarding the thyristor types completely. The supplies chosen 
(Appendix 1) were designed for use with industrial D.C. servomotors 
and feature the ability to store some regenerated energy from the 
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load. The switching frequency is 5kHz which permits more 

effective filtering of high frequency components of load current 

than was previously possible and renders the frequency response 

of the supplies near infinite as far as the control system is concerned. 

2.2.3 Position sensors 

The position sensing systems used in SUMSBS since its inception 
all fall into the category of simple analogue optical shadow move™ 
ment detectors and have never been capable of tracking a particularly 
wide range of model movement, perhaps * 15 degrees pitch/yaw, 
without mechanical realignment (11) . Other types of sensor exist, 
capable of tracking a wider range , notably the MIT electromagnetic 
position sensor (EPS) (10) , and other types are proposed with 
promise of still wider capabilities such as photodiode or CID array 
target trackers (5, 17 ) . SUMSBS requires 6-component sensors 
capable of tracking ideally 90° of pitch rotation with simple models. 
Despite long deliberations and much experimentation it was concluded 
that in the context of the current research none of the more advanced 
sensors could be incorporated into SUMSBS due to various time and/or 
physical space and/or financial constraints. Thus a new set of 
traditional sensors has been constructed with widerange pitch 
tracking capability effected by arranging that the appropriate sensors 
may be rotated about the model ' s pitch axis, conventional fixed 
sensors being used where extreme model attitudes are not required 
(Appendix 1) . The obvious should be stressed, that is that these 
sensors are quite inappropriate for use in a LMSBS and do not 
represent current thinking in this area. 

2.2.4 Control systems 

The simple all analogue control systems previously used with 
SUMSBS (11) were not capable of straightforward or realistic 
adaption to the requirements for suspension over wide ranges of 
attitude or for incorporation of more advanced control algorithms. 

The major part of the control system was thus abandoned and 
replaced by a mini-computer with local A/D and D/A conversion and 
limited analogue pre- and post-processing (Appendix 1) . 
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3. 


ROLLING MOMENT GENERATION SYSTEMS FOR MSBSs 


3 .1 Introduction, the need for roll control 

Historically, one of the key obstacles to the development of 
MSBSs for application to large scale wind tunnels was the 
development of a rolling moment generation system of adequate power. 
Nevertheless numerous small MSBSs have operated over many years 
without any active roll control system, indeed to date only two 
MSBSs have been demonstrated with such control (11, 18, 19) . It 
is pertinent to examine how this can be and why operation at large 
scales would be impractical without active roll control. 

Firstly, much useful testing could be, and has been, performed 
in small scale systems on bodies of revolution, notably sphere or 
missile drag measurements and Magnus force investigations , where 
active roll control was not required. Secondly, a considerable , 
percentage of the research effort expended with MSBSs has been in 
developing the systems themselves, particularly features such as 
position sensors, where suspension of models other than bodies of 
revolution is unnecessary. Where fixed roll attitude was required 
it frequently proved relatively easy to induce a preferred roll 
attitude, with slight positive stiffness about that attitude, by 
arranging the model's magnetic centre to lie above its centre of 
gravity . 

The overwhelming majority of testing in a LMSBS would certainly 
involve winged aircraft models. Small "tare" roll moments are 
likely to exist due to slight magnetic, aerodynamic or mass 
asyiranetries but presumably these could be countered by designing 
in a preferred roll attitude in the classical fashion. In general, 
however, whenever a winged model is tested with non-zero yaw angle 
(w.r.t, tunnel axes) significant roll moments will be generated, 
typically via the L^ aerodynamic derivative. In order to permit 
testing under these conditions a magnetic torque must be generated 
to oppose the net aerodynamic torque. Further, at high angles of 
attack the onset of stall is frequently asymmetric, producing steady 
roll torque of either sign or perhaps oscillatory torque , such as 
with asymmetric vortex shedding . Damping of roll oscillations is 
essential under all the above-mentioned conditions and is best 
introduced by appropriate manipulation of a method of positive roll 
control,, 
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It is difficult to estimate the magnitudes of torque that 
might be met in a future LMSBS because the aerodynamic character- 
istics of the models are of course unknown. However, Ref. 20 
includes the following data for an F-94A ; 

C (U.S. equivalent notation for L ) in the range 

Ln V 

-0.07 (high speed) to -0.05 (low speed) 

Span = 11.37m Wing area = 22.2m2 

Scaling to a 4-foot span model in an atmospheric 

pressure tunnel at M = 0.85 : 

L - 11 Nm /degree Yaw 

Data from Ref. 21 relating to an arbitrary design resembling a 
Hunter gives a somewhat larger result (43 Nm/deg) . It is clear 
that considerable roll torques are to be expected from this source. 

In the case of asymmetric stall we may take an arbitrary case of 
an aircraft with one wing completely stalled and the other com- 
pletely unstalled. A stalled C of approximately unity might be 
expected with the unstalled perhaps twice this value. Without 
reference to any specific geometry or test conditions it is 
immediately seen that in the (worst) case of an untapered uniformly 
loaded wing the resultant lift vector will act some one-sixth of 
the semi-span outboard of the aircraft centreline along the unstalled 
panel. 

There thus seems little possibility of realistic operation of 
a LMSBS without a powerful method of magnetic roll torque generation. 

3.2 Methods of magnetic roll torque generation 
3.2.1, Introduction 

Various methods of generating roll torque, falling into three 
main categories, have been proposed and evaluated over a number of 


years. These include: 

D.C. field systems 

1) "Bent" fuselage core (22) 

2) Shaped fuselage core (11, 23) 

3) Through wing magnetized wing 

cores (11, 18 ) 

4) Active model mounted coils (11 ) 

5) Passive model mounted coils (11) 

6) Spanwise magnets (24 ) 
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A.C. systems 

7) Planar conducting loop ( 19, 22 ) 

Aerodynamic systems 

8) Active aileron control ( 11) 

9) L (C ) control (11) 

Types 4, 5, 8, 9 are currently considered unsuitable for LMSBS 
application and types 1, 2, 3 have shown inadequate torque capa- 
bility in previous analyses. Type 7, the A.C. system, promises 
high torques but significant complication of overall MSBS design 
seems necessary to accommodate the strong A.C. fields required. 

An investigation of the two subclasses of type 6 is presented 
later. No studies of the current status of other methods have been 
undertaken so direct comparison between methods cannot be included. 

3.2.2. D.C. methods 

3. 2. 2.1 Theory of D.C. torque production 

At any point in a magnetized body, where magnetizations and 
applied fields are quasi-steady, the forces and torques acting on an 
incremental element of material can be written as follows s 

^ = (M . V ) H 6v ^ = (M X H) 6v 

If M represents permanent magnetization (polarization) there is 
no difficulty and H may be taken directly as the external or applied 
H field. Where M is composed wholly or partly of induced magnet- 
ization there is some ambiguity since M will be a function of applied 
Hi If the total force and torque on an isolated body are required 
then it is valid to consider the induced magnetization as permanent 
and the effective H field as that field that would exist if the body 
in question were removed. A body of magnetically soft material will 
require a low demagnetizing factor in at least one direction for the 
product (M X H) to be large. 

Resolving using the co-ordinate system of Fig 3.1 ; 

r_ = a^ + bj_ t c^ 

Fig. 3.1 

Axis system for magnetic 
calculations 
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( 3 . 1 ) 


Integrating roll torque over the volume of the core : 


Z {6t + (r X 6p ) } 

V X - — - X 


L =/(MH -MH) + (MH -t-MH +M H)b 
V y z z y X xz y yz z zz 


(MH +MH + MH)cdV 
X xy y yy z yz 


(3.2) 


3. 2. 2. 2 Review of some existing D.C. methods 

Shaped fuselage cores exploit the terms M H and M H . In 

X xy X xz 

normal suspension / M H dV and / M H dV represent the 
V X xz V X xy 

"lift" and "lateral" forces respectively and rolling moment can be 

generated essentially by redistributions of these forces over the 

volume of the core such that /(-M H c + M H b) dV is 

V X xy X xz 

non-zero. The core cannot be axisymmetric to satisfy this condition. 

Drawbacks of the method are that either (H ) or (H ) the 

xy z xz y 

vertical and lateral gradients of the conventional "sideforce" or 
"lift" fields respectively, must be of large magnitude with slender 
fuselages and that valuable fuselage core volume must normally be 
omitted to satisfy the condition of non-axisymmetry . 

Bent fuselage cores introduce transverse components of magnet- 
ization, M^, M^. Principal variations of the method proposed to 

date seek to exploit the terms M H or M H by arranging suitable 

yz z y 

distributions of transverse magnetization and applied fields to make 
the relevant integrals non- zero, for example the E coil method (22) . 
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Large transverse magnetization components will be required to 
generate high torques but this can only be achieved in general by 
severe camber of the core magnetization or by use of more complex 
multiple reverse camber core geometries which in turn necessitate 
rather awkward field distributions. 

The previous wing mounted magnet scheme used at Southampton 
University introduces asymmetrically disposed vertical components 
of magnetization , which, with a suitable applied field distribution , 
produces roll torque via the term above. Historically, the 

system was chosen to suit an existing E/M arrangement . The self 
demagnetizing factors of vertically magnetized wing cores are very 
high in all conventional aircraft geometries, making iron cores 
difficult to magnetize and permanent magnet cores susceptible to 
demagnetization . 


3 .2. 2. 3 General comments 

Moment generation via integration of fuselage forces (shaped cores) 
seems intuitively unattractive for conventional geometries due to the 
rather short moment arms available within a slender fuselage. 

However, there are at least three cases where this method may at 
least contribute usefully to an overall torque figure. These are : 

(1) Wide body transport. 

A wide body transport will exhibit simultaneously a high ratio 
of fuselage volume to first moment of wing area and a conservative 
test attitude/speed envelope. It may thus be argued that if 
sufficient forces and moments (apart from roll torques) can be 
generated with slender aircraft geometries then some of the 
fuselage core volume may be sacrificed, in the wide body case , 
for roll torque production. 

(2) P-15 style fighter. 

Several modern fighters are now exhibiting partial bifurcation 
of the rear fuselage to accommodate well separated twin engines. 

Such peculiarities of aerodynamic design create a non-axisymmetric 
fuselage of rather high moment of volume to volume ratio . 

Suitably exploited, therefore, such designs may yield useful roll 
torque without sacrifice of any core volume . 

(3) Underwing stores, engine pods, etc. 

These create a similar effect to that described under (2) above . 
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However, production of torque by exploitation of these character- 
istics could not be regarded as generally applicable for obvious 
reasons. 

Torque generation via transverse components of magnetization 
seems appropriate but attempts to introduce transverse magnetization 
in soft iron fuselage cores are hampered by the usual requirement 
for powerful axial magnetization. Permanent magnet fuselage cores 
would need to be carefully configured to avoid strong mutual de- 
magnetizing effects between ad j acent orthogonal magnetization 
components. 

It is logical to concentrate the spanwise magnetizations in the 
wings of the model , though parts of the fuselage may be used, since 
the self-demagnetizing factors of cores of the general proportions 
of conventional wings, magnetized along their long axis, are very 
low. 

Examination of equations 3.1 indicates that cross couplings 
between degrees of freedom are likely if transverse components of 
magnetization are introduced. Cross couplings exist in all 
practical MSBS designs but it has been standard practice to arrange 
for each degree of freedom to be controllable by relatively simple 
spatial distributions of applied field with cross couplings well 
suppressed at the datum model position. It is not thought that 
elimination of all significant cross couplings is a realistic goal 
for LMSBS design since a wide variation in model geometries, 
including magnetization distribution, is expected and a modestly 
sophisticated control system should be capable of handling quite 
severe cross couplings. There is considerable advantage in E/M 
array design, however , if primary forces and moments can be 
generated by simple field distributions. 

Simple unidirectional transverse magnetization cannot be ruled 
out simply on the basis of its severe cross coupling into pitch, 
nor because it destroys the mirror symmetry of the model, but it 
will exhibit a null torque line unless vertical fields can be 
generated relatively independent ly in the areas of axial (fuselage) 
and transverse (wing) magnetizations. 
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The null torque line exists in uniform 
applied field gradients. 


Symmetrically disposed transverse magnetization components 
will generate roll torque with application of a vertical field 
with transverse gradient as Fig. 3.3. 



Fig. 3.3 

Spanwise magnet 
rolling moment 
generation. 
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There are two subclasses of this system which will be dealt 
with separately and are referred to here as the Spanwise Permanent 
Magnet (SPM) and the Spanwise Iron Magnet (SIM) schemes. 

3.2.3. A.C. Methods 

The planar conducting loop method proposed by Stephens and 
since successfully demonstrated at M.I.T. is an ingenious approach 
to the problem and does have the advantage of minimal coupling 
between roll and other degrees of freedom. However, the strong 
A.C. fields required would be difficult to generate in a large 
scale system and may not be possible at all with current technology 
superconducting E/Ms. The predictions of roll torque capability 
appear promising but it is felt that the practical difficulties of 
operating an A.C. system at large scales would lead LMSBS designers 
to choose a D.C. system, if one were available with sufficient 
torque capability, in preference to the existing A.C. method . 
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4 


SPANWISE PERMANENT MAGNETS (SPMs) 


4.1 Introduction 

If transverse magnetization components are generated by 
installation of permanent magnet cores, predominantly in the 
model's wings, then roll torque is directly generated by appli- 
cation of the appropriate through-wing fields as Pig, 3.3. 

In all good permanent magnet materials at low applied fields 
the polarization (J) is practically independent of applied field 
strengths and may also be considered unifom over the volume of 
the core in the case of rare-earth cobalt (ReCo) materials, this 
class of permanent magnet exhibiting by far the highest known 
values of energy product and coercive force (25). A serious 
disadvantage of all known materials suitable for SPM application 
is poor mechanical properties, particularly high brittleness, 
that would preclude the use of permanent magnet cores as major 
load bearing members in the model's structure. Only some fraction 
of the available wing volume could thus be utilised for the SPM 
cores. 


4.2 Theory 

4.2.1 Ideal applied fields and simple spanwise magnetization 

Where the transverse magnetizations are entirely in the y 
direction (spanwise) then equations 3.1 reduce to: 


F 

X 

= / J 

y 

V 

H dV 

xy 
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y 

= / J 

y 

V 

H dV 

yy 
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z 

= / J 
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H dV 

yz 

L 

= / J 

y 

V 

H + J H b - J 

z y yz y 

M 

= / J 

y 

V 

H c ” J H a dV 
xy y yz 

N 

= /-J 

y 

V 

H + J H a - J 

X y yy y 


H c dV 

yy 


H b dV 
xy 


- (4.1) 


If it 
over the 

/ 

V 


is assumed that the applied field gradients are constant 
volume of the core then since : 

dV = O (by symmetry about the xz plane) 
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equations 4.1 collapse to: 


F =F = F = M = O 
X y z 

L = / J (H + H b) dV 

^ y 2 yz 

N = /-J (H + H b) dV 

^ y X xy 

further , if it is assumed that applied fields exhibiting large 
components, (axial, magnetizing, lift) will be symmetrical about 
the xz plane (Appendix 5), then : 

/ J H dV O 

y X 

V 

Thus we have: 

F = F = F = M= 0 
X y z 

L = / J (H + H b) dV - (4.2) 

y z yz 

V 

N =i f~J H b dV 
^ y xy 

It is immediately seen that, since is the principal applied 
field gradient for sideforce generation via conventional axial 
(fuselage) magnetization, there exists one fundamental cross coupling 
due to the presence of the spanwise magnetizations, that is coupling 
into yaw from a sideforce demand. It is not thought that this 
presents any serious difficulty since the magnitudes of aerodynamic 
forces and moments in the lateral plane are generally considerably 
smaller than in the vertical, whereas the magnetic force and moment 
capability due to the axial magnetization will often be approximate ly 
equal in' these two planes. 

Examining a SPM wing core with ideal applied fields we have (see 
Fig. 4.1) : 

From Eqns . 4.2 : 
s 

L = / J A g(b) (H + H b) db - (4.3) 

b=-s y o z yz ' ' 

Now H may be written as H b where H is the value of H at the 
z z^ — z z 

t s t 

tip. Assuming symmetry about the xz plane : 

L = 2/ JAg(b)H — db 

^ J y O z S 

O t 
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If J may be considered constant ; 

y 


L 



/ 


2A^g(b)b db 


(4.4) 



The integral is the first moment of volume of the core and depends 
only on its geometry. J is material dependent but does not vary 
dramatically amongst the best permanent magnet materials. The 
maximtut usable value of H depends on material properties and the 

internal demagnetizing influences of the core and can vary 
considerably between different materials and core geometries. 

If the cross sectional area of the core is constant with b then 
g(b) =1 and we have; 

L = 2H J A s - (4.5) 

z^ y o 

t 
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4,2.2 Theory ■ General applied fields and simple spanwise 
magnetization 

In realistic MSBS configurations , uniform applied fields or 
field gradients cannot be achieved . Some treatment of cross 
couplings arising due to the interaction between non-uniform applied 
fields and the transverse magnetizations is important and can be 
made using Eqns, 4.1. The general approach taken here is to break 
down the integrals of 4.1 into x, y and z forms and express all 
integral arguments as odd (O) or even (E) functions of x, y and z. 
The characte ri Stic geometrical symmetry (Fig. 4.1) immediately 
yields, writing an odd function of y as etc. ; 


/ O dy = O 

y — 

y (wing) 

Where the wing has a symmetrical section, similarly; 

/ O dz = O 
z (wing) 


(4.6) 


(4.7) 


Little error in fact occurs in using Eqn. 4.7 in most cases since 
the wing thickness is typically small with respect to other dimensions 
Applied fields may be represented by their characteristic symmetry 
in model axes, although this must be recognized as a considerable 
idealization, valid exactly with the model in the usual datum location 
but only approximately so otherwise. The characteristic symmetry 
of the applied through-wing roll field is as follows : 


H 

z 



- 26 - 



which may be represented by the matrix equations 



b = E o. E ; c = E E O and a = O E E 
xyz xyz xyz 


equations 4,1 become: 
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o 


N = - / O / E / O dz dy dx (J H term) 

X V z y X 

X y z 

+ / O / E / O dz dy dx (J H a term) 

X •' y •' z ^ y yy 
X y z 

- f O /E /Odzdydx(JHb term) 

x y z yxy 

X y z 


Thus equations 4.1 collapse to: 


= O 


= O 


F =F =F =M = N = 0 
xyz 

L =/JH + JHb-JHcdV 
„ y z y yz y yy 


(4.10) 


Of course c is small (Fig 4.1), thus equations 4.10 reduce 
further to: 


F =f =F = M = N = 0 
X y z 

L“/JH + JHbdV 
^ y z y yz 


- (4.11) 


It is immediately seen that there exist no primary couplings due to 
the transverse magnetizations with typical applied roll fields. 

There may, however, exist couplings due to the applied roll field 
with the conventional axial (fuselage) magnetization or due to 
other applied fields and the transverse (wing) magnetizations. 
Couplings between the applied roll field and conventional axial 
magnetization will be small with all realistic geometries, since the 
roll field exhibits a null along the x axis (Fig, 4.2), and will be 
zero with axisymmetric fuselage cores. 

The full calculations as shown above will not be carried through 
here for all cases but the important results are as follows; 

Table 4.1 Characteristic symmetries of applied fields 
(Appendix 5) 


Field classification Primary field or 

field gradient 
component 


Field component charac- 
teristic symmetry matrix 
(as Eqn. 4.8) 


Magnetizing 


E 

O 

O 


E 

O 

E 


E 

E 

O 


Drag 
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E 
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E 


E 

E 

0 
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Lift 


H 


xz 
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Yaw 


H 
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E E 0 

O O O 

O E E 

' O E o' 

E O O 

E E E 
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O E E 

, O 0 O , 

' O O E ' 

E E E 

1 E O O 


From equations 4.1 etc. the following conditions can be identified 
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These criteria can be related to the symmetry matrices in 
Table 4.1 either directly or by using; 
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The conditions above become: 
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The origin of many of these cross couplings is easily identified 
as core asymmetry in the yz plane. If this condition is abandoned 
and a fully symmetric core considered then further couplings dis- 
appear, leaving: 
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f O with magnetizing field 
F / O with yaw field 

y 

F^ f O with pitch field 

L =0 with all fields 

M O with lift field 

N f O with sideforce field 

4.2.3 Theory. General spanwise magnetization 



Fig. 4.3 Sweptback SPM core 


In the majority of cases of interest the model's wings will be 
sweptback as shown in Fig. 4.3 above. To preserve low demagne- 
tizing factors for the SPM cores the direction of magnetization 
must be sweptback also, perhaps to lie along the axis of centroids 
of core cross-sectional area. 

Full analysis of this case becomes complex since both and J 
components exist. However, if the SPM cores are assumed slender, 
whence applied fields may be considered uniform over any core 
cross-section, the fundamental effects of the x-wise magnetization 
components can be identified: 
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J = - J A g(b) Sin C 

X o 


J = J A g(b) Cos C 
y o ^ 



y 

o 






J 

z 

0 





From Eqns . 3.1; 
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/ J H + 
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X xz 
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(J H + J H ) c 
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y z 
V ^ 

+ b H 

yz 

c H ) 
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+ 

J (b H - c H ) 

X xz xy 

dV 

The J 

X 

terms may be 

broken down as before. 

using Eqns. 4.9; 
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But since c is typically small, we may write; 

L=;/j(h + bH)+J b(H ) dV > (4.12) 

y y X yz X xy 

The J term is identified as the sum of the moments of the "lift” 
X 

force (F^) components over the wing volume, the terms being as in 
Eqn. 4.11. 


4.3 Theoretical performance 

Direct evaluation of the integrals in Eqns. 4.1 etc. is not 
generally possible since field and field gradient components can 
seldom be represented as simple functions of x, y and z. Where 
the MSBS E/Ms are air cored, the applied field at any point within 
the model core may be calculated analytically or numerically 
(Appendix 2), permitting evaluation of the total forces and torques 
acting on the model, where the magnetization distribution is known, 
effectively by numerical integration of Eqns. 3.1 over the core 
volume. The magnetization distribution may be assumed known in 
the case of ReCo materials below the onset of irreversible de- 
magnetization, then typically being constant and uniform, otherwise 
tending to be difficult to compute. 

A key requirement for LMSBS application is the identification of 
the absolute maximum torque capability of a particular core, which 
in the SPM case will be determined by the onset of irreversible 
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demagnetization of the core. 

The total field within the core is the sura of the applied (E/M) 
field and the internal self-deraagnetizing field of the core. The 
self-demagnetizing field of cores of arbitrary shape and known 
magnetization distribution may be calculated numerically and certain 
idealised forms yield to direct analysis. Where the total field 
acts parallel to the local direction of magnetization the onset of 
irreversible demagnetization may be predicted directly from the 
materials BH (or JH) curve, ReCo materials generally exhibiting 
near straight-line curves in the 2nd quadrant (Fig. 4.4). 



Unfortunately, the total field will seldom act as described 
above, and whereas some data exists concerning the demagnetization 
of single crystals of certain materials with total fields not 
parallel to the direction of magnetization (for instance see Refs. 
26, 27, 28), it is not thought that this can be applied to bulk 
materials, demagnetization then being powerfully affected by 
processes principally dependent on the material's microstructure, 
such as domain wall movement (28, 29, 30). No relevant data 
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concerning bulk ReCo materials has been discovered (see Ref, 

31 for treatment of particular ferrites), rendering predictions 

of maximum torque capability of general SPM systems impossible at 

present. However, a crude estimate of the order of magnitude of 

available torque may be made by the following argument. 

In slender cores magnetized along their long axis, the internal 

self-demagnetizing field is small over the majority of the volume 

and may be neglected here. The effects of the magnetization 

reversal at the wing root (Fig. 4.1) and of the detail geometry of 

the cores is beyond the scope of this treatment. The magnetic 

field strengths internal to the cores may thus be regarded as pre- 

dominantly due to the applied field. The conventional parameters 

characterising the performance of permanent magnet materials, such 

as H , H . , J , BH etc. (25), are insufficient for identification 
c cx r max 

of incipient irreversible demagnetization. Several authors resort 
to use of the parameter H , representing the demagnetizing field 

JC 

strength that reduces the intrinsic polarization to 90% of the 
remanent value (25, 32, 33) , illustrated in Fig. 4.5. 



Fig, 4.5 Definition of 



Values of for relevant materials are sparsely quoted, 
reflecting the arbitrariness of the parameter's definition, some 
available data being reproduced in Table 4.2 below; 


Table 4.2 Some properties of bulk ReCo materials 


Material 

J 

r 

(T) 

H 

c 

«k 

(kA/m>— 

H . 
cx 

Reference 

Sintered 

0.85 

650 

900 

1200 

25 

(1977) 

Sm COg 







"Commercial" 

0.84 

560 

488 

800 

32 

(1971) 

Sm Cc)^ 
5 







Sm Co^ based 
5 

0.90 

700 

1000 

1200 

33 

(1976) 

"RECOMA 20" 








Current material developments are expected to result in 
improvements in H and J with H tending to comfortably exceed 

C IT K 

At this stage it would not appear unreasonable to assume a 

usable value of H of some 800 kA/m (U H -IT) and J of 0.85 T. 

K O JC ]T 

If this value of can be taken to represent the peak total field 
strength that can be applied to the material before the onset of 
irreversible demagnetization, with the working polarization taken 
as 90% of the value of quoted immediately above, then some 
progress may be made. 

The applied field is a sum of many contributory fields ("lift", 
"drag", etc.), but it should be noted that the two fields (neg- 
lecting roll) likely to involve the most powerful fields within 
the test section need not cause the exposure of the wing core 
to those high field strengths (Fig. 4.6). 

If it can thus be assumed that the principal component of applied 
field within the wing cores is the applied "roll" field then H 

^t 

in Eqns. 4.4, 4.5, etc. becomes H . Taking the values assumed 
above, we have from Eqn. 4.4; 

2.4 X 10® ^ 

L =! f A g(b) b db - (4.13) 

s o 

o 

and from Eqn. 4.5 for rectangular slab cores ; 

L - 1.2 X 10®A^s - (4.14) 
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Fig. 4.6 Components of lift and drag fields 
at the wing core 

Eqn. 4.14 corresponds to 600,000 Nm/m^ of core. For the represen- 
tative SPM core studied in Section 4.3.1 below, this would indicate 
a maximvim capability of approximately 500 Nm. It may be, however, 
that applied fields of sufficient strength to realise this capability 
cannot be generated with existing technology E/Ms. In fact, taking 
the roll torque calibration calculated below (Section 4.3.1), maximum 
torque would occur at around 16,000 A/cm^ in the specified E/Ms. 

This is well beyond current E/M technology (see Section 5.8.2) . 
Present apparent LMSBS roll torque requirements (5) may be 
approachable with more carefully configured E/M arrays . 

4.3.1 Performance of representative SPM configurations 

Since accurate analytic predictions of performance (below the 
onset of demagnetization) are not practical, numerical examination 
of representative SPM configurations is appropriate , the computer 
program FORCE (Appendix 2) being suitable for such studies. Some 
preliminary results are presented here. 
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The geometry of E/M array and model is chosen so as to be 
directly comparable to the baseline geometry used in the studies 
of the SIM system (Section 5) and readers should refer to Section 
5.4 for explanation of the evolution of these geometries. The 
wing core is chosen to occupy one half of the available cross 
section (Pig. 4,7), the aerodynamic and other loads being assumed 
as principally carried by the surrounding structure. A 
quadrupole field elegantly provides the required through-wing 
field, in turn being supplied by the appropriate quartet of E/Ms 
from the baseline SIM array (Figs. 4.8, 5.4), though the entire 
SIM array is used for examination of 360° roll capability. 

Confirmation of some of the predicted couplings between axial 
or spanwise magnetization and the applied roll or other fields is 
straightforward. Conventionally configured but arbitrarily 
dimensioned E/Ms are specified to provide lift, pitch, etc. fields 
(Figs. 4.9, 4.10). It should be noted that the overall E/M array 
is reaiListically proportioned although the axial and main E/Ms do 
(unintentionally) spatially conflict. 

The magnetic forces and torques generated on the model and their 
variation with simple displacements or rotations of the model are 
shown in Figs. 4.11 - 4.19. Confirmation of the predicted magnetic 
couplings due to the spanwise magnetizations can now be made. 

These are (from Section 4.2.2): 


Force/torque 

Applied field 

Figure (zero displacements) 

F 

X 

Magnetizing 

4.19 

F 

Yaw 

4.18 

y 

F 

Pitch 

4.17 

z 



L 

None 


M 

Lift 

4.14 

N 

Side force 

4.16 


The wing core's symmetry about the yz plane is conveniently 
destroyed by a displacement in the x direction. Consulting the 
relevant Figures therefore, the extra predicted couplings for a 
core lacking symmetry about the yz plane may be confirmed also: 


- 37 


SIM cores (Section 
5) this size 


V v_ Non-magnetic 2 
envelope (7.2cm ) 

Magnetic core(7.2cm^) 


Detail view on A showing direct comparison 
to baseline SIM geometry (Fig. 5. 5) 


Equivalent airfoil 


Fig, 4.7 Model magnetic configuration for SPM calculations 













Applied fields defined by the following relationships of currents ; 








Fig. 4.11 Couplings for wing in roll field 


E/Ms as Fig. 4.8, 4.10 
Model as Fig. 4.7 
Koll field defined by : 

+I =-i =+i =_i = 1,459,240 A (1000 A/cm^) 

11 12 13 14 



All forc:e and torque components not shown are exactly or 
effectively zero. 

Fig. 4.11a Fig. 4 . 11b 
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Pig. 4.13 Couplings for fuselage in roll field 


E/Ms as Fig. 4.8,4.10 
Model as Fig. 4.7 
Roll field defined by : 


4-T T sr—T 

11 12 13 14 


1,459,240 A 


"l " ^10 = ° 


Model rotations and torques as Pig. 4.11 

All force and torque components not shown are exactly or 
effectively zero. 


Fig. 4.13a 


Fig. 4.13b 


T — I — I — I — I — r— I — r~ 

0 0.1 0.2 0.3 0.4 

X (m) 
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Fig. 4.14 Couplings for wing in lift field 









Fig. 4.15 Couplings for wing in drag field 

E/Ms as Fig. 4.9,4.10 
Model as Fig, 4.7 
Drag field defined by ; 

+I„ ="I,_ = 2,500,000 A (1000 A/cm^) 

y xo 

I, - and I„ - I. . = 0 
1 8 11 14 


Model rotations and torques as Fig. 4.11 

All force and torque components not shown are exactly or 
effectively zero. 


Fig, 4,15a Fig. 4.15b 



X (m) y (m) 
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Fig^ 4.16 Couplings for wing in sideforce field 


E/Ms as Fig. 4.9,4.10 
Model as Pig. 4.7 
Sideforce field defined by : 

-1, = +1. = +1. = = 1,600,000 h (1000 A/cm^) 

^468 



Model rotations and torques as Fig. 4.11 

All force and torque components not shown are exactly or 
effectively zero. 

Fig, 4.16a Pig. 4.16b 






Fig. 4.17 Couplings for wing in pitch field 


E/Ms as Fig. 4.9,4.10 
Model as Fig. 4.7 
Pitch field defined by : 

-I, =+I, = +1^ = 1,600,000 A (1000 A/cm^) 

i 3 5 / 



Model rotations and torques as Fig. 4.11 

All force and torque components not specified are exactly or 
effectively zero. 

M is non-zero but less than 0.003 Nm in Figs. 4,17a“f 


Fig. 4.17a Fig. 4.17b 



0 O.X 0.2 0.3 0.4 ■ 0 0.1 0.2 0.3 0.4 


X (m) y (m) 
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Fig. 4,19 Couplings for wing in magnetizing field 


E/Ms as Pig. 4.9,4.10 

Model as Fig. 4.7 

Magnetizing field defined by ; 

+I„ =+I-- = 2,500,000 A (1000 A/can^) 
9 XO 


Ig and 1^^ - = 0 


Model rotations and torques as Fig. 4,11 

All force and torque components not shown are exactly or 
effectively zero. 

Fig. 4.19a Pig. 4.19b 




56 - 








Force/torque 

Applied field 

Figure (positive x 
displacement) 

F 

X 

Drag 

4.15a 

F 

Sideforce 

4.16a 

y 



F 

Lift 

4.14a 

z 



L 

None 

-- 

M 

Pitch 

4.17a 

N 

Yaw 

4.18a 


It should be noted that the two couplings involving generation 
of pitching moments (from applied lift or pitch fields) are very 
weak, due to the short moment arms of wing core elements about the 
y' axis in the chosen configuration. 

The full information presented in Figs. 4.11 - 4.19 indicates 
that further couplings will occur when the model is displaced or 
rotated from its datum position, for instance a lift field is seen 
to generate pitch and yaw torques when the wing core is rolled. 

The principal result of these will be the manifestation of magnetic 
stiffnesses (positive or negative) in various model degrees of 
freedom when the classical fields considered are applied. These 
effects will somehow disturb the magnetic stiffnesses and dampings 
generated by conventional control algorithms and may require action 
to counter adverse effects. 

However, study of the Figures reveals that these extra couplings 
occur in a relatively straightforward manner and few are of large 
magnitude. Whilst further treatment of this topic cannot be pre- 
sented here, it is felt that modest adjustments to and expansions 
of conventional control algorithms will restore full required 
performance . 

Since no attempt has been made to ensure representative relation- 
ships between the strengths of the classical applied fields, 
further interpretation of the data discussed above should be made 
with extreme caution. 

4.4 Discussion 

The performance of SPM systems is critically dependent on the 
precise magnetic properties of the core material, particularly the 
intrinsic polarization and useful coercivity (see Section 4.3). 
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There seems good reason to expect considerable improvements in 
the values of these parameters exhibited by the best available 
materials in the near future. Fig. 4.20 shows how the highest 
reported values of the classical parameter of magnetic 

perfonnance, have risen consistently over the past decades, the 
performance of ReCo materials still being well below the theoretical 
limits at least in the respect of intrinsic coercivity (28, 34) . 

The torque capabilities calculated heretofore are thus expected to 
rise as time passes within the limits of the assumptions made. 

Further, it is widely accepted that as the temperature of ReCo 
materials is lowered (below 300K, say) their magnetic properties 
tend to improve steadily. The parameters of interest here, that 
is intrinsic polarization and coercive force, are typically found 
to increase by a few per cent and many tens of per cent respec- 
tively, during cooling from room temperature to around 80K (28, 30,35, 
36, 37 and Fig. 4.21) . The force and torque capability of a ReCo 
cored MSBS model would thus be significantly greater at the typical 
lowest operating temperatures of a cryogenic wind tunnel than at 
"normal" temperatures. 

4.5 E'/M configurations for SPM roll control 

The E/M configuration chosen so far is a fully symmetric quadrupole 
(Fig. 4.8). There is no particular reason to regard this as being 
in any way optimum. 

From Fig. 4.11 it is seen that roll torque capability falls 
rapidly with increasing roll displacement from the selected datum 
orientation. This would be a serious limitation in practice, being 
due to the inability of the chosen array to develop ® 

y' axis displaced by 45° in the roll sense. This problem may be 
completely overcome by including extra E/Ms in the roll array. 

These may perhaps be shared with other force/torque generation duties 
but will be regarded here as entirely specialised for simplicity. 

The complete 8 E/M array used for the SIM calculations (the 
SPM 4 E/M array being, of course, derived from this) provides roll 
torque at any roll angle, the variation of torque capability with 
angle being shown in Fig. 4.22. 

If 360° roll angle capability is deemed unnecessary, then the 4 
E/M array may be adapted to increase the angular capability, certain 
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' (Sm,Pr)Co 


AlNiCoW o PtCo 


Fig. 4.20 Increase 
in highest reported 
values of energy 
product with time 


After Ref. 3 5, 


AlNiCo ®y 


MnBi/ •BSD 




Sr-ferrite 




^ FeCoV 
AlNiCo 


1900 1920 1940 1960 1980 




After Ref.28 


SmCo 16.84 At% Wt. 

^ Sm. (Ref. 37) 

SmCo_ mean of several 
® (Ref. 30,36) 






E/M geometry as 
Pig. 5.4 

E/Ms a as Fig. 4.8 

All currents are 1,459,240 A 
(1000 A/cm^) 

Current directions chosen 
to create positive torque 
at all angles 



z 


Contribution of ; 



Pig. 4.22 Variation of roll torque with roll displacement for 


complete SIM E/M array (SPM cores) of Fig. 5.4 
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array to enhance angular capabilit; 






Fig. 4.24a Variation of roll torque capability with roll angle for 
modified quadruplet geometries. 




E/Ms as Fig. 4.23 Model as Fig. 4.7 



Fig. 4.24b Variation of roll torque capability with roll angle for 
modified quadruplet geometries. 




5. SPANWISE IRON MAGNETS (SIMs) 


5.1 Introduction 

Synunetrically disposed transverse magnetization components can 
relatively easily be induced in a magnetically soft wing core by 
application of a symmetrical field as shown below: 



Rolling moment will be generated with application of through 
wing fields as in Fig 5.2., Section 3.2.2., and Section 4. 
Rotation of the induced spanwise magnetization vectors by the 
through wing fields will be inhibited by the high demagnetizing 
factors in the through wing direction. Calculation of the 
performance of SIM systems is not straightforward, the induced 
magnetizations not being directly analytic. 
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H (positive) 



Fig. 5.2 SIM through wing field 

Thus, whilst the SPM equations may be approximately applied if 
the induced magnetizations remain predominantly spanwise and the 
value of polarization is replaced by some mean value of induced 
magnetization, that value of magnetization and the detail perfor- 
mance of the SIM system cannot be adequately predicted without 
further analysis. 

However, many magnetically soft materials exhibit saturation 
inductions of approximately 2 Tesla and high permeabilities to around 
that value. Operation with higher levels of spanwise magnetization 
than those achievable with SPMs (approximately 1 Tesla) is thus 
possible, leading at least to higher torques per unit through wing 
field. Further, it is assumed (Section 5.4) that magnetically soft 
materials with good mechanical properties are available, allowing the 
entire wing volume to be utilised if required, and since no irrever- 
sible demagnetization effects exist in the SIM case it is to be 
expected that the absolute maximum available torque for any config- 
uration may be greater than with SPMs. 

Since soft magnetic materials generally behave entirely isotropically, 
available material data is sufficient to permit, in principle. 
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calculation of the performance of the SIM system at all applied 
field levels. 

5.2 Solution of iron-air-conductor magnetostatic systems 

This class of problem has attracted considerable interest over 
recent years with the application of superconductors to various 
fields such as nuclear research and with the drive to improve per- 
formance and efficiency of electrical machines. 

General magnetostatic problems involve the solution of Poisson's 
eqviation or, outside current carrying conductors, Laplace's equation. 
These equations may be solved in principle by a number of methods 
including (after 38) 

1) Analogue 2) Graphical 3) Analytic 4) Numerical 

Predictions of roll torque from SIMs require treatment of 
saturation effects, leading to non-linear solutions, and are inherently 
three diim’nsional. In most practical cases methods 1) and 2) above 
are not able to handle the non- linear problem, indeed graphical 
methods are generally restricted to two dimensions. Analytic methods 
are available for non-linear 3D problems but only for highly restricted 
geometries of conductor and iron. Application to general problems is 
currently quite impractical. Again following (38) existing numerical 
mcjthods in this field may be divided into four principal categories: 

1) Finite difference 2) Images 

3) Integral equations 4) Variational formulations 

although other classes of solution do exist, such as the Monte-Carlo 
method (39) . 

It would appear that image methods are inapplicable to non-linear 
problems within the bounds of the present formulations. Method 4) 
above may be considered an energy method, somewhat analogous to virtual 
v,'ork methods in structural problems, whereas 1) above tac)cles Poisson's 
or Laplace's equation directly but both generally require that the 
complete volume of the problem is meshed with a computation grid. 

The characteristic geometry of the SIM roll system is very "open" 

(Fig. 5.3), that is a small iron region separated from the conductors 
by large air gaps. Methods 1) and 4) above would thus require meshing 
of considerable volumes of air, at least enclosing all the conductors, 
leading to large computation times. Integral equation methods 
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(3 above) need only require meshing of iron regions and thus appear 
appropriate here . 

5 ■ 3 Program GFUN 

This program, developed at the Rutherford Laboratory, Didcot, 

Oxon, U..K. since 1970 by Trowbridge et al, is an example of the use 
of integral equation methods. GFUN has been applied to a wide 
range of problems over many years (38, 40, 41, 42, 43) and has 
exhibited consistently good accuracy, consequently gaining a high 
international reputation. A disadvantage inherent to the solution 
procedures employed is a somewhat awkward representation of field 
within iron regions, as the vector sum of the fields from external 
currents and induced magnetizations, expressed as constant within 
suitably shaped elements (such as tetrahedra) . In high permeability 
regions the two contributions to the field may be nearly equal and 
opposite, leading to rather poor resolution of the internal field. 

GFUN in fact handles the simultaneous integral equations as a single 
matrix equation, solution of the latter generating eigen solutions for 
the internal field. Clearly, prediction of the detail geometry of 

the internal field is affected by the choice of element distributions. 
The lack of precise information concerning the internal field of the 
SIM cores is not thought to be especially critical. It must be noted 
that the nature of the solution procedure implies that the foregoing 
adverse comments do not apply outside the iron region, however para- 
doxical that may seem, and need not apply to the resolution of forces 
and torques provided methods of field integration over control volumes 
external to the iron are chosen. 

Access to GFUN was granted by the British Science Research Council* 
under Grant No. GR/B/3691.5. Modifications to the program were carried 
out by Simkin to permit the full symmetry existing in many of the 
required SIM cases to be exploited (reducing computation time) and to 
provide a torque integration option. The torque integration scheme 
is conceptually similar to the well established methods for force cal- 
culations but the fact that previous users have apparently not required 
torque information must be seen to represent a major possible source 
of systematic error in data included hereafter, indeed some diffi- 
culties were experienced before a consistent integration scheme could 

* Now the Science and Engineering Research Council. 
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be identified. To counter this uncertainty, various low field 
torque measurements have been made and limited cross checking under- 
taken (Section 5.7). Saturation of iron regions, which could not be 
achieved in corroborative experiments, should not directly influence 
the reliability of torque predictions due to the nature of the 
methods used, 

5.4 Baseline GFUN geometry and configuration 

A baseline geometry and configuration is required from which the 
effects of variations of various parameters, such as wing aspect ratio, 
may be examined. At the time of commencement of this study there 
existed no clear specific choice of geometry of either model or E/Ms, 
necessitating evolution of a baseline geometry on the following, 
somewhat arbitrary, basis. 

The favoured scale for studies of LMSBSs has been a test section of 
approximately 8 ft x 8 ft cross section (such as NTF) . Calculations 
are therefore made at this scale, but the scaling of results to 
different tunnel sizes is quite straightforward (Appendix 3) . The 
test section is assumed square with no corner fillets. Clearance is 
allowed around the aerodynamic cross section for structure, plenum 
chamber etc., and is chosen to be one foot (5). A similar allowance 
is made for the thermal insulation and structure surrounding each E/M 
(5) . The most uncertain characteristics of the E/Ms are the winding 
shape and maximum usable current density. It appears (5) that manu- 
facturers prefer circular windings where possible due to reduction of 
difficulties associated with conductor stressing. An idealized 8 E/M 
configuration has been chosen (Fig. 5.4). The problem of optimizing 
the E/M array must be dealt with separately, being heavily influenced 
by particular requirements for forces and moments in other degrees of 

freedom. Maximum usable overall current density (J) for super- 

2 

conducting E/Ms varies from order 1500 A/cm for cryostable conductors 
2 

to order 15000 A/cm for adiabatically stable conductors, within the 
limits of present technology. Doubts exist (5) as to whether adia- 
batically stable conductors could be applied to the 8 ft scale case so 
generous winding cross sections are allowed in the E/Ms, permitting 
partial saturation of the wing cores at suitably low values of J. 

High J data is included for the purpose of identifying trends at high 
field levels. The current levels in the E/Ms are represented by 
and as defined by Figure 5.4. 
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The geometry of the baseline SIM wing core is simply chosen as a 
10% thick, 10; 1 aspect ratio rectangular hexahedron (rectangular "slab" 
core) spanning slightly under one half (49.2%) of the test section 
width (Fig. 5.5). 

The choice of core material (or its magnetic characteristics) 
presents some difficulty at this stage due to the necessity in a final 
design to achieve an optimum compromise between many magnetic, electrical 
and mechanical properties. The SIM wing cores are visualised as 
occupying nearly all the aerodynamic volume of the wings with only those 
detail fittings and features potentially subject to changes during test 
programs being added from non-magnetic material. It is recognised that 
this view may require qualification in the light of detailed results and 
experience in model design, but it is nonetheless clear that the core 
material requires good mechanical properties such as high yield point, 
low brittleness and high Young's Modulus, perhaps also at low tempera- 
tures. High electrical resistivity may be preferred in order to 
suppress eddy current flow in the cores but is not essential to funda- 
mental operation of the system. Low coercive force is necessary to 
avoid torque calibration being dependent on the recent past magnetic 
history of the cores, although if cores are operated far into saturation 
this feature becomes of lesser importance. Ultra high permeability is 
not necessary since in any airgap dominated system, iron regions of 
medium and higher permeabilities tend to behave as if infinitely 
permeable. The key magnetic parameter is undoubtedly saturation 
induction. Studies of LMSBSs (5) have indicated that valuable economies 
can be made in E/M size, hence cost, by utilising fuselage core materials 
with the highest available saturation induction. It is logical that 
similar criteria should apply to the SIM case since operation of the 
core beyond the saturation point is anticipated. 

One class of materials that appears promising is the cobalt iron 
alloys, classically represented by the 50:50 Iron-Cobalt alloy 
"Permendur". These materials have not found widespread industrial 
application due to their relatively high cost but this is unlikely to 
be a problem in LMSBS applications. Indeed the cost of a fuselage 
core for an 8ft tunnel has been estimated at U.S. ^ 50CX) (1981 ^) , 

(5). This is small compared to a typical manufacturing cost of a 
large wind tunnel model. 

Magnetic and mechanical properties of material of this type are 
dependent on the precise alloying constituents, heat treatment and 
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preparation of samples but typical properties for some commercially 
available materials are shown below. 


Material 
(Trade Name) 

Perraendur 

Vanadium Permendur 

Vacoflux 50 

Source 

Western Electric Co. Ltd. 

Vacuumschme Ize 
GmbH 

Reference 

44 

44 

45 

Density 
(kg/m ) 

8300 

8300 

8150 

Initial y 

800 

800 

lOOO 

Saturation Induction 
(T) 

2.45 

2.4 

2.35 

Static coercivity 
(A/m) 

159 

159 

110 

Resistivity 

0.7 X 10“^ 

— fi 

0.26 X 10 

0.35 X lo"® 

Youngs Modulus 
GPa 

- 

* 

230 

Yield strength 
MPa 

- 

- 

400 


It should be noted that the saturation inductions, in the range 
2.35 - 2.45T (exceeding 2.5T in laboratory specimens) are substantially 
higher than the corresponding value for high purity iron of some 
2.158T (room temperature. Ref. 44). 

The peaks in the permeability versus magnetizing field strength 
curves for the above materials (Fig. 5.6) present some problems to 
GFUN. The permeability of each iron element in GFUN is assumed 
constant throughout the element and is updated at each iteration. 

Sharp rises in the permeability of elements with relatively small 
increases in the magnetizing field acting on those elements (and vice 
versa) occurring from one iteration to the next may cause local 
oscillations of the iron's interior field and consequently slow 
convergence. These potential problems may be avoided in early work 
by choosing a permeability versus magnetizing field strength charac- 
teristic that falls monotonically . This has been done by arbitrarily 

fitting intermediate points between the initial constant permeability 
line and the terminal saturation boundary, thus establishing the base- 
line characteristic shown in Fig. 5.7. A slightly conservative value 
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Magnetizing field strength (A/m) 



The effects of variations of some 


of 2.3T is chosen for M 

sat 

material properties on the torque capability of the SIM system is 
studied later . 

The baseline wing element distribution is chosen along well 
established principles . Tetrahedral elements have proven to be the 
most reliable choice for GPUN, the total number of elements is the 
largest that can be handled by the convenient batch version of the 
program and the spatial distribution of the elements is chosen such 
that the elements are relatively numerous in the region of strongest 
anticipated magnetization . This non-uniform element distribution 
also yields superior convergence . The control volume for use with 

the torque integration schemes is dimensioned such that its surfaces 
lie close to the model core, thus yielding best accuracy. 



Fig. 5.7 GFUN baseline BH curve 
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5.5 Simple theoretical approximations 


5.5.1 Torque on ellipsoids 

Some insight into the behaviour of the SIM system may perhaps 
be gained by study of some magnetostatic torque producing system 
that is analytic. Unfortunately, as previously mentioned, few such 
systems exist and few of these yield to straightforward analysis. 

One suitable example is torque on magnetically soft uniform ellipsoids 
in uniform applied fields. The general equations applying are as 
follows 



H 

X. 
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N M 
X X 
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Combining (5.1) and (5,3) 


M = u H . . 

X X (-r— —— ) 

e 1 + Y N 

''■m X 


(5.4a) 


M = y H , 

y O Y„ ( 


TB 


e '1 + X N ' 
y 


(5.4b) 


The torque on the ellipsoid is as follows; 


T = (M H - M H ) V 
X y y X 

e e 


(5.5) 


Following Fig. 5.8 we have from (5.1) 


H. Cos6 = H Cos(b - N M COS0 
1 e x 


(5.6a) 


H. Sin0 = H Sin<i) ~ N M Sin0 
1 e y 


(5.6b) 


Combining (5.6a) and (5.6b) 


1 = 


Cos(j) 


H. + N M 

1 X 

I / 


SirxJ) 


H. + N M 

X 


(5.7) 


This equation gives from any given (j> and M. Below satura- 
tion M, H and X ate functions of each other (Eqn. 5.3). Solution 
of (5.7) thus involves an iterative process using a material's 
particular BH curve. When H. is found, X is known from the BH 
curve and torque can be found from (5.4) and (5.5). This procedure 

is somewhat inconvenient. Where y is large (permeable material 

m 

below saturation) we have H^-0 and Eqns. 5.1 collapse to; 


M 

X 


li H Cosd) 
o e ^ 


N 

X 


- (5.8a) 


M U H Sind) 

y - o e ^ 


- (5.8b) 


and from (5.5); 


u H Cos(J) Sin(j) V "(5.9) 
o e 


Directly, for maximum torque at any given applied field strength 
(j) = 45° and: 

T a applied field ^ (at any given 4>) “(5.10) 
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Where the material is far into saturation M = M , = constant, 

sat 

and Y - M ^ 
sat 


'^o X 


Prom (5.6) : 


H Cosd) 
_e 

~Cos0 


N M 
X sat 


Sir»(|) 

SinG 


N M ^ 
y sat 


Re-ax'ranging : 

H Sin ((j)-e) = ’sCN “ N ) M Sin 20 

6 V X 




(5.11) 


Giving : 




(j) = Sin 


And from (5.5) 


(N - N ) M ^ Sin 20 
y X sat 


2 y H 
^o e 


+ 0 


1 

V 


M H Sin (ij)-0) and using (5.11): 


(5.12) 


T 

V 


M (N - N ) Sin 20 
sat y X 

2 


(5.13) 


(5.13) shows that 


T 

VJ Max 


is independent of and has the value: 


V/Max 


M (N - N ) 

sat y X 

2 


(5.14) 


where 0 = 45 . Now (N - N ) can never exceed unity thus: 

y 3C 


M 


V / Max, 


'Max 


sat 

2 y. 


(5.15) 


Equation (5.12) becomes: 


. -M (N “ N ) M 


Sin 


X sat 1 + 45 


(5.16) 


/Max 


2 y H 
o e 


It should be noted that (5.16) and (5.12) break down where H is 

e 


Sin$ 


small. 

From (5.6) 

again: 
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TORQUE 

/VOLUME 


(Nm/m X t . 4 

10 -^ 


. PROLATE 


OBLATEN 


Long axis = a (x wise) 
Short axis = c (z wise) 
Applied field in 
xz plane only 
Torque about y axis 


Ellipsoid data 0.4 
from Ref. 47 


b/a=0.lN 


M ^ = 2.35 T 
sat 


Pig. 5.11a Torque per unit volume (maximum) for ellipsoids. 
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At this point, from 


(5.18) 


Fig. 5.9 shows the variation of the direction of induced 
magnetization, after Eqn. 5.12, for a variety of cores. It is seen 
that with a powerful applied field the induced magnetization lies 
nearly parallel to that field. 

Fig. 5.10 shows torque per unit volume for a particular core 
(Eqn. 5.13) and clearly illustrates the existence of a finite and 
approachable absolute maximum torque. E’or any given geometry and 
applied field direction it is seen that the torque per unit applied 
field is far from constant. In fact, at high values of ()> with 
powerful applied fields, increases in applied field strength may 
apparently result in reductions in torque. 

The variations of the maximum torque and torque per unit volume 
with ellipsoid geometry are of interest and are shown in Fig. 5.11, 
using data from Ref. 47. The significant feature here is that 
relatively slender ellipsoids are capable of generating high torques, 
perhaps implying that slender wing cores (SIMs) will remain relatively 
capable torque producers. 


5.5.2 Linear SIM theory 


H , = (k^ Cos2(|) + Cos4(f) + 

™ r r 

, = -1^ Sin2^ + Sin4(J) + —) 

p r "‘r 


(5.19a) 


where k and 1 are functions of r. For an identical quadrupole 
n n 

r r 

displaced by 45° about the x axis (giving the through-wing field 
quadrupole) : 

H , 


(kj^ Sin2<j) + k^ Sin4(() + — — ■ — ) 
r r 


(5.19b) 


(1^ Cos2<j) + 1^ Cos4(|) + 

^p r r 


•) 


If a linear magnetization characteristic is assumed for the wing 
core (permeability constant or »1) and a suitable radius chosen, we 
may write: 


M , pH, 

y ^ o y^ 

D 

y' 


M , pH, 
z' o z* 
=; R 


(5.20) 


where D , and D , are effective demagnetizing factors in the y' and 

y z 

z' directions respectively. If H^^, is assumed to be fairly constant 
along any particular radial we have: 


H , TI Tr 
z ' - yz ' 

r 


(5.21) 


Equation 4.4 becomes: 


L =! 4M ,H , 
y ' yz ' 


A g(b)b db + M , term 
o z 


b = o 

For constant geometry and neglecting the M^, term: 


L “ M ,H , - (5.22) 

y yz 

Applying a magnetizing field and a through-wing field simultaneously: 




Im(ki Cos2<J) + — ) + I^(k, Sin2<}) + — ) 


R 


r ' 1 


R 


Tl (1 Cos2(|) + ) - I (1 Sin24) + — -) 

L r ij^ m 1 r 

Neglecting terms above 2(p: 

L <i >5(1^ - I^)Sin4(j) + l^I^Cos4(|) 


r m 


(5.23) 
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Fig. 5. 13 Variation of roll 
torque with (|), 

Simplified linear theory. 






This equation is illustrated in Fig. 5.13. Where (^ = 0, it is 
seen from Eqn. 5.23 that: 

L « 1 1 “ (5.24) 

* r m 

Thus, in a fully symmetric configuration, the maximum torque for 
given total available ampere-turns will occur at; 

I = I (<t> = 0) 

5.6 GFUN results 

5.6.1 Baseline geometry - Effect of variation in material properties, 
low and intermediate applied fields 

Where the permeability of the core is constant (linear solution, 

typical of low applied field levels) the effects of the magnetizing 

(h^) and through-wing (J^) fields are independent. Since, under these 

conditions, the spanwise magnetization is everywhere proportional to 

and the through-wing field proportional to J^, the roll torque is 

expected to vary as the product J J . This is confirmed in Fig. 

m r 

5.14 - 5.16, being as predicted by Eqns. 5.9 and 5.24. 

The variation of torque with core permeability, permeability held 
constant within each solution, is of interest. Fig. 5.17 showing com- 
paratively low sensitivity to permeability variations at high values 
of permeability, the core then behaving as if infinitely permeable. 

Using the baseline BH curve it is clear that there exists some 
maximum level of applied field commensurate with the whole core lying 
in the initial constant permeability region of the BH curve. At 
higher applied field levels the permeability of certain strongly 
magnetized voliimes of the core will progressively fall, eventually 
the bulk of the core settling onto the terminal (saturated) region of 
the BH curve. During this process the magnetization is no longer 

proportional to J , rather reaching some limiting value. It might 
m 

therefore be expected that torque becomes proportional to alone. 

However at high values of through-wing field the induced magnetizations 

may no longer be predominantly spanwise, rather turning to lie more 

nearly parallel with the direction of local (applied) field (see 

Fig. 5.9). Torque may therefore reach some limiting value , or 

continue to rise (or fall, see Fig. 5.10) with rising applied fields 

as some function of J and J (Fig. 5.18). 

r m 
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Fig. 5.18 Schematic diagram showing possible 
terminal torque characteristics 

Figs. 5.19, 5.20 show the breakaway from the linear permeability 
cuirve onto a near straight line characteristic at moderate applied 
field levels. Data is included for a close approximation to the 
representative vanadium permendur BH curve in Fig. 5.6, showing a 
very weak influence of precise core magnetic properties. 

5.6.2 Variations of baseline geometry. Effect of wing AR, 
t/c and taper 

GPUN solves linear (constant permeability) cases directly, 
without recourse to iterative procedures, hence relatively economi- 
cally. The variations of the initial (constant permeability) 
torque capability of the system with various geometrical parameters 
as defined in Fig. 5.21 are shown in Figs. 5.22 - 5.24, the solid 
symbols representing the baseline core. 

These figures require some e3q>lanation. It is clear that the 
torque for a given applied field is far from being a constant per 
unit core volume and in fact does not obey any simple relation to 
geometry (such as first moment of volume about the x axis) . It 
is believed that this effect is due to the fact that the effective 
spanwise demagnetizing factors are predominantly determined by the 
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Fig. 5,22 Torque versus chord. Baseline E/Ms, Core as Fig. 5. 21. 
Linear solutions. 



Fig. 5.23 Torque versus thickness. Baseline E/Ms. Core as Fig. 5. 21 
Linear solutions . 
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Fig. 5.24 Torque versus taper. Baseline E/Ms. Core as Fig. 5. 21. 






slenderness of the core, being lower for more slender cases. Thus 
the removal of wing volume by reduction of chord or thickness lowers 
the effective spanwise demagnetizing factor, hence increasing the 
level of spanwise magnetization of the core for any given applied 
field level. This may be partly justified by inspection of Figs. 

5.25 - 5.26, showing the increase in the peak element spanwise 
magnetization with increasing slenderness. 

Tapering the core under the criteria chosen does not affect the 
core volume but displaces volume from the most effective regions 
near the tips to less effective regions near the roots. The slender- 
ness of the tip portions of the core is increased and since these lie 
in the regions of strongest applied fields the peak value of element 
magnetization tends to be increased as shown in Fig. 5.27. 

5.6.3 Effect of the presence of fuselage and axial magnetizing field 

It is to be expected that the presence of an unsaturated iron 
fuselage should act to increase the mean level of spanwise magneti- 
zation in the wing cores since it provides an easy flux path at the 
wing root (Fig. 5.28). A saturated fuselage may or may not act 
similarly. With the fuselage geometry as defined in Fig. 5.29, the 
wing core span remains unchanged, the root now being enveloped by the 
fuselage core. The table below shows some results for this geometry, 
confirming the expectation of augmentation of torque with fuselage 
present. The effect appears weak. 


Fuselage Dimensions (cms) 

Torque (Nm) 

Absent 

109.78 

50 X 5.2 X 5.2 

110.22 

50 X 7.4 X 7-4 

113.53 

100 X 7.4 X 7.4 

114.40 


If a soft iron fuselage core is used it will generally require an 
axial magnetizing field. The effect of an axial field on the wing 
cores will tend to be to rotate the spanwise magnetization vectors in 
the plane of the wing, in the sense of sweepback or sweepforward. 
Where the core permeability is constant, the magnitudes of the span- 
wise components will be unaffected. A reduction in the spanwise 
components, hence torque, will be expected where the wing core 
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permeability is falling with rising magnetization. With the fuse- 
lage dimensions fixed at 50 x 7.4 x 7.4 cms (Fig. 5.29) and a near 
uniform axial field of representative strength applied by a large 
Helmholtz E/M pair of otherwise arbitrary dimensions this expectation 
is confirmed by the results shown in Fig. 5.30. 

The magnitude of the axial magnetization components induced in the 
wing core by the axial field will depend to a large extent on the 
axial slenderness of the core, which will in turn be most strongly 
affected by the core's thickness to chord ratio, relatively thick, 
narrow chord cores being least powerfully affected. Very thin cores, 
such as the F-16 case studied later, may thus be seriously affected 
although no relevant data currently exists. 

5.6.4 Effect of sweepback 

Sweepback of the wing cores as defined in Fig. 5.31 rotates the 
easy axis of wing core magnetization away from the spanwise direction, 
but does not affect the core volume nor the position of the centroid 
of volume. The induced magnetization vectors with a purely spanwise 
magnetizing field are thus expected to themselves be sweptback, at 
some angle probably less than the geometrical sweep angle. This 
would, in general, lead to lower torque per unit applied field. 

However, at low sweep angles it would appear that the increasing 
slenderness of the core, caused by the chosen geometry, dominates, 
leading to slightly augmented torques at modest sweep angles, as shown 
in Fig. 5.32. 

When an axial magnetizing field is applied, components of that field 
act along the easy axis of magnetization of the sweptback core, in- 
creasing or decreasing the spanwise magnetization components (hence 
torque) depending on the field polarity. For a sweep angle of 30 
degrees Figs. 5.33, 5.34 show the effect to be significant. At 
relatively low spanwise fields it is seen that powerful axial fields 
of either polarity reduce the torque for particular spanwise and 
through-wing field levels. 

5.6.5 Behaviour at high levels of roll torque generating field 

At high applied roll field levels the core becomes saturated over 
most of its volume. The induced magnetization components then behave 
as vectors of constant strength but variable direction, as predicted 


97 



















Fig. 5.33 Effect of axial field on sweptback core. Baseline E/Ms. 
Core as Fig. 5.31. Sweep = 30°. 5 point BH curve. 



Fig. 5.34 Effect of axial field on sweptback core. Baseline E/Ms, 
As Fig. 5.33. 
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for ellipsoids in Section 5.5.1. It is not immediately obvious 
whether, under these conditions, increases in applied field should 
lead to increases in torque generated and, indeed, whether or not 
some absolute maximum torque capability will exist in this case for 
each particular core geometry. 

The studies in Section 5.5.1 suggest that an absolute maximum 
should exist and that this maximum can be realised with finite 
applied fields of appropriate sense and direction. 

Various high applied field cases have been computed with GFUN 
(Figs. 5.35, 5.36) with somewhat inconclusive results, a clear limit 
or limiting trend not being identified. It should be noted here 
that the torque integration schemes in GFUN are expected to become 
progressively less reliable at increasing applied field levels. 

This is due to the difficulty of accurately resolving on the surface 
of the external control volumes (Section 5.3, Fig. 5.5 and (43) ) 
the "model" field (due to the core's induced magnetization and now of 
essentially fixed magnitude) from the total field, which becomes 
mostly due to the applied field from the E/Ms. 

The peak torque levels achieved in Figs. 5. 35, 5.36 are, however, 
at least one order of magnitude in excess of apparent existing LMSBS 
requirements (5) . 

5.6.6 Pseudo F-16 wing core performance 

It is understood that a representative aircraft type for use in 
LMSBS design studies has previously been chosen to be the F-16 fighter. 
This type presents a considerable challenge to the SIM roll scheme 
since its wing thickness (hence volume) is very low, the blockage 
effects of the fuselage would necessitate choice of model wingspan 
considcirably below the 50% of test section width used heretofore and 
the extreme taper both in chord and thickness leaves relatively little 
volume in the magnetically most effective regions of the wing (the 
tips) . 

Each wing panel has been crudely represented with GFUN as a hexa- 
hedral slab, uniformly tapered in both thickness and chord, with 
approxi.mately the same span, total core volume and moment of volume 
about the chordwise centroids of volume as a typical F-16 model (Fig. 
5.37). The element distribution, particularly the element aspect 
ratio, within the core is at the limits of what is generally 
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Fig, 5.35 Torque versus constant. High applied field levels. 


Baseline E/Ms. Baseline core. 5 point BH curve. 



Fig. 5.36 Torque versus Baseline E/Ms and core, 5 point BH curve. 
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acceptable with GFUN, hence the results must be regarded as subject 
to increased uncertainty (see Section 5.7) but again comfortably 
exceed apparent LMSBS requirements (5) at moderate applied field 
levels (Figs. 5.38 - 5.41) . Ref. 5 was, in fact, based upon a 
requirement for 140 Nm roll torque at the scale computed herein, this 
in turn being based upon a requirement for testing at atmospheric 
stagnation pressure, Mach. 0.85. 

5.6.7 Effect of variations of E/M geometry 

Non-circular E/Ms, for instance those in Fig. 5.42, exhibit 
improved packing around the wind tunnel test section compared to the 
circular baseline E/Ms. The field distribution in the region of the 
model will be altered also, although the effect tends to be slight 
due to the relative remoteness of the E/Ms. The performance of 
different E/M configurations may be approximately normalized by an 
appropriate measure of the field in the region of the model, but 
calculated with the model absent. This has been done by computing 
the mean field level along the y axis, taken over the span of the 
model. It is seen in the tables below that the effect of E/M 
geometry on the generated torque under these conditions is small. 


E/M geometry 

Normalized J 

Torque ( y = 1000) 

r 

Fig. 5.4 

1000.0 

120.51 

Fig. 5.42 

894.97 

119.78 


For reasons other than production of rolling moments, the classical 
array of 8 E/Ms distributed in the yz plane may not be preferred. A 
16 E/M system, shown in Fig. 5.43 has therefore been computed as an 
example with E/M performance normalized as above: 


E/M geometry 

Normalized J 

r,m 

Torque ( y = 1000) 

Fig. 5.43 

869.65 

112.82 


It is clear that the performance of the SIM system is not strongly 
affected by the detail geometry of the E/Ms, hence permitting con- 
sideration of alternative E/M geometries and configurations with the 
existing GFUN results being approximate ly applicable provided E/M 
performance is normalized by the model's near field using Fig. 5.44. 
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Fiq. 5.38 Torque versus J . Baseline E/Ms. "F-16" core geometry. 
— E — r.m 



Fig. 5.39 Torque versus J^, constant. Baseline E/Ms. "F-16" core 
geometry . 
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= 894.714 A/cm 



equivalent 


J in non-circular E/Ms 
corrected to approximate 
baseline field in region 
of model 


Fig. 5.42 C omparison of 
baseline and equivalent 
non-circular E/M geometry, 

All dimensions in metres. 
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Fig. 5.44a Near field of GFUN baseline magnetizing field 
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Fig, 5,44b Near field of GPUN baseline through-wing field 



5.7 Verification of GFUN data 


5.7.1 Alternative torque computations 

No direct alternative computations of the performance of any 
representative SIM system have been attempted since this would 
necessarily require access to an alternative computer program of 
comparable power and sophistication to GFUN, preferably solving 
the magnetostatic system by an entirely different method. Such 
programs probably do exist. 

However, as mentioned earlier, the reliability of GFUN's pre- 
diction of magnetization of iron regions has been considered good 
for many years and there seems no particular reason to regard the 
computed wing core magnetizations as being subject to any more 
than the usual levels of uncertainty (43) . Were this the only 
source of uncertainty, the computed torques could be regarded as 
likely to be accurate, within the limits of geometrical and other 
assumptions, to plus or minus very few percentage points at all 
computed applied field levels for the baseline and similar geometries, 
with somewhat greater uncertainty in certain cases, such as the F-16 
core, where the iron element distribution was sparser than desirable 
(see Section 5.7.3). 

The main potential source of uncertainty must be regarded as the 
torque integration schemes themselves, since these have been specially 
developed for the SIM computations and hence not subjected to such 
extensive testing and verification as the rest of GFUN. 

It is possible with GFUN to arbitrarily fix the induced magneti- 
zation in iron elements, thus effectively converting them to permanent 
magnet material. If this is done in such a way as to approximately 
preserve the typical spanwise magnetization in the SIM cores (Fig. 
5.45) then direct and representative verification of the torque 
integration schemes is possible. Alternative calculations have been 
made using the computer program FORCE (Appendix 2) which calculates 
forces and torques by elementary numerical integration of the relevant 
vector products of applied field and core magnetization over the 
volume of the core. This method differs fundamentally from the 
methods used in GFUN. Results for the geometry of Pig. 5.45 are as 
follows: 
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Fig. 5.45 Approximate spanwise magnetization distri bution (along y 
axis) for comparison of GFUN and FORCE torque integration 
schemes. 
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Torque integration scheme 

Predicted Roll Torque 


(Nm) 

GFUN by Maxwell field stress 
integration over surfaces of 
external control volume. 

140 X 40 X 20 cms. overall 

113.94 

FORCE by vector product 
integration over volume of 

116.82 

core 



The discrepancy of approximately 2.5% is considered acceptable 
though leaving scope for future improvement. 

Certain comparisons have been made of predicted torques over 
differing GFUN control volumes with typical results as follows: 


GFUN control volume overall dimen- 
sions (cm) 

Relative torque to 
standard volume 

140 X 40 X 20 (assuming symmetry 

1.0 

in the yz plane) 


140 X 40 X 20 (no symmetry) 

0.999993 

180 X 70 X 50 (yz symmetry) 

0.9856 


These results are similarly considered acceptable. 

As mentioned in Section 5.6.5 it is thought that the accuracy of 
the GFUN torque integration schemes will fall with rising applied 
field level when the SIM cores are well into saturation. No direct 
high field computations were made with permanent magnet cores 
specified and this was a serious omission but has been partially 
rectified by more detailed analysis and computation using existing 
data. Specifically, since the publication of Ref. 48, the induced 
magnetization distribution of two high field GFUN cases has been 
used as input data to FORCE, each GFUN element now representing a 
uniformly magnetized tetrahedra of permanent magnet material. 

The vector product integration procedures carried out by FORCE 
are not likely to yield accurate estimates of the "true" torque pro~ 
duced in these particular cases, indeed such procedures would other- 
wise be used in GFUN, being relatively straightforward. The reasons 
for this cannot be fully explained herein, but are principally 
founded on the fact that GFUN’s prediction of local external field 
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is markedly superior to that of local., internal field, hence magneti- 
zations . In these cases the predicted internal magnetization dis- 
tribution will be somewhat erroneous, but the field external to the 
core should be well resolved (see (41, 42) ). Pig. 5.46 shows that 
discrepancies do exist between GFUN's torque predictions and the 
FORCE estimates from GFUN magnetization data, and that the magnitude 
of these discrepancies rises with rising applied field level. In 
the principal region of interest (J <10,000 A/cm^) it is considered 
that the agreement is sufficiently close, bearing in mind the comments 
above, to manifest significant extra confidence in the GFUN torque 
integration scheme, although further verification would clearly be 
advantageous . 

5.7.2 Experimental measurement of torque with low applied fields 
and correlation with GFUN predictions 

Experimental verifications of the bulk of GFUN's predictions are 
not possible without an array of powerful (high field) E/Ms. Such 
an array was not available. However, the 8 main E/Ms from SUMSBS 
became available during the reconstruction of that system, enabling 
some representative low field (approximately constant permeability) 
torque measurements to be made. These E/Ms are not, unfortunately, 
axisymmetric but it is believed that the SIM system is not particu- 
larly sensitive to detail E/M geometry, rather to the mean applied 
field levels in the region of the SIM cores (see Section 5.6.7). 

The geometry and characteristics of the experimental SIM cores is 
shown in Fig. 5.47, the E/M layout in Pig. 5.48 and the ejqierimental 
layout in Fig. 5.49. Torque could be applied to the wing cores via 
a fine thread, pulley, scale pan and weights (Fig. 5.49), the 
experimental procedure being to measure the equilibrium angular 
positions of the wing with known E/M currents and applied torque. 

There are, as predicted by Eqn. 5.23, generally two such positions 
in each octant, one stable, the other not. Pig. 5.50 shows a 
typical calibration curve for one octant, the general form of which 
is repeated for all combinations of currents. Taking the (measured) 
reference angle for zero torque at zero roll current as 91.6°, Eqn. 
5.23 may be used to generate a set of theoretical calibration curves 
with insertion of an appropriate scaling factor. This factor is 
chosen so as to zeroise the cumulative torque error over all measured 
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Fig. 5,47 Geometry and characteristics of SIM cores for experimental 


verification of GFUN predictions 







Magnetizing" E/Ms marked A 
Through-wing" field E/Ms marked B 



Field probe support 










Reference angle (degrees) ( ({) ) 


Pig. 5.50 Typical curve of roll torque versus 
angular position for 8 E/M experimental SIM 
system. 


25 A cable current in E/Ms 
corresponds to 10^000 A-turns 



data points, taking the measured angular positions as correct. 

Fig. 5.51 shows the comparison between these theoretical curves 
and all measured data. 

Certain irregularities in the experimental data are inevitable. 

The non-axisymmetry of the E/Ms is expected to result in slight 
loss of left-right symmetry of the calibration curves, taken about 
their nominal peaks, in Fig. 5.51. Some evidence of this is 
apparent but the effect appears weak. Friction and irregularities 
in the wing assembly support bearings, together with the usual 
uncertainties in measurements, introduce significant random and 
systematic errors, which are considered to be of comparable magnitude 
to the observed residual discrepancies in Fig. 5.51. It is thus 
concluded that this Figure effectively verifies the simplified 
theory of Section 5.5.2. 

The GFUN representation of the E/Ms was, of necessity, somewhat 
idealised (Fig. 5.52), and a small correction to E/M current density 
proved necessary to achieve correct predicted field strengths in the 
region of the model. Complete measurement of the field of the 
experimental E/M array and comparison with the GFUN representation 
proved to be too laborious, therefore the current density correction 
was calculated by measurement of the through-wing field components 
along a datum spanwise axis of the experimental SIM cores, with 
the cores absent. This axis corresponds to the reference axis of 
the E/Ms (p = 90°, see Fig. 5.51) rather than the spanwise axis of 
the SIM cores with the cores in their zero roll torque, zero roll 
current orientation (p = 91.6°). The comparison between experimental 
and predicted (corrected current density) fields is shown in Fig. 

5.53. Since the magnetizing field E/Ms are nominally identical 
to those creating the through-wing field, the data shown effectively 
verifies the complete E/M array, although detail discrepancies in 
the field distribution are still possible. GFUN's predicted field 
distribution (corrected) is shown more fully in Fig. 5.54. 

Comparison between selected performance curves from Fig, 5.51, 
two GFUN predictions and fitted linear theory (Section 5.5.2) is 
shown in Figs, 5.55 and 5.56. Only two GFUN points were computed 
due to the lengthening of the calculations caused by the complex E/M 
geometry. Taking into consideration the extensive computational 
idealisations and various possible sources of experimental error. 
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Fig. 5.53 Comparison of measured field from experimental 
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Fig. 5.54 GFUN field distribution for SIM magnetizing field 
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Fig. 5.56 Comparison between measured 
and predicted roll torque 
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the agreement is considered to be excellent. The apparent dis- 
crepancies may, in fact, be unrepresentatively small (see Section 
5.7.3), being less than the estimated experimental error alone. 

It is concluded that major systematic errors in the predictions 
of GFUN cit low or moderate field levels are unlikely. 

5 .7.3 An assessment of likely accuracies of GFUN results 

Considering previous results with GFUN (40, 41, 42, 43) regarding 
prediction of magnetization levels, the agreement between the GFUN 
torque integration scheme and an alternative method applied to a 
representative low field case (Section 5.7.1.) and the performance 
of GFUN in predicting torque for the experimental system (Section 
5.7.2.), it is thought that the error magnitude in any GFUN prediction 
for the baseline or similar geometries at low applied field levels 
should not exceed 10%. Typical error may be significantly less than 
this figure, perhaps 5%. The error will tend to be mostly systematic 
in nature, available data showing low random content. This implies 
that trends in performance should be reliably identified provided 
fairly consistent formulations of the problem (for instance element 
distributions) are used. This was in fact done. The effects of 
geometrical and other idealisations are not included in the above 
figures. It being anticipated that increased inaccuracy may occur 
with rising applied field levels and a predominantly saturated core, 
it would seem appropriate to increase the figure for peak anticipated 
error to perhaps 20% for intermediate applied field levels (arbi- 
trarily 2000 lOOOOA/cm^). At still higher applied field 

levels the predictions become progressively less relevant to immediate 
requirements (5) since the torques predicted are high and the peak 
E/M fields required to achieve those torques are outside the limits 
of existing technology (Section 5.8.2.). The high applied field 
results (J^ ^ > 10000 A/citi^) should therefore, perhaps, be regarded 
for the time being as somewhat speculative and requiring further 
verification, such as more detailed testing of the GFUN torque 
integration scheme at high field levels. 

Significant departures from the baseline core geometry (sweep, 
etc.) are achieved only by utilising iron element distributions that 
are undesirably sparse. A version of GFUN exists (using the sub- 
program GETM 400) that can deal with up to 400 independent iron 
elements, rather than the 100 in the standard batch program, which 
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could tackle most cases herein with nearly ideal element distri- 
butions. This facility did not become available during the period 
of this study and anyhow would require much greater CPU time for a 
solution, CPU time for solving the magnetizations being limited and 
varying approximately as the square of the number of iron elements. 
Thus, cases such as those involving core sweepback, the effect of 
axial fields and the F-16 cores, are subject to increased uncertainty. 

Estimates of the likely peak error magnitudes in the above cases 
are principally based on judgement but realistic estimates are 
thought to be; 


Case 

Peak antici- 
pated error. 

Comments 

Baseline 

low fields 

10% 


high fields 

20% 

Rising with applied field 

very high fields 
(J > 10000) 

50% 

Requires further veri- 
fication. Ditto. 

Sweepback < 20° 
low fields 

20% 

Standard integration volume. 
No usable symmetry. 

Sweepback > 20° 
low fields 

25% 

Increased integration 
volume necessary 

high fields (positive 
Axial) 

30% 


high fields (negative 
Axial) 

40% 

Poor convergence of 
solutions. 

Axial field with fuselage 
low fields 

20% 

Large integration volume . 

high fields 

30% 

No usable symmetry. 

F-16 

low fields 
high fields 

40% 

60% 

Undesirably sparse element 
distribution. 


The principal idealisation inherent in all cases herein is the 
representation of wing chordwise cross sections (airfoil sections) as 
rectangles, since it has not (Section 5.4) been thought that slab 
cores buried inside non-magnetic aerodynamic envelopes would be used 
in practice. In order that the existing GFUN predictions be appli- 
cable to MSBS cases with true airfoil sections some appraisal of the 
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effects of this idealisation is necessary. However, it is not 
obvious on what basis, apart from cross sectional area, GPUN data 
can be matched to true sections. In the computed ;'-16 cases the 
moment of volume about tho centroid of area was chosen arbitrarily 
to generate a GFUN representation. The thickness no chord ratio, 
defined conventionally, could also be used, leading to generally 
similar results with classical streamline 'wing sect ions. Typical 
comparison between a GFUN section and a true sectic a is shown in 
Fig. 5.5 7 : 


GFUN Slab core 



Fig. 5.57 Comparison of GFUN and true c ore cross 

sections 

Since the spanwise slenderness is scarcely affected, it is 
expected that the mean spanwise magnetizations sho.ild be near 
identical in both cases, though the detail distributions must be 
different. The chordwise and through wing slenderness will be 
affected by the change of section, though the pn^c'.se value of 
through wing slenderness is not thought to be critical at low or 
moderate fields. The (increased) chordwise slera ^rness with the 
true section will lecid to somewhat greater sensitl vity to axial 
fields, the amount by v/hich the effects of the a>,i al field are 
ampli, fif'd being perhaps in the range O - 50%. 

5.8 Disc ussion an d conclusions 

5.8. 1 F urther use of GFUN in SIM computations 

GFUN'fi creators could not have envisaged that i would see use 
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in this type of problem. Because of this# several detail features 
are inconvenient for SIM computations. Some of these, such as the 
lack of a proven torque integration scheme, have been partially 
rectified (by Simkin at Rutherford Laboratory) during the course of 
this study, but some remain. Examples are the inability to exploit 
all the symmetry existing in all problems, the lack of a torque 
integration scheme with the high element number version of GFUN, and 
the doubts concerning torque predictions at high applied field levels. 
It is believed that these difficulties can be overcome relatively 
easily by further improvements or expansions of the program code. 

If this were done it seems likely that GFUN could provide torque (and 
force) predictions to a much higher level of accuracy than those 
estimated for the results herein. Definition of E/M and model core 
requirements for specific performance demands (all degrees of freedom) 
for specific LMSBSs should then be possible, to adequate accuracy for 
preliminary and intermediate LMSBS design and cost studies. It is 
difficult to imagine GFUN or similar programs becoming sufficiently 
accurate for precise predictions of performance (say better than 1% 
accuracy) in all cases of interest, due principally to the geometrical 
idealisations required in the formulations of the problem, though this 
level of accuracy should be attainable for certain simple cases. 
Empirical calibration of practical systems would therefore appear 
mandatory. Universally high accuracy may, however, be considered 
unnecessary, for instance where significant over-capacity is incor- 
porated in LMSBS E/Ms. 

5.8.2 Application of the SIM system to LMSBSs 

Even taking account of the relatively large errors considered 
likely to exist in the computed results, it is clear that the apparent 
torque capability of the SIM system considerably exceeds those pre- 
dicted heretofore for other rolling moment generating systems. For 
instance, with the baseline geometry (5 point BH curve), the value of 
600,000 Nm/m^ of (permanent magnet) core predicted in Section 4.3, 
perhaps 300,000 Nm/m^ of wing volume, is exceeded at approximately 
4000 in all E/Ms. 

Subject to the provision of adequately powerful E/Ms the SIM scheme 
must be considered a viable contender for LMSBS application. The 
matter of E/M design requires further analysis, but data is included 
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(Fig. 5.58), showing the baseline magnetizing field at 1000 A/cm . 

The baseline through-wing field is similar but indexed by 45°. 

Examining this data it is difficult to imagine the specified E/Ms 

2 

being operated above, perhaps, 3000 A/cm , with existing technology 
superconductors, whence the peak field within their bores becomes 
approximately 6.6 T. That value of field may not be increased by 
geometrical adjustments, only by improvements in superconducting 
E/M technology. Whilst certain optimisations of E/M geometry 
remains possible, it is seen that a fundamental limit to the available 
torque from a particular SIM system exists. 

Wliere peak performance is required it is easily seen that the E/Ms 
should generally be located as close to the model as possible 
(Fig. 5.59). 

The matter of cross couplings between roll and other degrees of 
freedom due to the presence of the spanwise magnetizations, or the 
(applied) wing core magnetizing or through wing fields, has not been 
directly addressed here. It is felt that the analysis presented in 
Section 4 will remain approximately applicable though the presence of 
axial magnetizations and/or the lack of fore-and-aft symmetry in 
swept wing (such as the F-16) cases will result in relatively complex 
couplings, mostly second order in magnitude. 

This study has indicated that at least with respect to available 
roll torque, the SIM system is viable for application to LMSBSs. 
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SIV\i<i’WISE PEMANENT MAGNET ROLL CONTROL WITH SUMSBS 


6 . 

"A Poduction 

Eince the Spanwise Permanent Magnet roll control technique was 
viev/ed as potentially the most suitable method for use with SUMSBS 
in the foreseeable future, preliminary suspension and minimum cali- 
bration of a pilot model has been attempted in order that this view 
might be verified. 

6_.2 Special fe atures of SUMSBS for roll con trol 

SPM roll control requires that SUMSBS be used in its X con- 
figuration, achieved by rotating the model's datum axes by in 

the roll sense (Fig. 6.1). The E/M geometry thus generated doess 
not directly correspond to previously discussed EPM E/M arrays 
(Section 4) due to the fore-and-aft separation of the two quadrupole 
groups of main E/Ms, but the required through-wing fields are still 
genercited relatively effectively (Fig. 6.2). 

The pilot SPM model merely consisted of a standard axisyrnmetric 
model (Appendix 1) , crudely adapted to carry vaguely representative 
permanent magnet wing cores (Fig. 6.3 and Appendix 1). Al], magnetic 
cores were Alnico V. The awkward appearance of this model must be 
excused but it should be noted that the wing core cross sections and 
volumes are in approximately representative proportion to the fuselage 
core for an aircraft of relatively high wing volume, such as an A-10. 
The measured torque capability may, of course, be satisfact<.:rily 
scaled to differing wing sizes (Section 6.4) . 

The roll attitude sensing system, described more fully it) Appendix 
1, consisted of a low power laser directed onto a model-mounted 
mirror, the reflected beam falling onto a light spot position 
deteci:ing device. Some difficulties were expei'ienced v/ith stray 
reflections of the laser beam, but the sensing system is not neces- 
sarily seen as representing anything more than a short term solution. 

All-digital control was employed (Appendix 1) v;ith the roll 
stabilisation loop incorporating, purely for convenience, only one 
phase advance network (see 11, 59), instead of tiie two (in series) 
used in all other degi'ees of freedom. This corresponds to previous 
practice with SUMSBS. 

The magnetic couplings occurring due to the presence of t^he si^an- 
wise magnetizations were expected to require adjustments to the 
control algorithms to alleviate their effects, but accurate estimates 
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Fig. 6.1 Operation of SUMSBS in an X configviration (SPM roll control) 
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of coupling effects were not possible for two principal reasons; 

1) The iron cores of the main E/Ms preclude accurate 
representation of these E/Ms with FORCE (Appendix 2) 
and measurement of the actual field distributions in 
SUMSBS was deemed to be too lengthy. Therefore, the 
line conductor approximation to SUMSBS' configuration 
and geometry used in Section 7 (Figs. 7.5, 7.18) was 
employed for this study, with the current levels in each 
E/M adjusted so as to match the calculated field with 

the measured field at the origin of balance axes (Fig. 6.5). 
The field distributions around this origin are not expected 
to be particularly well represented. 

2) The polarization distribution of Alnico V cores will 
not be constant and uniform, as is assumed in all FORCE 
computations. The relatively low coercivity of this 
material typically results in weakening and splaying of 

the polarization vectors, particularly around the extremities 
of the core, partly due to self demagnetizing effects and 
partly due to local demagnetization by applied fields, as 
illustrated in Fig. 6.4. 



N S T3 S 


Idealized Practical 


Fig. 6.4 Idealised and practical polarization 
distributions for Alnico V cores . 

For FORCE calculations the core polarization was assumed constant 
and uniform at 1 Tesla. This figure is below the saturation polar- 
ization of Alnico V (=^1.25T, Ref. 25) but may represent an over- 
estimate of the effective polarization level. The forces and 
torques generated due to various classical demands (calculated) 
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are shown in Tables 6.1 and 6.2, with the E/M geometry and current 
directions clarified in Fig. 6.5. Pure roll torque or drag force 
components are generated by classical roll or drag demands respec- 
tively and there exist no couplings into roll torque or drag 
force from other demands (lift etc.). It is concluded that de- 
coupling of the controller is unnecessary for preliminary work 
since no coupling terms are of the same order of magnitude as the 
relevant primary terms (Table 6.2), with only one coupling 
appearing as significant, that being the expected coupling into 
Yaw from a Sideforce demand. Further, generous stability margins 
are available . 

Table 6.1 Required current directions for generation of classical 
demands with X configuration (SUMSBS) 

E/Ms as Fig. 6.5. 

Current directions shown as positive (+) or negative (-) as Fig. 6.5 



Table 6.2 Coupling terms for SUMSBS classical applied fields 
with SPM pilot model 


Fields as Table 6.1. 
Model as Fig. 6.3. 
FORCE calculations 


Demand 
(model axes) 


F . 
y' 

M' 

N' 

Lift (F ,) 
z 

0.394 

0 

10“^ 

0 

Sideforce (F^,) 

0 

0.394 

0 

-0.01 

Pitch (M' ) 

-0.03 

o 

0.03 

O 

Yaw (N') 

0 

-0.08 

0 

0.03 

Primary terms (e.c 

j. lift force 

from lift 

demand) 

shown underlined 

Forces given as N 

per ampere 

(relevant E/Ms » in 

SUMSBS) 


Torques given as Nm per ampere (similarly) . 
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6.3 Suspension and calibration 


Stable suspension was achieved very easily, following triiraning 
of D.C. offsets within the controller (to provide standing currents 
to support the model's weight) and adjustments to various loop 
gains. This latter requirement was due to the large increase in 
the deadweight and moments of inertia of the model compared to the 
unwinged case. Suspension quality was good. 

Time only permitted the absolute minimum calibration to be 
attempted, that being measurements of roll torque versus E/M 
currents for a model at one position and attitude (the usual datum) . 
For convenience, roll torque was applied via a lightweight torque 
rod attached to the extreme nose of the model (Fig. 6.6). This 
technique has not been widely used in the past and has certain dis- 
advantages, perhaps principally the likelihood of applying axial 
forces to the model as well as roll torques. 



Only one E/M current could be accurately monitored with available 
equipment, but it is assumed that the current increment in all E/Ms 
due to the applied roll torque will be nearly equal. The restriction 
of the model to one fixed position and attitude (position/orientation 
error integrators operating in all control loops) further validates 
this approach. The line conductor approximation to SUMSBS (Fig. 6. 5) 
was again used to generate an estimate of the roll torque per ampere 
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(8 E/Ms) with the specified model. 

Results are presented in Pig. 6.7 with Figs. 6.8 - 6.10 
illustrating certain aspects of performance. 

6.4 Discussion 

The torque capability demonstrated is at least ten times higher 
than has been possible previously, under comparable conditions, 
with other SUMSBS roll control systems (11, 23). 

The calibration curve of Fig. 6.7 exhibits slight nonlinearities 
and a slope approximately 50% less than the computed estimate shown. 
However, the shortcomings in the experimental technique, the gross 
simplification of E/M (hence field) geometry in the computed E/M 
configuration and the partly unknown model magnetization distri- 
butions are thought to be sufficient to account for these discre- 
pancies and neither point is regarded as a matter for any concern 
at this stage. 

Scaling the measured calibration curve, using Eqn. 4.4, to a 
realistic aircraft geometry, taken to be the Cessna Citation I (a 
convenient choice with zero sweepback) indicates a torque capability 
of the order of : 

8.2 gm. cm. /Ampere ; perhaps + 80 gm.cms. normally usable: 

- the model span being taken as 4*j inches with representative core 
volumes. This is lower than the pilot SPM model due to the lower 
thickness and span, also the taper of the Citation wing. 

SPM roll control is concluded to be entirely practical for use 
with SUMSBS. 
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Fig,, 6.7 Experimental measurements of roll torque versus E/M current 
increment for pilot SPM model. 
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7. EXTREME ATTITUDE TESTING WITH MAGNETIC SUSPENSION AND 
BALANCE SYSTEMS 

7.1 Introduction 

Contemporary small MSBSs have not previously been capable of 
suspending models over a particularly wide range of attitudes, the 
axis of the model's fuselage typically being limited to approxi- 
mately 30 degrees excursion from the major axis (wind axis) of the 
MSBS (MIT(B)) or less. Whilst restricted ranges of sensitivity 
of position sensing systems represents the primary limit to such 
excursions in most cases, it is undoubtedly true that none of the 
existing MSBSs (apart from the modified SUMSBS) have been configured 
with the intention of suspending models significantly outside this 
range, indeed in the majority of cases the MSBSs would be magnetically 
incapable of realistic operation over a substantially wider range. 

It would be highly desirable for a LMSBS to be capable of 
supporting and restraining typical test models over a wide range of 
test attitudes under representative test conditions. This would 
enable rapid exploration of the aerodynamic characteristics of aircraft 
or missile models over their complete attitude envelope, including of 
course the regimes where the model may be thought of as being wholly 
or partially stalled. Such regimes are currently of great interest 
in military work and tend to be rather intractable analytically. 

Current test methods in this class generally require a range of model 
supports, each functional over a restricted range of attitudes. 

Support interference corrections can be virtually impossible since 
the supports often lie in regions of grossly turbulent wake. 

In the context of this report "extreme" attitudes are classified 
as any attitude outside the model axis excursion limit (30 deg.) 
mentioned above. It is thought that extreme attitude capability for 
a MSBS should include model axis angular excursions up to and 
including 90 degrees in at least one plane, permitting, for instance, 
model testing over the range of angle of attack of -90 to +90 degrees. 
Yaw and roll excursions need not be over such wide ranges but ideally 
would be so. Since the model's magnetic characteristics are com- 
pletely insensitive to the wind direction, combined pitch, yaw and 
roll excursions over a -90 to +90 degree range would encompass, with 
"reversal" of the model in the MSBS axes, the complete range of 
possible attitudes. 
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Several fundamental difficulties arise when considering model 
suspension at extreme attitudes. These include : 

1) Identification of E/M array geometries and configura- 

tions capable of generating, via field and field 
gradient components, forces and torques on the 
model in the required senses and magnitudes over the 
full range of model attitudes. 

2) Synthesis of control algorithms capable of accommodating 

large changes in model aerodynamic characteristics and 
magnetic couplings to the E/Ms. 

3) Design of position, attitude and other sensors to 

monitor wide ranges of model motion. 

Part of 1) above is addressed here, that is, the inclusion of 
adequate versatility into the E/M array configuration. Sizing the 
E/Ms thus specified to satisfy particular absolute force and torque 
requirements must be performed separately. 

7.2 Theoretical background 

7.2.1 Required field and field gradient components 


Forces and torques with conventional slender, axially magnetized 
models are assumed to be predominantly created by the following 
field and field gradient components (Appendix 5) : 


Table 7.1 Required field and field gradient components 
Force/torque Component 


^X- 

(axial force) 

H , 

o 


(sideforce) 

H , 

xy ' 
o 

F , 

(normal force) 

H , 

xz ' 
o 

L' 

(rolling torque) 

H , 

y^o 

M' 

(pitching torque) 

»Z- 

o 

N • 

(yawing torque) 

^o 

- 

(magnetizing field) 

»X- 

o 

- 

(magnetizing field) 

H 1 

yy^ 


Axes as Fig 7.1 


(spanwise magnets) 


(soft iron fuselage 
core) 

(soft iron spanwise 
magnets ) 
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In "normal" suspension model and balance axes coincide and 
these components correspond, neglecting the cases of roll and 
magnetizing fields, to: 

H H H H H 

XX , xy , xz , z , y 
o o o o o 


Fig. 7.1 Model and tunnel axis 
systems 


y' 



Tunnel axes — 
x,y,z 

Model axes — 

x',y',z' 


Sequence of rotations 
Yaw, pitch, roll 


o subscripts indicate 
the origin of axes. 


Pitching or yawing the model through 90 degrees translates 
these components into: 


H - H 

o o 


- H 


xz , 
o 


H 


H 


and 


H - H 




H - H , 

z^, x^ respectively. 


It may be immediately noticed that all nine primary field 
components : 

H 


,H ,H ,H ,H 

X y z XX xy XZ yy yz 


zz 
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are required independently at the origin for the full range of 
model attitudes to be usable 

The generalized couplings between field and field gradient 
components in tunnel and model axes may be calculated as follows: 

'Jhe transformation matrices for Yaw ('}'), Pitch ( 0 ) and Roll 


(^>) , in 

tunne 1 

axes are 

as follows: 





Costfj 

-Sintjj 

0 


CosG 

0 

Sin 0 


1 

0 

0 

Siniji 

Co sip 

0 


0 

1 

0 


0 

Cos(l) 

“Sin<j> 

0 

0 

1 


-Sin 0 

0 

CosG 


0 

Sincj) 

Cos<{) 


The couplings for a vector in model axes are given by the product 
of the reverse coupling matrices : 


V , 

X ' 


bosGCosi|j 

CosGSiniJ; 

-Sin 0 


:= 

Sin 0 Sin<{)Cosiji 

Cos(j)CosiJi 

Cos 0 Sin(|) 

y 


-Cos(J)SiniJj 

tSincfiSinGSiniJt 




SinctiSinij) 

Cos(l)Sin 0 Sin 4 ) 

Cos(j)Cos 0 



+Cos(t)Sin 0 Cosiji 

-Sin(t)Cosij; 



Which may be written as: 



The couplings for field components are given straightforwardly by: 


H’= AH - ( 7 . 2 ) 

Whereas for field gradient components : 

V = aV and using ( 7 . 2 ) 

V'.H' = aV.AH = bVh - ( 7 . 3 ) 

where the reduced form of B is given by : 


Ih aj^ 2aia2 2aias aa^ 

AX 

\y* aib2+ aib3+ a2b2 

aabi aabi 

\z' ^1*^1 ^1^2+ aiC3+ a2C2 

~ 32C1 asci 

2bjb2 2bib3 b2^ 

^’vz' biC2+ biC3+ b2C2 

b2Ci bsci 

2ciC2 2 ciC3 C2^ 


2 a 2 a 3 

32 b 3 + 
S 3 b 2 

32^3 + 

asC 2 

2 bab 3 

b 2 C 3 + 
t> 3 C 2 

2 c 2 C 3 



“( 7 . 3 a) 
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The matrix A being written as ; 



ai 

32 

33 


A 

= 

bi 

b2 

b3 




Cl 

C2 

C3 


It is seen that field components behave as vectors during axis 
rotations whereas field gradient components do not and that the 
couplings from applied field gradients to field gradients in model 
axes are extremely complex. Certain simplifications are possible 
where rotation is restricted to one plane, say pitch, where ; 


if) , ^ = 0 ; 


and; 


1 H , 

XX ' 


xy ' 


H 1 

xz' 


H 1 

yy' 


yz* 


^ 1 

zz* 



Cos^G 

0 

CosGS: 

0 

0 

Sin^G 


F = dV 


|h .1 


CosG 0 

-SinG 










— 

0 1 

0 


1 

H 

y' 





y 



SinG 0 

CosG j 


H 

z'l 



1 


i ^ 


0 

-2CosGSinG 

0 

0 


CosG 

0 

0 -SinG 


nG 

0 

Cos^e-Sin^e 

0 

0 

-1 


0 

0 

1 

0 


SinG 

0 

0 CosG 



0 

2Sin6CosG 

0 

0 

( 




- 

(7. 

3b) 

d torques 

can be given 

by any 

of 

the 

and 


T = /M X H + 
— V — — 

■ £ x(M.Vh) 

dv 


- (7.2a) 


0 
0 

Cos^G 


H 

XX 


xy 


xz 


yy 


yz 


zz 


- following Eqns. 3.1 


The force and torque couplings will differ from 7.2, 7.3 etc., due 
to the effects of the spatial variations of the applied fields over 
the volume of the model's core. Where the applied fields are 
relatively uniform over the core these effects will be of second 
order and the force and torque couplings will be approximated by 
7.2, 7.3 etc. 

It is considerably more convenient in many cases to represent the 
couplings by transforming model magnetizations into tunnel axes 
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rather than field and field gradient components into model axes : 


F 


H 

H 

H 


M 

X 


XX 

xy 

XZ 


X 

F 

=: V 

H 

H 

H 


M 

y 

xy 

yy 

yz 


y 

F 


H 

H 

H 


M 

z 


xz 

yz 

ZZ 


z 

L 


0 

H 

-H 


M 




z 



X 

M 

~ V 

-H^ 

0 

H 


M 



X 


y 

N 


H 

-H 

0 J 


M 



y 

X 

' 


^1 

magnetizations 

being 

given by : 


M 


$ 0 

4' 

M' 



- (7.4) 


+ / r x6f term 

V — — 


7.2.2 E/M configurations for multiple independent field and 
field gradient component generation . 

The requirement to generate 9 field or field gradient components 
independently necessitates at least 9 independent E/Ms. The desire 
for symmetry in the E/M array (Section 1.2.1) acts to increase this 
figure. 

A symmetric quadruplet of E/Ms as shown in Fig. 7.2 can generate 

four field gradient components at the origin;H (Ii , I 4 = -I 2 / Is)* 

^o 

di, I2* 13= I4)* H di, I2, I3 = I4), and 
o ^o 

H (Ii * I? I Is= l4)»but it is immediately seen that H , H 
zz XX yy 

o 0 0 

and H are not independent. Two field components at the origin; 
"o 


H d;i, I 2 = “l 3 f I 4 ) and H (Ii, I3 = -l 2 * l 4 ),may also be 
o o 

generated. Modified geometries (Fig 7.3), in fact corresponding 

to the "vertical" or "lateral" E/Ms in SUMSBS, are found to be 

relatively weak in H as shown. If H (or H ) were regarded 

XX zz yy 

o 0 0 

as a prime component of field for this sub-configuration and the 

' stray ' components H and H (H ) were countered by some other 
XX yy zz 

o o o 

means, the modified quadruplet would be a useful generator of four 

independent field or field gradient components jH , H , H , 

00 o 

H , (H ) . 

zz,- yy^ 

Ll O 

A symmetric pair of E/Ms can generate one field and one field 

gradient component independently at the origin H , and H as 

o o 

shown in Fig. 7.4, with stray components H and H 

yy zz 

o o 
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Conventional E/M configurations can be considered as an assembly 
of quadruplets and pairs as defined above and it is straightforward 
to arrange for all required field and field gradient components to 


be somehow provided. 

In fact, current MSBS configurations (eg. MIT(B), SUMSBS) are 

frequently capable of generating 8 field or field gradient components 

(H , H I, H , H , H , H , H , H ) , effectively though not 
X y z XX xy xz yy zz 
o o o o o o •' o o 

completely independently and typically lack capability only in H 

^^o 

It may thus be expected that conventional configurations are 
capable of sensible operation over a wider range of attitudes than 
has been exploited previously. 


7.3 Identification of maximum force and torque capabilities of 
a MSBS 

Any extreme attitude capable MSBS will incorporate a number of 
independent E/Ms, perhaps between 10 and 20. Each E/M tends to 
create all nine independent field and field gradient components at 
general points. Thus, with a suspended model at a general attitude 
and position, each E/M creates force and torque components along or 
about each axis (x, y, z) . It is clear that there will seldom exist 
a unique solution for E/M currents with any particular force and 
torque requirement. The absolute maximum force and torque capa- 
bility of a particular system may not, therefore, be directly 
analytic. 

With 6 independent E/Ms however, there is no serious difficulty 
since, neglecting demagnetization or saturation effects (see Section 
7.7), the required senses for force and torque specify variables 
kj-ks such that : 

kiF = = kjL = k^M = kjN = . (7.5J 


F 

X 


F 

y 


F 

z 


L 


M 


N / 



= 

kgF 

z 

II 

L = 

k4M 

of 

the 

6 E/Ms may 1 

^1 

^2 

^3 

^4 

as 

hi 

b2 

h >3 

^4 

hs 

Cl 

C2 

C 3 

C4 

C5 

di 

d2 

dj 

d 4 

ds 

Cl 

C2 

®3 

64 

65 

fl 

f2 

f3 

£4 

fs 


where ai-feare dependent on the systems geometry and magnetic 
characteristics and are functions of model position and attitude; 

Il“ le representing the current levels in each E/M. It is immediately 
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seen that where aj-feare mostly non-zero, Iz-Ie (say) are dependent 
variables with only Ii unspecified. Since all E/M currents are now 
fixed relative to each other, the absolute maximum capability is 
found by increasing all currents until one reaches its limiting 
value, that value being determined by the E/M design. 

Useful results may be obtained for the case of 7 E/Ms, starting 

with the equations for a 6 E/M system operating with one current 

P 

(say Ii) at its limiting value (Ii ), as above : 


1 F 1 


lai 

32 

as 

ai. 

as 

ae \ 



F 

y 


1 bi 

b2 

b3 

b4 

bs 

be 


la 


F 


Cl 

C2 

C3 

C4 

C5 

C6 


1 

I 3 


Z 

L 


di 

d2 

d3 

d4 

ds 

de 


I 4 

- (7.7) 

M 1 


ei 

S2 

63 

©4 

65 

ee , 




r / 



f2 

fa 

f4 

fs 



(i. 1 


1 6 being determined 

using Eqns. 

7.7 

and 

7.5. 

Introducing 


another E/M (I 7 ), the equations may be reformed : 


/f ' 

X 




la 2 

as 

34 

as 

ae 

a? 1 

1 I'M 

F 

y 


(bl 



bs 

b4 

bs 

be 

by 


1 I 3 

F 

z 

-I^p 

Cl 

=r 

C2 

C 3 

C4 

cs 

ce 

cy 


I 4 

L 


dl 


dz 

d3 

d4 

ds 

de 

dy 


Is 

M 1 


eij 


62 

63 

64 

65 

6e 

6,1 



n/ 


I'V i 


f3 

£4 

fs 

fe 

^’1 


U’l 


Since the coefficients aj-fi and are constants for any 
particular model attitude and position, resolution of 7.8 follows 
a similar procedure to 7.6. This leads to a second current 

p 

reaching its limiting value, say I 2 • 

If the number of E/Ms is represented by n, then the cases above 
yield the criteria that where n =6 or n=7, n-5 E/Ms should be 
operated at their peak design currents to achieve a particular 
maximum force and torque capability. It is now argued that this 
criteria holds for all n^ 6 . It is not, however, obvious which 
n-5 E/Ms should be chosen to be at their respective peak currents, 
but this can be resolved by computing all permutations. 

Cases where the matrix in 7.6 etc. is sparse may be dealt with 
by reduced forms of the analysis above. 
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7.4 


+" E/M configuration study 


7.4.1 Geometry 

The revised configuration for SUMSBS can fall into this category 
(see Sections 1.2.1 and 2.2.1). The performance of this system 
cannot be explicitly computed by FORCE (Appendix 2) since the main 
E/Ms are iron cored, however a system having similar proportions 
but arbitrary exact dimensions has been computed to yield an 
indication of the usable attitude range of SUMSBS (Fig. 7.5). 

This representative system is chosen to be at ten times the scale 
of SUMSBS for convenience. Single loop E/Ms are used to simplify 
calculations but are expected to provide adequate approximations of 
the performance of SUMSBS' E/Ms, the current levels being adjusted 
where appropriate to approximate the relative field capabilities of 
the E/Ms in SUMSBS. 

7.4.2 Pe r f ormance 

As might be expected from Eqns. 7.3 etc., the generalized couplings 
between E/Ms and model core are likely to be extremely complex. 

Since the inherent limitations of the optical position sensing system 
used in SUMSBS would, at least initially, restrict large angular 
displacements to one plane, specifically the xz (pitch) plane, all 
computations presented here are made with the model restricted to 
this plane, with its centroid fixed at the origin of balance axes. 

Thus x=y=z=ij;=(t)=0. 

With or without the above restriction, the magnetic couplings 
between representative E/Ms (Fig. 7.5) and a typical model (Fig. 7.6) 
are easily found for given model attitudes and positions using FORCE. 

It is not, however, clear in what form the couplings should be 
presented; as couplings from individual E/Ms (Fig. 7.7a-d) or as 
couplings from conventional groupings of E/Ms (Fig. 7.8a-j). In 

this case, the latter approach fully separates the creation of 
magnetic forces and torques but is less general, partly because the 
groups of E/Ms are not mutually exclusive. 

The effects of the spatial variations of the applied fields around 
the origin of axes may be found by recomputing Figs. 7.8 with model 
polarization concentrated into a single dipole at the origin (Figs. 
7.9a-j). It is found that these results are identical to direct 
computation from Eqns. 7.3 etc. with the field capabilities of the 
E/M groupings (calculated) as: 
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All E/Ms are line conductors 
All currents are positive as shown 


All E/M curregt 
limits are 10 amps. 




< 


E/Ms 

approximatq^^^ 
as 32 straight- 
line lementls 



<r. (2) 



^ Current limits of E/Ms 9,10 

are reduced to 598,400 amps 
for calculation of 7.13 etc 

Fig. 7.5 Geometry of + configuration 

Dimensions in metres. This config- 
uration represents an enlarged 
representation of SUMSBS 



Angle\ 


of (0) 
attack / 



Fig. 7.6 Geometry of model 

for + and X studies 

Dimensions in metres 
Polarization = 1 Tesla 
20 model elements 



Fig.. 7.7 Couplings from + configuration individual E/Ms 


E/Ms as Fig. 7.5 
Model as Fig. 7.6 

Data is shown specifically for only 4 E/Ms. Remaining E/Ms may be 
deduced from this data by symmetry, resulting in the sign table below: 

Signs of force /torque components in: 

Fig. 7. 7a Fig. 7.7b 



Fig. 7. 7c Fig.7.7d 

E/M F . F„ L M N E/M F F M 

' X Y Z X z 

2 (shown) - - + - - - 9 (shown) + + 

4 “ + + + — + 10 + + + 

6 — .j. « 4. — 

8 -- + + + + 



Fig. 7.7b 

+ Configuration 
E/M 3 (5) 


Fig. 7.7c(i) 

+ configuration 
E/M 2 (4,6,8) 

Forces 

FORCE 

(N) 


a 





Fig. 7.8 Couplings from + configuration classical applied fields 


E/Ms as Fig. 7.5 
Model as Fig. 7.6 


Fields are defined by the relationship of E/M currents as shown below: 


Field 

II 

l2 

I3 

l4 

l5 

l6 

I7 

1 8 

l9 

IlO 

Lift 

+ 

N/A 

- 

N/A 

- 

N/A 

+ 

N/A 

N/A 

N/A 

Drag 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 


•f 

Auxiliary 
drag No.l 

- 

N/A 

- 

N/A 

— 

N/A 

— 

N/A 

N/A 

N/A 

Auxiliary 
drag No. 2 

N/A 

— 

N/A 

- 

N/A 

— 

N/A 

— 

N/A 

N/A 

Sideforce 

N/A 

- 

N/A 

+ 

N/A 

+ 

N/A 

- 

N/A 

N/A 

Pitch 

- 

N/A 

+ 

N/A 

- 

N/A 

+ 

N/A 

N/A 

N/A 

Auxiliary 
pitch No.l 

- 

N/A 

— 

N/A 

+ 

N/A 

+ 

N/A 

N/A 

N/A 

Auxiliary 
pitch No. 2 

N/A 

— 

N/A 

— 

N/A 

+ 

N/A 

+ 

N/A 

N/A 

Axial 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

- 

+ 

Yaw 

N/A 

- 

N/A 

+ 

N/A 

- 

N/A 

+ 

N/A 

N/A 


- following Fig. 7.5 


This representation of applied fields is derived from the fields 
classically applied with the model at its datum position (lift, pitch 
etc.), augmented so as to encompass all possible combinations of 
currents. 

The separation of Auxiliary Drag into two constituents (Figs. 7. 8c 
s d) is done since simultaneous application of these two fields with 
one in the reverse sense produces the through-wing field used in SUMSBS 
for SPM roll control. 
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Fig. 7.8g 


3000 


Auxiliary pitch 
No.l field 



Fig. 7.8h 

Auxiliary pitch 
No. 2 field 
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standard model, and are 
identical to the corre- 
sponding Figs. 7.8 


Fig. 7.9b 

+ configuration 
Drag field 
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Fig. 7.9c 


2100 


+ configuration 
Auxiliary drag 
No.l field 



Fig. 7.9 d 

+ configuration 
Auxiliary drag 
No, 2 field 
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L68 



Fig. 7.9g 

+ configuration 
Auxiliary pitch 
No.l field 


Fig. 7.9h 

+ configuration 
Auxiliary pitch 
No,. 2 field 




ANSLE OF ATTACK CDESS 




0 


Table 7.2 Field capabilities for + configuration 


Field components (T) / lO^AT /E/M 


E/M grouping 

B 

XX 

o 

B 

^^o 

B 

XZo 

B 

yvo 

B 

yzo 

B 

ZZo 

B 

x,y,zo 

Lift 

0 

0 

0,1969 

o 

0 

O 

0 

Drag 

0.0676 

0 

0 

-0.0338 

0 

-0.0338 

0 

Auxiliary 

0.002 

0 

0 

0.1252 

0 

-0.1272 

0 

drag 1 
Auxiliary 

0.002 

0 

0 

-0.1272 

0 

0.1252 

0 

drag 2 
Sideforce 

0 

0,1969 

0 

0 

0 

0 

0 


B 

^o 

B 

^o 

B 

^O 

All gradients 



Pitch 

0 

0 

-0.1267 

0 




Auxiliary 

-0.1147 

0 

O 

0 




pitch 1 
Auxiliary 

-0.1147 

0 

0 

0 




pitch 2 
Axial 

-0.0808 

0 

0 

0 




Yaw 

0 

0.1267 

0 

0 





The correlation between the couplings for the representative model 
and the zero volume model, the latter being essentially theoretical 
results, is generally fair, the cases where the E/Ms are relatively 
remote from the model (Figs. 7.9b, 7.9h) exhibiting the best agreement. 
The "theoretical" curves may thus be adequate for synthesis of experi- 
mental control algorithms, provided generous stability margins are 
available and such algorithms would not then require updating 
following minor model changes. Peak performance of the system is, 
however, unlikely to be approached unless the computed couplings of 
Fig. 7.8 or similar are used directly. 

The variations of the maximum force and torque capability can be 
found from the data in Fig 7.7 using the methods of Section 7.3. 

Again there is difficulty due to the range of possible combinations 
of model forces and torques. However magnetic torque capability, 
apart from roll torque, is seldom a critical factor in MSBSs and 
with model excursions restricted to the pitch plane, sideforce 
capability may be largely neglected. Figs. 7.l0(a-b) thus illustrate 
the mciximum force capability for the representative + system with the 
force vector constrained to lie in the plane of pitching and with all 
torques zero. The sense of the force vector is clarified in Fig. 7. 11. 
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FORCE CN> FORCE <N> 


Fig. 7.10 Maximum performance of the + configuration 
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The direction of the magnetic force vector required to oppose the 
net aerodynamic force would normally be expected to lie in the quadrant 
defined by 0$3<90°. Transient forces may be outside this range, as 
may be the model ' s weight vector . However, for the purpose of 

identifying the absolute maximum force capabilities of a MSBS, and 
since the curves of Fig. 7.10 are repeated for 3*‘3±180° the 0^3^90° 
quadrant is seen to be of principal interest. Some of the content 
of Fig. 7.10 may now be simplified by choosing three components of 
aerodynamic significance, namely pure lift force, pure drag force and 
normal force, all defined in Fig. 7.12, and adding the envelopes 
within which all curves for 0^3^90° must lie (Fig. 7.13). The axial 
E/M current limits (E/Ms 9 and 10) have been factored by 0,5984 to 
approximate the measured relative performance of the main and axial 
E/Ms in SUMSBS. Computation for the zero volume model is illustrated 
in Fig. 7.14 and exhibits significant differences in certain areas, 
such as drag forces at low angles of attack. Effects of variations 
in the E/M array geometry are well illustrated with this form of 
presentation. For instance, doubling the current limit applied to 
the axial E/Ms (E/Ms 9 and 10) increases, as would be expected, the 
drag force capability at small angles of attack but also attenuates 
the otherwise rapid fall off in normal force capability with angle of 
attack increasing from small angles (Fig. 7.15). 

The discontinuities in slope of the curves in Figs. 7.13 - 7.15 
are caused by changes in the distribution of E/M current limiting, for 
example at some particular angle of attack E/Ms 1 and 7 (say) may be 
operating below their current limits, with all other E/Ms at their 
respective limits (see Section 7.3). Small increases in angle of 
attack may result in the currents in E/Ms 1 and 7 increasing, other 
currents remaining fixed, and at some angle all E/Ms will be at their 
limiting currents. Further increases in angle of attack will 
generally require E/Ms other than E/Ms 1 and 7 to be reduced in 
current, overall performance then being delineated effectively by 
a different performance curve. This may be illustrated for one case, 
chosen to be the normal force curve from Fig. 7.13, Fig. 7.17 clearly 
showing tne changing E/M current distributions. 

7.4.3 Discussion 

Reserving treatment of the significance of the absolute magnitudes 
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Fig. 7.13 Maximum performance of the + configuration 

E/Ms as Fig. 7.5 Model as Fig. 7.6 Force vectors as Fig. 7.12 



'4 





Pig. 7.14 Maximum performance of the + configuration with zero 



Pig. 7.15 Maximum performance of the + configuration with E/Ms 


9,10 doubled in strength. Otherwise as Fig. 7. 13 
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This represents one possible method of alleviating the poor 
performance of the + configuration at 90 angle of attack, but 
has not been fully investigated. ^ 

Skewed lateral E/Ms produce sideforce at 90° 









of forces and moments in Figs, 7.7 - 7.10 and 7.13 - 7.15 for 
Section 7.7, the computed couplings do nevertheless indicate the 
practicality of the particular + configuration geometry studied 
for high angle of attack testing. , 

Figs. 7.10 and 7.13 clearly illustrate that the principal aero- 
dynamic forces (appearing in the plane of pitching) can be 
countered relatively effectively at most angles of attack with 
only the normal and drag force minima around 30° angle of attack 
presenting any serious difficultyi Augmentation of the main 
drag field partially alleviates this problem, the severity of which 
is heavily dependent on the precise geometry (see (49)). Without 
indulging in further analysis it is concluded from Figs. 7.8f-j 
that sufficient capability for the purposes of generating pitch and 
yaw torques (roll not being treated here) exists at all angles of 
attack, leaving only one fundamental flaw in the E/M configuration, 
that is the inability to generate sideforce at 90° angle of attack 
(Fig. 7.8e). This is due to the lack of capability in 
was shown in Section 7.2.1 to be required for this purpose. 

The particular geometry studied appears to be practical for 
operation up to at least 60° angle of attack, then being limited 
principally by falling sideforce capability. It is thought that 
this view remains broadly valid for all + configurations of con- 
ventional proportions (49) , though detail optimisation of perfor- 
mance is clearly possible. One example of this might be the in- 
corporation of "skew" into the lateral E/Ms (Fig. 7.16), although 
this technique does partially destroy the symmetry of the 
configuration. 

It is important to recall at this stage that the field component 
required for generation of roll torque using the SPM or SIM schemes 
is array were augmented in such a way so as to 

include capability for this component, then operation to 90° angle 
of attack and, in fact, through 360° pitch attitude range would 
appear possible. 

The somewhat bizarre fluctuations in E/M currents apparently 
required to achieve maximum force capability with rising angles of 
attack (Fig. 7.17) are almost certainly impractical due to the low 
E/M current slewing rates available in realistic MSBSs. However, 
the complete reversals of current occurring with small increases in 
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angle of attack at certain points in Fig. 7.17 generally indicate 
that the E/M in question is relatively ineffective in generating 
force in the required sense at that particular attitude. Approxi- 
mation to the absolute maximvim performance with appropriate 
restrictions on current fluctuations are thus expected to be possible 
with relatively small performance degradation, though no analysis 
has yet been attempted. 

7.5 "X" E/M configuration study 

7.5.1 Geometry 

The lack of capability in of the + configuration can be 
rectified by simply rotating the model's axis system (in the roll 
sense) so as to convert the + configuration to an X configuration. 

A representation of this version of SUMSBS has been computed, using 
identical E/M and model geometry and characteristics as the represen- 
tative + system of Section 7.4 (Fig. 7.18). 

7.5.2 Performance 

Following identical analytic procedures to those employed for the 
+ system. Figs. 7.19a-b show the main E/M couplings, those from the 
axial E/Ms (E/Ms 9 and 10) remaining as shown in Fig. 7.7d. Figs. 
7.20a-g show the couplings from classical E/M groupings, with 
couplings to the zero volume model presented in Figs. 7.21a-g. 
Immediately it can be seen that sideforce capability exists at 90° 
angle of attack (Figs. 7.20/7.21d). 

Computing the maximum force capabilities in the plane of pitching 
yields a disappointing result (Fig. 7.22). Although the lift force 
capability at zero angle of attack is augmented, compared to the + 
configuration, as would be expected since all 8 main E/Ms now contri- 
bute to lift force, the performance around 45° angle of attack is 
disastrous. The zero volume model exhibits similar results 
(Fig. 7„23), whilst doubling the axial E/M strength effects a 
significant but probably inadequate improvement (Fig. 7.24). 

Eqns. 7.3 etc. may be used to illustrate the reason for the 
poor performance around 45°, using the (calculated) field 
capabilities of the E/M groupings : 
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+ geometry as shown 
in Fig. 7.5 


E/Ms (3) & 7 



& 6 


X geometry 



/Ms (3) & 7 


E/Ms (1) 
& 5 


E/Ms (2) & 6 


z 

Fig. 7.18 Derivation of X configuration geometry 

Both geometries shown schematically and in Section through the 

yz plane 
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Fig. 7.19 Couplings from X configuration individual E/Ms 


E/Ms as Fig. 7.18 
Model as Fig. 7.6 

Data is shown specifically for only 2 E/Ms. E/Ms 9 and 10 remain 
as shown in Fig. 7.7d. Remaining E/Ms may be deduced from the data 
shown by synunetry, resulting in the sign table below; 


Signs of force/torque components in Fig. 7.19a 


E/M 

F 

F 

F 

L 

M 

N 

X 

y 

z 




1 

- 

+ 

+ 

+ 

- 

4 - 

2 

- 

- 

4 

- 

- 

- 

7 

- 

+ 

4 

- 

+ 

- 

8 

- 

- 

4 

+ 

+ 

4- 

Signs 

of force/torque 

components in 

Pig. 

7.19b 

E/M 

1 F 

F 

F 

L 

M 

N 

1 X 

y 

z 




3 

+ 

- 

- 

_ 

- 

- 

4 

+ 

+ 

- 

4- 

- 

+ 

5 

+ 

- 

- 

+ 

+ 

4- 

6 

+ 

+ 

- 

- 

4 - 

- 


( shown ) 


(shown) 
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Fig. 7.19a (i) 

X configuration 
E/M 1 (2,7,8) 
Forces 


Fig. 7.19a(ii) 

X configuration 
E/M 1 (2,7,8) 
Torques 
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6 i0se90 4aEa«o7a6aoa 
ANGLE OF attack CDEG5 
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Fig. 7.19b(i) 

X configuration 
E/K 3 (4,S,6) 
Forces 


Fig. 7.19b(ii) 

X configuration 
E/M 3 (4,5,6) 
Torques 
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Fig. 7.20 Couplings from X configuration classical applied fields 


E/Ms as Fig. 7.18 
Model as Fig. 7.6 

Fields defined by relationships of currents as below: 



Ii 

l2 

I3 

l4 

Is 

Is I? 

Is 

I9 

IlO 

Lift 

+ 

+ 

- 

- 

- 

+ 

+ 

N/A 

N/A 

Auxiliary drag 

- 

- 

- 

- 

+ 

+ + 

+ 

N/A 

N/A 

Sideforce 

+ 

- 

- 

+ 

- 

+ + 

- 

N/A 

N/A 

Auxiliary sideforce 

+ 

- 

+ 

- 

+ 

+ 

- 

N/A 

N/A 

Pitch 

+ 

+ 

- 

- 

- 

+ 

+ 

N/A 

N/A 

Auxiliary pitch 

+ 

+ 

+ 

+ 

- 

- 

- 

N/A 

N/A 

Yaw 

+ 

- 

- 

+ 

+ 

- 

+ 

N/A 

N/A 




- 

following Figs. 7.5 

7.: 

18 


Axial fields (E/Ms ! 

9 & 10) are 

! as 

Figs. 

7.8 

b & i 





Fig. 7.20a 
Lift field 
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Fig. 7. 20b 
Auxiliary drag 


Fig. 7.20c 
Sideforce field 
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Fig. 7.20g 


Yaw field 
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Fig. 7.21 Couplings from X configuration classical applied fields 
- showing zer o volume model 

These Figures are directly comparable with I'igs. 7.20, which 
should be consulted for details of model, E/Ms and current relation- 
ships . 

SOLID lines represent the zero volume model. 

BROKEN lines show the standard model and are identical to the 
corresponding Figs. 7.20. 
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F ig. 7.21b 

X configuration 
Auxiliary drag 


Fig. 7.21c 

X configuration 
Sideforce field 
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Fig. 7.21£ 

X configuration 
Auxiliary pitch 


Fig, 7.21g 

X configuration 
Yaw field 
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Fig. 7 <22 Maximum performance of the X configuration E/Ms as Fig. 7. 18 



Fig. 7. 23 Maximum performance of X configuration with zero volume model 

Otherwise as Fig. 7.22 
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FORCE CN> 


Fig. 7.24 Maximum performance of the X configuration. Axial E/Ms 
doubled in strength. Otherwise as Fig. 7. 22 



“ 193 



Table 7.3 Field capabilities for X configuration 
Field components (T) / lO^AT / E/M 


E/M grouping 

B 

^^o 

B 

XV 

■^o 

B B 

xz yy 

o ■'■'o 

B 

yz 
■' o 

B 

^^o 

B 

x,y,z^ 

Lift 

o 

0 

0.2786 0 

0 

0 

0 

Drag 

0.0676 

0 

O -0.0338 

0 

-0.0338 

0 

Auxiliary drag 

0.0039 

0 

0 -0.0019 

0 

-0.0019 

0 

Sideforce 

0 

0.2786 

0 0 

0 

0 

0 

Auxiliary 

0 

0 

0 0 

-0.2525 

0 

0 

sideforce 








B 

^o 

o 

>1 

in 

B All gradients 

^o 



Pitch 

0 

0 

-0.1793 0 




Auxiliary 

-0.2296 

0 

O 0 




pitch 







Axial 

-0.0808 

o 

O O 




Yaw 

0 

0.1793 

O 0 





Eqn. 7.4 gives, for the x and z components of force in this case; 


F = V(HM + HM) 

X XX X xz z 

F = V(HM + HM) 

Z xz X zz z 

and since M =M Cos9 and M =-MSin0 we have, at 45° angle of attack; 
X z 


' s: 

k 

(H 

- H 

X 


XX 

XZ 

1 23 

k 

(H 

- H 

z 


XZ 

ZZ 


{ k = MV/ /2) 


- (7.9) 


These equations are heavily coupled, directly via the term 

and via the lack of independent capability in or shown in 

the table above. The requirement is generally to develop positive 

F^ and positive F^ simultaneously in some proportion and it is seen 

that since H contributes to F and F in opposite senses and that 
xz X z 

the (coupled) and components are relatively small in the 
chosen configuration the available magnitudes of forces will tend 
to be low. 

7.5.3. Discussion 

Referring to Tables 7.2 and 7.3 it is seen that where H , H 

XX yy 

and components only are developed at the origin, the sum of 
these components is always zero. This corresponds to Maxwell's 
equations for a conservative field with no free poles : 
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7.6 E/M configurations for LMSBSs 


The examples of + and X configurations studied herein cannot be 
assumed to be fully representative of the general classes of such 
configurations, since large variations in overall proportions and 


relative E/M capabilities are possible. However, the inability to 


generate 

configuration. 


hence sideforce at 90° angle of attack with the + 
and the difficulties around 45° angle of attack with 


the symmetric X configuration are thought to be fundamental . 

There is no objection in principle to combining the two configura- 
tions in an attempt to combine their advantages and negate their 
disadvantages, indeed this was envisaged (but never carried out) in 
Ref. 49. Numerous detail variations are possible, such as : 




and may itself be subject to considerable variation) 
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Configurations a) and b) above correspond closely to proposed 
SIM roll arrays (Section 5.8). Indeed, although roll control has 
generally not been treated in this Section it is the case that the 
X array as specified can effectively generate the required through 
wing fields (H ,) for all model angles of attack. Large roll dis- 

placements from datum would not be permissible with the X array, 
however, since H , cannot be generated at (1)=45° (Section 4,3.1). 

The hybridised arrays above (a-c) do not suffer this limitation. 

The employment of model core geometries other than the conventional 
slender, axially magnetized configuration is possible, particularly 
with permanent magnet core material, but would generally lead to other 
complications, including increases in the complexity and/or severity 
of cross couplings between model degrees of freedom. 


7. , 7 Discussion 

The absolute magnitudes of forces (or torques) shown in the various 
Figures in this chapter are of somewhat limited significance. They 
do not represent the maximum attainable forces for the chosen con- 
figuration, merely the forces attainable within the constraints of 
given E/M current limits and model characteristics. The absolute 
maximum sittainable forces with a permanent magnet core are principally 
set by the magnetic behaviour of the core (demagnetization and usable 
polarization) but these limits could not be identified at the time of 
writing due to the lack of data concerning high coercivity permanent 
magnet materials (see Section 4.3). The corresponding limits for a 
soft iron core will be set by the behaviour of the induced magneti- 
zations as the core material enters its saturation region, presenting 
a challenging computational problem (see Section 5) . 

The studies made herein are principally intended to indicate the 
general practicality of the + and X configurations, partly by 
identification of any magnetic shortcomings in these configurations, 
and the magnitude of the control difficulties caused by the changing 
couplings between model and E/Ms. Further investigation of the latter 
area appears appropriate but cannot be undertaken here. 

In the design of an extreme attitude capable MSBS it is probably 
more appropriate to consider the E/M array simply as a generator of 
field and field gradient components in the test section, rather than 
an assembly of "lift", "drag", etc., E/Ms, as has been common 
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previously. Usable model attitudes are unlimited if all components 
can be generated effectively and independently. It would be to 
advantage if the applied fields were relatively pure and uniform 
in the region of the model in order to minimise cross coupling 
effects and this requirement may perhaps be met by the use of 
Helmholtz pairs and corresponding optimised quadruplet geometries. 

7.8 Experimental demonstration of suspension with SUMSBS at high 
angles of attack 

7.8.1 Introduction 

Section 7.4 concluded that SUMSBS should be usable, in its + 
configuration, to around 60° angle of attack with conventional models. 
Preliminary attempts at suspending models up to this angle of attack 
have been made and are reported briefly here. 

7.8.2 Special features of SUMSBS for high angle of attack operation 

The position sensing system, fully described in Appendix 1, 
represents only an expedient and very short term solution to the 
difficult problem of accurately monitoring model position and 
attitude over wide ranges of attitude. The system devised achieves 
high angle of attack capability only by mechanical rotation and 
realignment of the four main sensors, axial position being monitored 
by a single fixed sensor (Appendix 1, Fig. 7.26). 

Without modification, SUMSBS was capable of suspending conventional 
axisymmetric models to around 30° angle of attack. 

It was decided to attempt initial high angle of attack suspension 
at around 50°, requiring suitable relocation of position sensors of 
course, and decoupling of the now grossly coupled and quite inoperable 
controller. Decoupling proceeded using the computed couplings of 
Section 7.4, and was performed in model axes, since the position 
sensing system operates essentially in this axis system. The 
precise methods used are described below, but it should be carefully 
noted that they represent a considerable simplification from the 
ideal, made entirely for convenience, it being believed that 
successful stable suspension would be sufficient to verify the 
fundamental theory of this Section, even if achieved on a less than 
wholly representative basis. 
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Fig. 7„8j shows that yaw torque {model axes of course) is only 
weakly affected by pitch rotation, with Fig. 7.8e furthur indicating 
that sideforce is only affected by falling gain, easily restored by 
appropriate increase in controller loop gain. Consulting Fig. 7.7a 
it is seen that if E/Ms 1 and 7 are operated in conjunction, but with 
opposite signs of currents, then only pitch torque is generated and 
relatively effectively. The full pitch torque capability previously 
enjoyed is restored by supplementing the effect of E'/Ms 1 and 7 
by the use of E/Ms 2,4,6 and 8 via the Auxiliary Pitch No. 2 field, as 
Fig. 7.8h. 

The residual difficulty lies in the generation of forces in the 
required senses and magnitudes in the xz plane. Now Fig. 7.7d 
indicates that the use of E/Ms 9 and 10 (the axial E/Ms) in conjunc- 
tion produces a net force approximately perpendicular to the x' axis 
of the model, with the model at 50° angle of attack (Fig. 7.27). 

Forces approximately parallel to the model's x' axis are conveniently 
provided by E/Ms 3 and 5, with these E/Ms operated in conjunction, as 
Fig. 7.7b and Fig. 7.27. It was found that the currents required in 
E/Ms 9 and 10 required to oppose the relevant component of the model's 
weight proved to be very close to the available current limits 
(20 A). A standing current was thus demanded in E/Ms 2,4,6 and 8 
corresponding to the Auxiliary Drag No. 2 field (Fig.7.8d), effectively 
reducing the residual model weight required to be supported otherwise. 

Careful adjustment of controller loop gains was necessary before 
stable suspension could be achieved, but the "50°'' controller 
subsequently proved capable of suspension over the range of angle of 
attack of 40° to 60°. No calibration or detail examination or 
optimisation of performance was attempted but the Figs. 7.28a-g 
illustrate the full range of angle of attack explored. 

Suspension quality at the higher angles of attack was poor, perhaps 
not surprisingly in view of the highly simplified experimental 
approach, and considerable further study is thought appropriate, some 
aspects of which are discussed in Section 9.2.2. 

However Figs. 7.28 are thought to validate the fundamental and key 
conclusion of this Section, that is that the use of MSBSs for high 
angle of attack or extreme attitude wind tunnel testing is absolutely 
viable . 


199 - 




Fig. 7.27 Schematic diagram showing creation of forces in the xz 
plane at 50° angle of attack 























8 . 


RELIABILITY OF LARGE MAGNETIC SUSPENSION AND BALANCE SYSTEMS 


8.1 Introduction, the need for reliability 

MSBSs are inherently complex systems, incorporating a wide 
variety of sub-systems such as E/Ms, perhaps with associated 
cryogenic plant, power supplies, at least one set of model position 
sensors, analogue or digital closed-loop control systems and, of 
course, a wind tunnel and its ancillaries. Due largely to this 
complexity, all existing MSBSs are rather "unreliable" in the sense 
that control over model position is lost relatively frequently, 
perhaps after every few hours of operation. Model "flyaway" 
inevitably follows, since MSBSs of the types considered to date are 
open-loop unstable in at least one model degree of freedom (Earnshaw's 
Theorem) . Flyaway is seldom more than inconvenient with small 
systems since models tend to be small, relatively simple and 
inexpensive and operated with low tunnel dynamic pressures. 

In the case of a LMSBS, however, high dynamic pressures may be 
required to satisfy demands for Reynolds number and Mach number, 
and winged models are likely to be operated in high lift configura- 
tions at these dynamic pressures, leading to large net forces. 

Despite the fact that models will be relatively massive, since it is 
essential to occupy a high proportion of model volume with magnetic 
matericil, the typical maximum aerodynamic force will exceed the 
model's weight by a large factor (approximately 10 in Ref. 5). 

Model accelerations following a worst-case control loss may thus be 
high and it is clear that enormous impact loads would be generated 
if a breakaway model struck a tunnel wall. Whilst it would be 
mandatory to ensure that the tunnel itself could withstand such 
impact, extensive damage to the model would seem inevitable. Con- 
ventional large wind tunnel models represent considerable financial 
and temporal investment and it is thus certain that model control 
loss in a LMSBS cannot be permitted to be anything other than a very 
rare event. 

The prime reason for building a LMSBS is to acquire aerodynamic 
data concerning the test model. Precisely the same goal motivates 
the construction of most large wind tunnels of any type. Test 
programs in existing large tunnels tend to be interrupted relatively 
frequently by equipment malfunctions affecting data acquisition. 

This situation is tolerated since aerodynamic test programs can 
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usually accommodate delays or slippages without serious consequence 
and testing is frequently modified and new tests devised on the 
basis of early results. There seems no particular reason to assume 
that testing in a LMSBS should be any different in this respect. 

It would seem appropriate, therefore, to conceive of the "control" 
segments of a LMSBS being highly reliable, with data acquisition and 
ancillary functions remaining relatively unreliable. 

At this juncture the usage of the word "reliable" must be clarified. 
"Failures" of general systems are often classified into two types, 
these being "unreliability" where hardware mechanically fractures or 
electrically shorts etc; or "lack of integrity" where no identifiable 
physical failure occurs but nevertheless the system fails to perform 
its design task due to some other factor, perhaps adverse environ- 
mental conditions, such as electrical interference. Often the use 
of these two terms overlaps considerably and precise classification 
of particular failures is not always possible. The uncertainty over 
the terms seems especially severe in the case of MSBSs where the 
causes of loss of model control frequently do not involve simple 
mechanical or electrical failure. Herein, therefore, the reader 
should be aware of the use of the terms "reliability" and "integrity" 
somewhat interchangeably, the former usually encompassing the latter, 
though not vice versa, and implying the probabilistic likelihood of 
continuous control of the model being maintained within predetermined 
operational constraints. 

The inherent reliability of key hardware in contemporary MSBSs is 
not especially high and the large number of subsystems required leads 
to low overall reliability unless special measures are taken. 
Development of a LMSBS will involve high technical risk since at 
least the first LMSBS will be unique in overall architecture and 
most of its subsystems can be expected to incorporate significant 
technical innovations and exhibit considerable differences from any 
contemporary hardware in design, duty and environment. It is thus 
expected that considerable design effort is both essential and worth- 
while, and that the total system configuration must be especially 
formulated, in order to ensure that a LMSBS be capable of reliable 
operation. 

This Section seeks to explore the impact on overall MSBS archi- 
tecture of a requirement for highly reliable operation. Detailed 
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definition of hardware features or performance is not attempted. 

8.2 Reliability engineering in LMSBSs 
8.2.1 Modular design in LMSBSs 

The simplified architecture of a MSBS shown below suggests a 
natural division of the overall system into five major subsystems. 
There may be certain overlap in function between these subsystems, 
for instance, the system controller may include processing of model 
sensor outputs, but the division is valid insofar as the interfaces 
between subsystems are principally unidirectional and either 
fundamentally fixed or entirely optional. 



Module 

Primary interface to 
next subsystem 

Comments 

Model 

Fixed; Dynamics 

- 

Model and 
system sensors 

Optional: probably 

electronic 

May be sensitive 
to environment 

System con- 
troller 

Optional : probably 

electronic 

May include sensor 
processing 

Power supplies 

Fixed: Current, Voltage 

Closely coupled to 
•B/Ms 

E/Ms 

Fixed: Magnetic fields 

Mutual coupling 
between E/Ms 


Fig. 8.1 Inherent modularity of MSBSs 
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A highly modular design approach is thought to be particularly 
applicable in LMSBS development where the complete system will be 
novel, highly complex and composed of devices and systems from widely 
differing branches of engineering. Modularisation can be advan- 
tageously applied at the subsystem as well as the system level and 
the potential benefits of this approach include 

(1) Simplified system development. Devices or subsystems 
from several different working groups may be brought together 
in an agreed and considered overall architecture with minimum 
difficulty. 

(2) Easy upgrading of subsystems. A major advantage in 
LMSBS design since very few will be built and hence each will 
tend to be rather experimental in nature. As operational 
experience is gained and design or construction techniques 
improve, enhancements to overall performance may be made 
most economically by subsystem upgrading. 

(3) Simplified development and maintenance of subsystems. 

LMSBSs will certainly include many similar components within 
many modules or subsystems (such as several power supplies) . 
Considerable benefits accrue if groups of similar components 
can be arranged to be fundamentally identical and the number 
of different groups of components reduced. 

(4) Reduced overall cost. If (1) and (3) can be success- 
fully implemented the overall manufacturing costs will be 
reduced due to elimination of the duplication of component 
or subsystem development. 

8.2.2. Operational task classification and general reliability 
considerations 

It has already been noted that the tasks undertaken by a MSBS may 
be divided into two major classes, that is "Suspension" and "Balance". 
These tasks are of quite different criticality to the system and 
merit different approaches to their reliability. 

The suspension task is critical and failure to prevent model 
flyaway win be regarded as catastrophic failure of the system. No 
particular constraint need be placed on the quality of model suspen- 
sion for fulfilment of this task. 

The balance task is of lesser criticality and may in principle 
be unconditionally aborted at any time and with any frequency 
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without triggering catastrophic failure. The quality of sus~ 
pension should be high but useful fulfilment of this task merely 
requires that the overall system’s mean time between any event 
causing unacceptable degradation of suspension quality be reasonably 
long, with a correspondingly short mean downtime. These criteria 
should be satisfactorily met by a well-designed system without 
recourse to special measures. 

The suspension task requires that the mean time between f lyaway 
be very long indeed, perhaps in the order of thousands of hours. 
Taking into account the fact that each LMSBS will be of complex and 
largely novel design there seems no possibility of realising this 
objective without recourse to special design features. Further, if' 
component failure rates in a LMSBS follow the classical "bath-tub" 

•k 

curve then comparatively large numbers of failures would be 
expected in the early life of that system. Unless massive and 
lengthy test programs, such as those carried out for manned space 
missions, are contemplated, it would be unreasonable to expect that 
all latent design or manufacturing errors or "infant mortality" 
failures could be detected. This assertion leads directly to 
consideration of a design strategy where the system may somehow 
survive at least one random single point failure - a "fault 
tolerant" design strategy. 

Fault tolerant systems may be conveniently subdivided into two 
classes: 

(1) Fail operational 

(2) Fail soft 

The fail operational strategy implies that sufficient spare 
capacity is somehow incorporated so that the system may continue 
operation following a failure, usually at undiminished performance. 
Hardware redundancy is essential if this strategy is to be successful 
in surviving any single point failure. The fail soft strategy 
abandons the requirement for continued operation of the system at 
full performance but seeks to ensure that the consequences of any 


* 

Classical 
bathtub 
curve for 
component 


Failure 

Rate 



failure rates 
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failure are sufficiently attenuated so as to permit planned and 
graceful degradation of system performance. The requirements 
placed on LMSBSs suggest something of a hybridisation of these 
two philosophies. Critical suspension modules must fail opera- 
tional although certain degradation of performance may be per- 
missible. Balance hardware need not fail operational but must 
fail soft insofar as disruption of the suspension task must be 
prevented. Both fail operational and fail soft strategies 
generally require that propagation of failures through the system 
("cascade" failures) be prevented. Graceful degradation of 
overall system performance following a failure is quite acceptable, 
perhaps desirable, and could be manifested as an operational 
strategy where the suspension task is rapidly aborted, such as 
by initiation of a rapid shut-down sequence (RSS) , on detection 
of any failure in key hardware. It is interesting to note that 
this option is not open to, for instance, the designer of an 
advanced flight control system. In order to achieve good mean time 
between RSSs either the agglomerated mean time between failure of 
all key components and systems must be long or sufficient redundancy 
must be incorporated to permit continued operation at low risk 
following one or more component failures. 

Good serviceability requires that the fraction: 

Mean downtime following RSSs 

Mean time between RSSs 

be small. It is frequently proposed in other contexts that this 
criterion be met in complex systems by holding spare components or 
subsystems in reserve, rapidly to replace failed units. Since it 
has been anticipated that a substantial proportion of the causes of 
RSS initiation might not involve classical hardware failure, 
imposition of requirements for extensive spares holdings is inappro- 
priate here. 

It happens that even following loss of model control there will 
exist an identifiable optimum strategy, involving suspension E/Ms 
and pcr’^aps wind tunnel controls, which would "guide" the resultant 
flyaway, so as to cause minimum damage to the model and/or tunnel. 

One type of failure more serious than that described above as 
catastrophic may be contemplated, that is, a sequence of events 
triggering cascade mechanical failures of the suspension E/Ms. 
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This type of failure is technically credible, such as through helium 
supply disruption to an array of certain types of superconducting 
E/Ms, but can and should be rendered virtually inconceivable by 
careful design of relevant systems. Such a gross failure could 
otherwise be justifiably regarded as cataclysmic. 

8.3 Configuration of reliable LMSBSs by incorporation of 
hardware redundancy and some other techniques 

8.3.1 Incorporation of hardware redundancy into the E/Ms 
and some further considerations 

Introduction of some form of redundancy into the E/M array is 
necessary if sufficient force/moment capability to maintain the 
model under control is to remain following an E/M failure. However, 
overheads are incurred on several counts, including : 

(a) Provision of surplus redundant ampere-turns. An E/M 
could incorporate redundancy by being composed of several 
independent windings, each winding creating substantially the 
same useful field, configured such that the E/M is capable of 
supplying its maximum design field with one winding failed. 

The redundcint E/M must incorporate more ampere -turns that the 
equivalent non-redundant E/M and the factor : 

Total ampere-turns in redundant E/M 

imipere -turns in equivalent non-redundant E/M 
represents an E/M cost factor which depends on the number of 
independent windings per E/M. 

(b) Increased complexity in the design of the E/M array. 

This is particularly important since provision of multiple 
separate dewars (superconducting E/Ms) may seriously compromise 
close packing of E/Ms around the working section. Some segmen- 
tation of dewars is likely to be required to prevent certain 
serious failures, such as LHe supply loss, cascading through more 
than one winding. Failure to pack E/Ms closely causes larger 
E/Ms to be required since the operational effectiveness of each 
is depleted. 

(c) Lower thermal efficiency of superconducting E/Ms. A 
major path for heat leakage into the LHe environment is along the 
current supply leads to each winding where they enter the inner 
devrar. Further, since the leads cannot be superconducting, 
resistive power is dissipated along their length. Proliferation 
of supply leads will be expected to increase power losses on both 
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of the above counts . 

The effect of (a) is to decrease the E/M cost factor with 
increasing number of windings per E/M whereas (b) increases the 
cost factor with increasing number of windings per E/M. It is 
thus seen that there may be an optimum number of windings per 
E/M (Fig. 8.2). (c) demands that the number be low. 

It is thought that the E/M cost factor will be relatively small 
with the optimum number of windings per E/M and that that optimum 
number will be small due principally to three factors : 

(1) Use of Automatic Power Off (APO) . Without APO one worst 
case failure is a power supply runaway. In this event effectively 
two windings of a particular E/M are lost after one failure since 
the equivalent of one winding must simply be used to neutralise the 
runaway winding. The reduction in the idealised cost factor with 
APO varies with the number of windings per E/M as shown in Fig. 8.3. 
The system must be designed to cope with the current transient in 
the chosen APO procedure . 

(2) "Passive" Redundancy (PR). For reasons other than 
provision of redundancy it seems certain that a major architectural 
feature of any LMSBS will be the positioning of E/Ms all around the 
tunnel working section with extensive symmetry. Generally, at 
least the mirror image of any E/M suffering failure may take all or 
part of the load previously taken by that E/M. If this fact is 
exploited the complexity of each E/M may be reduced (Fig, 8.4). 

(3) State of the art of E/M dewar design. Continuous progress 
in this area is being made by many teams engaged in the design and 
construction of cryogenic systems for a wide variety of applications. 
It appears that the overall thickness of dewar and thermal insu- 
lation required to achieve satisfactorily low heat leakage into the 
LHe environment can be quite small, of the order of 1.5 inches 

(50, 51, 52). Provision of separate dewars on this basis would 
not seriously compromise the packing of E/Ms with the largest con- 
templated LMSBSs (perhaps 8 feet test sections) . Smaller systems 
would be more seriously affected, the thickness of insulation 
required not being a function of physical scale. However, full 
insulation is not necessarily required between the closest packed 
E/Ms (Fig. 8.5). 

In order to raise the fundamental reliability of each E/M winding 
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Fig. 8.3 The effect of APO on E/M cost factors (No passive 


redundancy) 
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Fig. 8.4 The concept of passive redundancy 
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Fig. 8.5 Suggested insulation requirements for close-packed E/M 
windings 
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it may be decided to incorporate redundancy into the associ ated 
ancillary systems. These are particularly numerous in the case 
of superconducting E/Ms, typically including cryogenic plant, 
external dump resistors for discharging a winding following failure, 
and instrumentation such as helium level detectors. 

It is relatively straightforward to incorporate redundant plant, 
sensors, etc., into the E/Ms but this would only be beneficial where 
the probability of, an E/M quench being initiated by means external 
to the E/M (perhaps mechanical shock) was low compared to the 
probability due to the ancillaries themselves. 

8 .3.2. Incorporation of hardware redundancy into the E/M power 
supplies 

Each winding of each E/M requires some form of power supply. 

The majority, perhaps all, of the E/Ms in a LMSBS would require 
current to be controlled with large amplitude frequency response 
to a few Hz. Terminal voltages of either sign are essential and 
it is likely that various factors would necessitate bipolar 
current capability, forcing choice of a four quadrant supply. 

Precise current and voltage requirements would depend critically 
on a system's detail design and specification but there seems no 
possibility at present of operation of a LMSBS without individual 
power supply capabilities of the order of several hundred amps and 
volts respectively. No supply is currently known to exist which 
simultaneously approaches the above requirements. 

The difficulties associated with switching powerful electrical 
supplies are severe, particularly where large quantities of 
energy (millions of Joules in Ref. 5) can be stored in each E/M. 
Large cur^rent and voltage transients might be expected during 
switching processes in this case. Therefore, whereas it is techni- 
cally feasible to incorporate redundancy in power supplies by 
straightforward duplication of supplies to each winding, with appro- 
priate means to switch to backup units following primary supply 
failure, it may be thought an unattractive possibility on practical 
grounds. The financial penalty associated with duplication of 
power supplies is, of course, a significant further consideration. 

If each E/M winding is operated with a single non-redundant power 
supply, failure of that supply necessitates immediate withdrawal 
of the winding from suspension duties since the winding current is 
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no longer controllable. With superconducting E/Ms, it is not 
necessary to power down a winding following supply failure, since 
the winding's time constant is effectively infinite. However, as 
previously mentioned, the possibility of supply runaway suggests 
immediate application of an APO procedure. The conventional means 
of achieving rapid power down with existing superconducting E/Ms 
is simply to short circuit the E/M terminals with an external dump 
resistor of suitable value. The hardware necessary for this type 
of procedure is easily duplicated, where required, to achieve 
highly reliable operation. 

It may be possible to exploit the inductive coupling that exists 
between windings, particularly between windings of the same E/M, 
which will, following Lenz's Law, tend to maintain the net magnetic 
field of any E/M as constant. Controlled shutdown of one E/M 
winding might thus tend to induce a natural increase of current in 
the other windings of that E/M, with little external intervention. 
This, of course, broadly corresponds to the required system response 
following winding failure. 

Interruptions of the power utility are not out of the question 
and some buffering would seem to be essential. Conventional means 
of buffering supply interruptions include : 

Emergency generators. 

Motor-generator sets in main power feeds. 

Batteries with static inverters. 

Considering the high transient powers involved, the use of any of 
the above techniques would place a considerable cost overhead on a 
LMSBS. However, since the power dissipation in the load E/Ms will 
be nominally zero, the power supplies may be regarded as operating 
by supplying energy to the E/Ms where current increase is required 
and retrieving it where decrease is required. It then becomes 
highly attractive from the point of view of energy efficiency to 
arrange for a quantity of energy to be somehow stored, perhaps 
within each supply, at least equal to the magnetostatic energy 
stored in the associated E/M windings at peak design currents. 

Energy is then simply transferred between the energy store and E/Ms 
with the utility demand being devoted merely to losses and 
ancillaries. As previously discussed, the supplies used with the 
UVa superconducting E/M MSBS featured energy storage (on internal 
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capacitors) in precisely the fashion described above (7), as do 
the supplies now used with SUMSBS (Appendix 1). It iS immediately 
seen tiiat supplies of this type could be made capable of executing 
a controlled RSS following loss of utility power, with only 
ancillaries such as control systems reguiring supply buffering. 

8.3.3 Reliability and hardware redundancy in the system 
controller 

The system controller is seen as the key module of any advanced 
LMSBS, since it must encompass a wide variety of functions, which 
together directly influence the configuration and operation of all 
other parts of the total system. These functions include : 

a) Model stabilisation and control. The principal traditional 
feature of MSBS control systems. lmsBSs are certain to be operated 
with tost models covering a wide range of variations in aerodynamic 
and/or magnetic configurations and characteristics. The control 
algorithms must, therefore, exhibit considerable versatility and 
robustness if they are to operate successfully, even over short time 
scales, since the aerodynamic characteristics of the test models 
will not be known in detail or to any accuracy before a test program 
commences. 

b) Prevention of model demagnetization. Permanent magnet cores 
and, to a lesser extent, magnetically soft cores, would need to be 
protected from exposure to combinations of applied fields that would 
result in partial (loss of calibration) or extensive (loss of model 
control) demagnetization of these cores. Demagnetization of 
permanent magnet cores will be possible in all realistic LMSBS 
designs and loss of magnetization of magnetically soft cores is 
credible in all MSBSs. 

c) System monitoring. The controller must incorporate 
extensive subsystem health monitoring if it is to respond intelli- 
gently to significant failures, such as by initiation of a RSS. 

It is obvious that the monitoring processes themselves should be 
carried out reliably. 

Very highly reliable system control hardware, implementing 
very dependable algorithms is seen as the foundation of a successful 
LMSBS design. 

Since it seems certain that the major part of the control system 
of a LMSBS will be realised with digital computers, there would be 
no particular difficulty in the incorporation of hardware redundancy. 
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Several options are available to achieve this, including : 

(1) Overall redundancy (Fig. 8.6) 

A small number, perhaps three, mainframe or mini- 
computers could relatively easily be configured as parallel 
independent channels such that the complete control (suspension) 
task is implemented simultaneously on each channel. Failure 
of one channel could presumably be detected by some form of 
majority voting procedure, with the failed channel being 
switched out by command signals from serviceable channels. 

This approach permits relatively straightforward execution 
in hardware, though at the cost of a considerable surplus of 
computing power (200% in a triplicated system) . 

(2) Localised parallel redundancy (Fig. 8.7) 

Series connection of several stages configured as (1) above, 
produces a configuration exhibiting improved tolerance of 
multiple failures but at the expense of increased hardware com- 
plexity and no reduction in the required surplus of computing 
power. 

(3) Multiprocessing (Fig. 8.8) 

A relatively large number of versatile processing units can 
be connected by multiple independent data buses as shown, to 
effectively form a single processing system. This system can 
exhibit high tolerance to multiple failures where the task 
previously performed by failed units is rapidly re-allocated, 
either by software or hardware actions, amongst surviving units. 
This architectural approach appears to be favoured currently 
for ultra-reliable flight control systems, and is consequently 
under intensive development for this application (53, 54) . 

The hardware and software complexity can be extreme however, 
although typically only a modest surplus of computing power 
(perhaps 50%) is required. 

The MSBS control task is broadly comparable to that of a compre- 
hensive flight control system in respect of the nature, complexity 
and frequency of the calculations required. However, the possi- 
bility of rapidly aborting the control task, such as by use of some 
form of RSS (Section 8.2.2), following detection of subsystem 
failure, limits the number of independent failures that need be 
tolerated. The development of a specialised multi-processing 
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controller specifically for this application may not, therefore, 
be justified. 

Further, whilst certain series segmentation of the control 
task appears possible, for instance : 

(a) Model state estimation, including sensor 

processing. 

(b) Control response calculation. 

(c) Translation to required power supply 

demands. 

the extent to which this segmentation is possible may be too 
limited to favour implementation of (2) above. 

The fact that the control computers are unlikely to represent 
the dominant cost in a LMSBS and since frequent upgrades or 
revisions of software are likely, particularly in the early stages 
of LMSBS development, appears to mitigate towards (1) above. 

Considerable attention will need to be paid to the software used 
in the control computers in order to ensure its own reliability, 
particularly since parallel operation of identical software is 
envisaged (parallel redundancy) . This is, however, a very difficult 
area, the science of software testing and reliability being very much 
in its infancy, and cannot be further addressed herein. 

8.3.4 Hardware and other redundancy in the model position sensors 

All model position (also velocity etc.) sensing systems so far 
considered for use with LMSBSs suffer a serious lack of integrity, 
even where the hardware used is reliable. Specifically : 

1) Optical systems. The type of system used to date with 
SUMSBS ( (11) and Appendix 1) is inoperable where rapid variations 
of the transmittivity of the optical paths are encountered. Such 
variations must be expected to occur relatively frequently in a 
cryogenic wind tunnel, for instance with free stream condensation, 
and in all tunnels when flow visualization is undertaken. 

Methods which measure in either an analogue (e.g. T.V. or laser 
scanning) or in a pseudo-digital fashion (solid state arrays) the 
absolute position of either a target affixed to, or of a suitable 
part of, the model, may be designed to tolerate the conditions 
mentioned above (24), but all optical systems must fail in conditions 
of severe degradation of the optical path. 

2) Electromagnetic position sensor (EPS). Any high intensity 
burst of electromagnetic radiation with suitable frequency content 
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will "white-out" an EPS. Such a burst might occur, for instance, 
during electrical storms or catastrophic failure of E/M windings 
or power supplies. At present, considering the problems experienced 
by the known users of EPS systems in filtering E/M power supply noise 
(10, 55) it seems unreasonable to assume that the problem of white- 
out may be easily overcome. 

It is seen that a powerful argument exists, quite separately from 
consideration of the reliability of position sensing hardware, 
favouring the simultaneous use of at least two fundamentally 
different position sensing techniques. As an example, this may 
perhaps be achieved by use of some form of optical system together 
with an EPS, since there is no particular reason to expect simul- 
taneous obscuration of the optical paths across the test section and 
electromagnetic disturbance. Without such provision a LMSBS would 
be unprotected against the common-mode loss of all position sensing 
systems arising, from the adverse environmental conditions mentioned 
above . 

Subsystem redundancy may still be attractive in cases where hard- 
ware reliability tends to be low (eg. light sources) but the inclusion 
of totally redundant sensing systems of similar type may not be 
fruitful in the context of the previous discussion, a particular case 
in point being the EPS. 

The only method by which sensible overall redundancy could be 
incorporated into the EPS is by use of systems operating at different 
frequencies since, for instance, duplicated sets of sensing coils 
would otherwise pick up the same (possibly erroneous) signals. It 
is not, however, clear whether wide spacing of frequencies is tech- 
nically feasible, the available frequency band being limited at low 
frequencies by the inherent electromagnetic noise spectra of the MSBS 
and by EPS coil and other system inductances at high frequencies. 

Now, the EPS hardware physically located within the E/M array of a 
MSBS is limited, merely an assembly of excitation and sensing coils. 
These can be made very reliable indeed. The driving, sensing and 
signal processing electronics will typically be located remotely to 
the MSBS and may relatively easily be duplicated, perhaps with fast 
switching out of failed units. There may, therefore, be no require- 
ment for overall redundancy within the EPS. 

It would seem reasonable to consider the application of special 
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aneilytic techniques to position sensor outputs in order that 
spurious information be rejected. The considerable redundancy 
of information that will exist, even with limited sensor redun- 
dancy, recommends the techniques of analytic redundancy, that is 
the detection of false sensor outputs by comparison of all 
available data with predictions derived from a real-time mathe- 
matical model. Such techniques are under intensive development 
for other purposes (56, 57). 

8.4 The impact of redundancy requirements on the sizing of the 
E/M array for a representative + configuration 

8.4.1 Theoretical background 

If each E/M is divided into a number of windings, say three, 
with the intention of tolerating the failure of one winding, then 
the E/M does not necessarily require 50% (in this case, assuming 
APO) surplus ampere-turns. The effect of passive redundancy 
(Section 8.3.1) is to reduce the requirement for surplus E/M 
capacity by distributing the load previously taken by the failed 
winding over a number of E/Ms, including, but not exclusively, the 
remaining windings of the partially failed E/M. 

It happens that the methods derived in Section 7.3 may be used 
to examine the effects of E/M failures on the maximum force and 
torque capability of specific MSBSs simply by adjusting (downwards) 
the current limit applied to a specific E/M, in order to simulate 
failure of a winding within that E/M. These methods seek to find 
the absolute optimum distributions of currents in all available (in 
this case, remaining) E/Ms so as to manifest the maximum capability 
for a particular combination of forces and torques. Here, this 
will tend to yield estimates of the minimum conceivable effect of the 
failure in question. Shortcomings in the re-distribution and 
manipulation of E/M currents will result in further performance 
degradation. 

Since there exists an infinite number of combinations of model 
attitude, position and force and torque requirements it is clearly 
impractical to attempt to present a comprehensive study of the 
effects of E/M failure, even within a single specific geometry. 

The approach taken here, therefore, is to examine the effects of E/M 
failure on some of the data from Section 7, specifically the maximxm 
force capabilities of the + configuration, with the model simply 
pitching, the resultant force vector constrained to lie in the plane 
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of pitching, and with zero magnetic torques. To recap, this 
system is of essentially similar overall configuration as SUMSBS, 
albeit geometrically and physically simplified to enable FORCE to 
be used to calculate the system's performance, and at ten times 
the true scale of SUMSBS for convenience. The absolute magnitudes 
of forces emerging from this analysis are not relevant here, but 
the relative effects of various failures, shown later, will be 
closely representative of the effects of the corresponding failures 
in SUMSBS. These studies are of limited significance insofar as 
they may convey somewhat optimistic or pessimistic assessments of 
the overall influences of the relevant failures, sideforce and 
torque generation not being studied, but serve to illustrate the 
powerful effect of passive redundancy. 

Certain simplifications have proved to be necessary in order to 
prevent the computation of each set of performance curves becoming 
excessively laborious. In fact the envelopes within which 
relevant overall performance curves in Section 7 must lie (see 
Figs. 7.13 - 7.15 and accompanying text) are omitted. 

8.4.2 Results 

The configuration here is identical to that used in Section 7, 
and is specified in Figs. 7.5 and 7.6. Fig. 8.9 defines a 
numbering sequence for the E/Ms. With the sense of the three 
force vectors of interest defined in Fig. 8.10, Fig. 8.11 shows 
the system performance with all E/Ms fully operational. This 
Figure is identical to the relevant parts of Fig. 7.13. Figs. 

8.12 to 8.19 contrast this performance to that with one-third or 
two-thirds of each E/M failed. One-third failure corresponds to 
a failure, with APO, of single windings of the main E/Ms of SUMSBS; 
two-thirds to failure, with APO, of two windings in the same E/M, 
or to worst-case failure of single windings without APO. SUMSBS' 
axial E/Ms have only dual windings, but one-third and two-thirds 
failures are computed here for consistency. 

8.5 Discussion 

The content of Figs. 8.12 - 8.19 has great significance to the 
cost of a LMSBS. With one-third of an arbitrary E/M failed, the 
reduction in force capability for the cases computed is nowhere 
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Fig. 8.10 Sense of force vectors for Figs. 8.11 - 8.19 
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Fig. 8.11 Performance of representative + system with all E/Ms 
operational 

Sense of forces as Fig. 8. 10, Geometry as Figs. 7.5, 7.6. 

This Figure shows the MAXIMUM capability for the forces as shown, 
with all other forces and torques zero. 
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Fig. 8.16/8.17 Performance of + system with failure of E/M 9 (or 10) 
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Fig. 8.18/8.19 Performance of + system with failure of E/M 2 (or 4,6,8) 
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greater than one-sixth and is typically around half that value. 
Should this data prove to be representative of all cases, it would 
follow that the surplus E/M capacity required to permit tolerance 
of single arbitrary winding failures with triple E/M windings is 
not greater than 20%. This should be contrasted to the value of 
50% shown in Fig. 8.3, where passive redundancy was not considered. 
Dependent upon which parts of the performance curves, such as Fig. 
8.11, prove to be critical, the requirement for surplus capacity 
may be even less than stated above. 

It is currently thought that the further cost overheads incurred 
due to the extra complexity of, say, triple windings, as compared 
to single windings of similar total ampere-turns, will prove to be 
modest and acceptable. Further, it is seen that dual windings may 
be viable, incurring a basic cost factor increment of not more than 
33%, and that large numbers of windings per E/M are unlikely to be 
worthwhile (see Fig. 8.20). 

Whilst no corroborative computed data is available, it is easily 
seen that increases in the number of individual E/Ms in a MSBS will 
tend to increase the effect of passive redundancy, thus decreasing 
the cost factors for given numbers of windings per E/M, since there 
would then exist greater opportunity for redistribution of the load 
of a failed winding. The number of E/Ms in SUMSBS is 10 and is 
thought to represent a minimum conceivable number for LMSBS designs. 
LMSBS cost factors are thus likely to be lower than those quoted 
above for comparable cases. 

Briefly to summarise this Section, it has been shown that the 
increase in cost, complexity and sophistication of LMSBS hardware 
caused by a requirement for highly reliable operation need not be 
outrageous. E/M cost factors can be kept satisfactorily low with 
apparently modest hardware complexity, appropriate use of existing 
types of position sensor should provide adequate redundancy and 
independence of information, comparatively straightforward control 
system hardware may be used and the overall power consumption of the 
complete system need not be materially affected. 
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Fig. 8.20 The effect of passive redundancy (PR) and APO on E/M 
cost factors, PR factor as Figs. 8.12-8.19 
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9„ DISCUSSION AND CONCLUSIONS 


Since the content of this thesis divides into three principal 
subject areas (roll control, extreme attitude testing and reliability), 
the first part of this Section is divided likewise. These sub- 
sections will commence with a brief summary of material from the main 
text, followed by specific discussion and conclusions. 

9..1 Roll control 
9 1 . 1 Summary 

A new method of roll control has been proposed, based on the 
mcignetization, in the spanwise sense and with port-starboard symmetry, 
of the wings of winged aircraft models. The system exhibits two 
principal variations, since the spanwise magnetizations may be 
permanent (SPMs) or induced (SIMs) . 

The roll torque generated by SPMs in a suitable uniform applied 
field gradient is shown to be proportional, at low applied field 
levels, to the first moment of core volume about the model's plane 
of symmetry, also, of course, to the applied field strength and the 
intensity of the spanwise magnetizations. The maximum torque capa- 
bility of SPMs is of major interest and is expected to depend princi- 
pally on the magnetic properties of the core material (almost 
certainly ReCo) but, due to the lack of available data, has proved 
impossible to identify at this stage. However a crude estimate 
indicates that substantial torques may be available. 

The existence of cross couplings between model degrees of freedom 
is recognized but it is shown that with unswept wings there exists 
only one primary coupling and relatively few other couplings of 
significance. 

A pilot model with SPM wings has been suspended in full 6-component 
control (with an all-digital control system, in fact) with no 
difficulty. The demonstrated roll torque capability and the achieved 
levels of stiffness and damping of rolling motion were adequate to 
recommend the system as a first choice for future use in SUMSBS. 

Subject to the provision of adequate E/M power, the SIM system will 
provide greater roll torque capability than SPMs, due partly to the 
higher levels of spanwise magnetization that are attainable. No full 
demonstreitions of the system were possible (that is, with 6-component 
model suspension) due to the unsuitability of the available E/M 
configuration. However, comprehensive numerical computations (GFUN) 
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clearly indicate that high torques will be available, restricted 
principally by current E/M technology. Limited experimental 
torque measurements at low applied field levels exhibit excellent 
agreement with corresponding calculations. 

9.1.2 Discussion and conclusions 

It is thought that development of a powerful method of active 
roll control is absolutely essential before construction of a LMSBS 
could be attempted. 

The variant of the Spanwise Magnet scheme that could be fully 
demonstrated (SPMs) may be practical for LMSBS application, subject 
to accurate identification of adequate maximum torque capabilities. 
Some complication of model design is necessary to incorporate 
permanent magnet cores into load bearing wing structure and, princi- 
pally for this reason, the system is regarded with less favour than 
the alternative variant of the Spanwise Magnet scheme, that is, SIMs. 

The limited practical experience so far gained with SPMs indicates 
that the system is an excellent choice for SUMSBS and demonstrates 
the general practicality of Spanwise Magnet roll control (SPMs or 
SIMs) . Cross couplings between model degrees of freedom due to 
the presence of the spanwise magnetizations are easily handled by 
suitable control algorithms and present no serious difficulty. 
Suspension quality with the SPM winged model was good. The spanwise 
magnetizations and the applied "roll" and other fields do alter 
(mostly adversely) the natural magnetic stiffness in several degrees 
of freedom but modest adjustments to the controller will be sufficient 
to restore the full stiffness and high stability currently achieved 
with unwinged models. Further computations of these effects would 
be appropriate, though the data presented in Section 4.3.1 contains 
the essential elements for a preliminary examination. 

The SIM system is currently viewed as the most promising system 
for LMSBS application. Model design is very straightforward, since 
wings may be machined from the solid with suitably chosen materials 
and the predicted torque capabilities are very high, limited princi- 
pally by current E/M technology. 

The level of confidence placed on the available calculations, of 
crucial importance of course, is currently high. However, extensive 
further and alternative calculations and a comprehensive demonstration 
at small scale are seen as essential before a LMSBS is committed to 
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the system. Such verification would necessarily include treatment 
of cross coupling effects, which are ejg)ected to be somewhat more 
pronounced than with SPMs. 

Subject to such computation and verification being satisfactory 
it is believed that careful study and optimisation of the roll E/M 
array will prove fruitful. The configurations used herein were 
chosen to be representative of likely LMSBS E/M arrays but were 
dimensioned entirely arbitrarily and are clearly deficient in many 
respects. If the ratio of useful field at, say, the model wing tip 
to the peak field level within or around the individual E/Ms is taken 
as a measure of performance, then the performance of the chosen SIM 
arrays is low. Adjustments to the E/M sizes, locations and winding 
cross-sections are expected to yield significant improvements. 

9.2 Extreme attitude testing 
9 . 2 „ 1 Summary 

The largely intuitive realisation that the suspended model is 
influenced merely by the local applied magnetic fields and is other- 
wise unaffected, apart from the direction of the gravity vector, by 
its location or orientation within the suspension system, is easily 
illustrated by analysis. 

This realisation leads directly to the inherent ability of MSBSs 
to suspend and control models at any attitude and location, provided 
that fields and field gradients in the required senses and magnitudes 
local to the model can be somehow developed. 

It happens that the transformations between applied fields in 
balance axes and those seen in model axes are straightforward, but 
corresponding couplings for field gradients are relatively complex. 
Further difficulties are presented by the effects of the spatial 
variations of the applied fields and field gradients about the model's 
centroid. These are quite pronounced in systems such as SUMSBS 
where the E/M faces are in relatively close proximity to the model, 
but will tend to be of lesser magnitude in LMSBSs. Such effects 
render the adequate separation of model degrees of freedom, essential 
for the operation of contemporary control algorithms, quite complex. 

The prediction that SUMSBS would be usable, in the + configuration, 
from nominally 0° to 60° angle of attack (5-component control) was 
verified by achieving suspension of an axisymmetric model over this 
range. Th«j inability of SUMSBS to operate around 90° angle of attack 
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in the + configuration and the difficulties apparent with the X 
configuration were traced to straightforward deficiencies in the 
E/M array . 

9.2.2 Discussion and conclusions 

It has long been claimed that one major advantage of wind tunnel 
MSBSs was the possibility of rapidly selecting near-arbitrary model 
attitudes. However, it is believed that no analytic or practical 
verification had previously been attempted. 

The problems presented by the gross changes in magnetic coupling 
between E/Ms and model, occurring during large changes of model 
attitude, are difficult but surmountable, the couplings being, of 
course, analytic or near-analytic and repeatable. 

Digital control is essential for sensible operation over wide ranges 
of angle of attack in order to permit adaption of the control loops 
to the changes in magnetic couplings and full digital control was, in 
fact, used for all demonstrations. 

Suspension quality with SUMSBS at the higher angles of attack was 
poor due to three major factors: 

(i) The decoupling of the axial position sensor was 
inadequate, resulting in coupling between pitch, heave 
and axial motions. 

(ii) The already approximate computed magnetic couplings 
(using FORCE) , essential for adequate decoupling of model 
degrees of freedom in the controller, were further simplified 
before incorporation into that controller, for convenience. 

(iii) The natural magnetic stiffness of certain model 
degrees of freedom in the quasi-steady applied fields 
supporting the model's deadweight apparently became strongly 
negative at the higher angles of attack. Analysis and com- 
putation of this effect should be undertaken before further 
attempts at high angle of attack suspension are made. 

The inability of SUMSBS to operate above about 60° angle of attack 
in the configuration and the poor performance predicted for the X 
configuration indicate that this simple fully symmetric 10 E/M array 
is unlikely to represent a good choice for a LMSBS. However, the 
shortcomings of the configuration, identified theoretically, are 
easily rectified by augmentation of the E/M array, indeed such 
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augmentation would be likely if SIM or SPM roll control was employed, 
and it is strongly believed that unlimited model attitudes, with full 
6-component control, will be easily achievable in realistic LMSBS 
configurations. The most difficult problem in such a system may, 
indeed, prove to be the provision of model position and other sensors 
capable of satisfactorily monitoring wide ranges of model attitudes 
and motion. 

9.3 Reliability 

9.3.1 Summary 

The widespread concern that the capital cost overhead incurred by 
a requirement for highly reliable operation would drive already 
expensive LMSBSs out of reach, has been partly dispelled. The two 
areas of greatest expenditure, that is the E/Ms (with associated 
plant) and the E/M power supplies, need only carry modest levels of 
redundancy and the complication of design, such as segmentation of 
E/Ms into several separate windings, is not severe. Indeed, 
segmentation of E/Ms and corresponding duplication of power supplies 
was featured in the only comprehensive LMSBS design study made to 
date (5), for reasons apparently unconnected with reliability. 

It is shown that the reliability of the remaining MSBS subsystems 
such as model position sensors and control systems may be adequately 
raised, albeit at the expense of some design complication, such as 
by the incorporation of at least two fundamentally different position 
sensing systems. 

9.3.2 Discussion and conclusions 

An unreliable LMSBS would be an expensive burden rather than a 
benefit. Any attempts to argue that a relatively high attrition rate 
of models in a "cheap" LMSBS could be justified is undermined by 
simple application of Murphy *s Law, which indicates that the most 
expensive and intricate models would be destroyed before any useful 
aerodynamic data had been obtained. 

Though full technical presentation would be inappropriate here, 
preliminary demonstrations of hardware (E/M and E/M power supply) 
redundancy have already been made with SUMSBS. Specifically, with 
an axisymmetric model in steady suspension (wind-off and low angles 
of attack), ANY power supply, including those contributing towards 
the support of the model’s weight, may be switched off with impunity. 
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The transient motion following the switch-off (or switch-on) is 
naturally controlled by the existing control algorithms, the model's 
weight, for instance, simply being redistributed amongst remaining 
operational E/Ms, and motion could be virtually eliminated if the 
controller was able to detect the "failure" of the power supply in 
real time, such as by monitoring the E/M currents. These demon- 
strations clearly confirm the existence and the powerful effect of 
Passive Redundancy, in this case not depending in any way on special 
configuration or adaption of a relatively simple MSBS. 

The key to reliable and economical LMSBSs is the EXPLOITATION 
rather than the SUPPRESSION of the many unique features of such 
systems . 

9.4 SUMSBS 

The commissioning of the new features, hardware or otherwise, of 
SUMSBS, has proved remarkably trouble free. 

Incorporation of digital control hardware and digitized control 
algorithms has been straightforward, with no particular difficulty 
in interfacing with the rest of the system. This success owes much 
to the care taken in the hardware design phase to create a system of 
high integrity, such as by the incorporation of complete electrical 
isolation between the control computer and the MSBS rig itself and 
the incorporation of sensible analogue filtering in position sensor 
preamplifiers and power supplies. The general performance of the 
controller is now as good as had ever been achieved previously with 
analogue systems. Changes or adjustments of control parameters may 
be made with precision and repeatability in seconds, rather than 
minutes or hours, and there now seems no reason to consider anything 
other than full digital control for a LMSBS or, indeed, any new MSBS. 

Considerable scope for improvements to the performance of the 
controller certainly remains, including the quality of suspension 
achievable with wind-on, the adaption of the controller to changing 
test conditions and the performance at extreme model attitudes. 

The new transistor switching power supplies, radically different 
from anything known to have been used previously in this application, 
are steadily amassing many hours of operation, to date completely 
faultlessly. The E/M loads (high inductance, low resistance) appear 
to be well suited to this type of supply- The level of radiated 
electrical noise, due to the high frequency switching, is much less 
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than had been feared, and has presented no serious difficulty. 
Accurate monitoring of E/M currents is not straightforward, but 
can be achieved with appropriate techniques. The independent 
and bipolar operation of each E/M has considerably enhanced the 
performance and versatility of the system, with no significant 
disadvantage yet apparent. 

The principal weakness in SUMSBS remains the position sensors. 
Significant changes in model geometry or gross changes of model 
attitude are only possible with mechanical realignment of the 
sensors. This is a serious inconvenience but in no way indicates 
a fundamental difficulty. Superior position sensing systems, 
optical or otherwise, are available and may be gradually incor- 
porated into SUMSBS in future years. 
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APPENDIX 1 


THE SOUTHAMPTON UNIVERSITY MAGNETIC SUSPENSION AND BALANCE SYSTEM 
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Al.l Ifitroduction 

The present status of SUMSBS hardware, represents only an interim 
phase in a more or less continuous program of development. This 
Appendix seeks merely to acquaint the reader with the current major 
features of SUMSBS, particularly those relatively recently intro- 
duced or updated. Fuller descriptions of certain hardware may be 
found in some earlier material, such as Refs. 11, 15. 

A1.2 E/M g eometry and configuration , 

The approximate dimensions of the main E/Ms and associated iron 
cores are shown in Figs. Al.l/2. All 8 E/Ms are nominally identical, 

6 in fact being carried over from earlier versions of SUMSBS. 
Previously, the iron cores had been configured rather as yokes (Fig. 
A1.2) with one "upstream" and one "downstream" E/M per assembly. 

The cores have , been modified by removal of the centre portions (shaded 
in Fig. A1.2) in order to reduce the fore-and-aft inductive couplings 
between associated E/Ms. This coupling had not previously been unduly 
troublesome, in fact the presence of the link augmented, for instance, 
the conventional "lift" field by some 30% due to the provision of 
return flux paths, but was considered undesirable for suspension at 
high angles of attack, or with spanwise magnet roll control, where 
more or less independent control of individual E/Ms is mandatory. 

Each E/M is wound with approximately 400 turns of 12 gauge enamelled 
copper wire, uncooled, and split into three sections normally 
connected in series. 
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The twin axial E/Ms (Fig. A1.3) are somewhat similar in con- 
struction but are wound with thinner wire (14 gauge) split into 
two sections of approximately 500 turns each, normally connected 
in parallel. 

The E/M configuration is fully symmetric and is illustrated 
in Figs. A1.4 - A1.8. 

A1.3 Power Supplies 

All 10 E/Ms are currently independently controlled by CSR 
Contraves NC 421 servo controllers (Ref. 58 and Figs. A1.9 - Al.lO). 
These controllers, rated at 21A continuous, 45A peak and lOOV for 
four quadrant operation, were, in fact, designed for industrial 
D.C. servomotor drives but are well suited to this application. 

The fundamental mode of operation is pulse width modulation, with 
the ability to be driven unconditionally to 0% or 100% modulation. 
Since the switching frequency is nominally 5kHz the forward path 
voltage gain, neglecting the effects of current feedback, can be 
set flat to around IkHz, representing a considerable advance over 
previous SUMSBS supplies (11). Fixed internal current feedback is 
in fact provided (Fig. Al.ll) and i.s considered to be extremely 
advantageous. Synthesis of model control algorithms is eased, 
countering the lengthening of the E/M time constants caused by 
removal of the series ballast resistors previously used (Section 
2.2.2 and Ref. 11), the linearity of response is improved, output 
drift with temperature and time reduced and the effects of voltages 
induced in E/Ms by magnetic coupling from others is significantly 
attenuated. 

The power regulators which feed the controllers (Figs. A1.12, 
A1.13) themselves feature internal energy storage (on capacitors) 
totalling some 280 Joules overall, this being slightly greater than 
the total magnetostatic energy retained in the E/M array with all 
E/Ms at peak continuous current. Despite the fact that only some 
20% of this capacity is normally dynamically available, it is seen 
that since regeneration of load energy is allowed (four quadrant 
operation) the net power consumption of the complete system under 
normal conditions will tend to be low, being mostly due to resistive 
losses in E/Ms and power stages. Some 25 kW (continuous rating) 
of power is available from transformers and has been found to be 
more than adequate for all requirements. 
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Fig. A1.7 Schematic diagram of 


SUMSBS' E/.M array. 
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A1.4 Position sensors 

The basis of the position sensing system used currently is an 
array of four analogue detectors of traditional form (Ref. 11 and 
Pigs. Al.14 - A1.21). Certain detail improvements have been made, 
most significantly to the light pickup assemblies where the first 
stage of amplification and filtering is now incorporated directly 
behind the sensor (Fig. A1.18). This is done in order to reduce 
the pickup of airborne electrical noise along the cabling between 
the MSBS rig and the main control system, this cabling now carrying 
relatively high level signals. 

High angle of attack capability is effected simply by appropriately 
rotating the aforementioned sensors broadly in the pitch sense (Figs. 
A1.21, A1.22). The range of sensitivity is such that sensing of a 
range of angle of attack from O to 60 degrees is practicable (see 
Section 7.8). 

Two sensors differ significantly from those described above, namely 
the axial position and roll attitude sensors. The axial position 
sensor was required to be capable of monitoring relatively large 
apparent axial motions, these being induced during large pitch 
rotations in the chosen configuration (Figs. A1.20, A1.21) . A long, 
slender photodiode provides the necessary pickup range but would tend 
to be unduly sensitive to ambient light unless further measures are 
taken. Utilisation of a laser light source (ImW He - Ne) enables a 
narrow band optical filter of suitable wavelength to be installed in 
front of the photodiode, reducing the ambient light pickup by a factor 
of approximately 50 (measured) . 

The roll attitude sensor comprises a laser (again ImW He - Ne) , 
directed onto a mirror suitably located on the model, the reflected 
beam falling onto a UDT PIN-SC/50 two axis position sensing detector 
(Refs. 60, 61, and Figs. A1.20, A1.21), This device is essentially 
analogue in nature, providing four output signals which, when suitably 
processed, yield information as to the location of the centroid of 
the total incident light on the device's surface. The output can be 
made relatively insensitive to ambient light, partly by use of optical 
filters, and to variations in the intensity of the signal beam. The 
two axis capability of the device may be more fully exploited if 6- 
component high angle of attack suspension is attempted. 
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Detail C 





rotated in sense D 
for high angle of 
attack operation, 
(see Fig. 7.26) 














A1.5 Models 


Recent recommissioning of SUMSBS has been performed almost 
entirely with simple cylindrical models of arbitrary aerodynamic 
shape (Fig. A1.23) fitted with existing cylindrical Alnico V 
magnetic cores. It has been found that these cores can now be 
demagnetized following model flyaway by certain combinations of 
applied fields, necessitating automatic shut-off of all E/Ms 
whenever position sensor signals indicate that the model is departing 
from their range. Even with this protection, the long term cali- 
bration of such a core must be regarded as dubious, recommending 
future changes in the core material. 

Exploratory 6-component suspension with spanwise magnet roll 
control has been achieved with a model basically as described above, 
but crudely adapted to carry arbitrarily dimensioned Alnico V wing 
cores (Fig. A1.24), The initial configuration of the roll angle 
sensor (Figs. A1.20, A1.21) necessitated a somewhat awkward location 
of the model-mounted mirror but this is, of course, not a fundamental 
shortcoming of this general technique. 

A1.6 E/M current sensing 

Since the output voltage from the CSR power supplies appears as a 
pulse width modulated, 5kHz, lOO V peak-to-peak square wave, with 
neither output terminal near ground potential (Fig. A1.25), special 
measures are necessary to acquire accurate signals from conventional 
shunts in the E/M supply leads. A shunt signal typically consists, 
under these conditions, of a small (- IV) differential component 
representing load current, together with a common-mode AC signal of 
approximately 100 V peak-to-peak (Fig. Al.25). Fully isolated signal 
amplifiers are necessary to permit safe handling of the peak common- 
mode voltages and high rejection of this signal component is essential 
for accurate resolution of load current. 

Configuring dual series shunts as signal sources for a linear 
isolation amplifier (Burr-Brown 3650 HG, Ref, 62), as shown in Fig, 
A1.2G, causes the "common-mode" signal to appear, in fact, as an 
"isolation" signal, that is, wholly between the input common and true 
ground. Adequate rejection capability of this signal is available 
from the chosen class of device, specifically some 140 dB at D.C. 
falling to approximately 90 dB at 5kHz, all figures given for a device 
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voltage gain of lO, from Ref. 62. In fact, the electrically noisy 
environment in which the devices are operated causes some degradation 
of this performance, mostly with regard to output noise, but modest 
filtering of the output provides a satisfactorily clean and effec- 
tively real-time signal for data acquisition purposes. 

Al.7 Control systems 

The control loop is implemented wholly on a digital computer 
(PDF - 11/34) with only limited signal conditioning performed 
externally. Local 12 bit A/D and D/A converters communicate with 
the host computer via two fully isolated, high speed, uni-directional 
16 (+ control) bit data buses, the computer in fact being located 
some 90 feet from SUMSBS itself (Figs. Al.27 - A1.28). Treatment of 
the control algorithms used is outside the scope of this report but 
all suspension necessary for the work reported herein has been accom- 
plished with relatively simple digital versions (59) of traditional 
analogue algorithms (11), absolute performance of such controllers 
thus being comparable to that achieved by the previous analogue 
systems. 
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Fig. A1.27 Schematic diagram of digital control system hardware 
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A/D subsystem 


Fig. A1.28 Schematic diagrams of A/D and D/A 
converter subsystems. 
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THE PROGRAM "FORCE" 
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A2.1 Introduction 

The forces and torques acting on isolated bodies whose 
magnetization distribution is known, are amenable to approximate 
numerical calculation where the applied field can be somehow found. 

The program FORCE is an example of one simple method, in fact 
carrying out elementary numerical integration of the relevant 
vector products (Eqns. 3.1) over the volume of the core. 

The core is converted to an array of dipoles, each dipole 
representing the magnetization of suitably chosen model elements. 

Such representation of the core is only straightforward if the core 
magnetization is essentially uniform, whence each dipole can be 
placed at the centroid of the corresponding element. Uniform 
magnetization is frequently assumed in the case of ReCo materials, 
where the material's high coercive force precludes self-demagnet- 
ization effects. However, computations may be approximately 
applied to other cases, such as slender axially magnetized iron cores. 

Computation of applied E/M fields is generally only straight- 
forward where the E/Ms are air-cored. In the case of MSBSs where 
the E/Ms are relatively remote from the model it is usually adequate 
to represent the E/Ms as assemblies of straight line wire elements, 
whence the external field is found by simple summation of the contri- 
bution of each element. The field of individual elements is 
easily found, such as by use of the Biot-Savart law (63 and various) . 

A program segment that had been shown to reliably compute the field 

from straight line wire elements was extracted from the program 
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TABLE, courtesy Massachusetts Institute of Technology (64, 65), 
and is used in FORCE essentially unmodified. 

A2 . 2 Program structure 

The fundamental program structure is shown in Fig. A2.1. The 
program is heavily modularised to permit implementation on a 
minicomputer (PDF 11-34) with relatively limited available memory. 
Main model and E/M data storage is in virtual arrays. Memory 
requirements are comfortably below the available 28K words (main 
code) and 32K words (virtual arrays) with 1000 model and 1000 E/M 
elements (see Fig. A2.2). 

The main (root) segment of the program simply performs the 
function of calling the operational subroutines in sequence, this 
approach allowing these subroutines to be overlaid in memory, as 
Fig. A2.2. 


A2.3 Subroutine descriptions 

MOD IN controls the filling of the model core dipole array via calls 
to user named datafiles or the subroutines FUSIN or WING. Model 
core data can be output to the console device or to user named data- 
files. The model can be located arbitrarily in space and looping 
can be called whereby one locating parameter, say pitch attitude, is 
stepped by chosen increments over a chosen range, with force and torque 
calculations repeated at each step. 

The fundamental equations controlling the creation of the dipole 
array are as follows; 

* 

M - / M 6v where the integration is over 

volume of a single element. 


leading to modified forms of Eqns. 3.1 : 



Z (M* . V ) H. 
1 —1 — 1 

* 

Em, X h, 

1 —1 —1 


(A2.1) 


FUSIN permits representation of simple axisymmetric fuselage geometries 
as shown in Fig A2.3 Where FORCE is used simply to calculate fields 
due to particular E/M arrays, FUSIN fills the CORE array with dummy 
model elements of zero magnetization (Fig. A2.4) . 
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Fig. A2.1 Fundamental block structure of program FORCE 
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Fig. A2.2 Memory overlay structure of program FORCE 
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Sectioii on AA 


Section on 


As shown. 


— showing elemer 
NL = 5, NA = 4, NR = 


Fig, A2.3 Axisynunetric fuselage 








Fig, A2.4 Field grid. Fuselage input option = 2. Zero magnetization 
in all elements 
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WIMG permits representation of simple slab wings, magnetized 
symmetrically (SPMs, see Section 4), as shown in Fig. A2.5 

COIL creates an array of straight line wire elements representing 
the assemhily of E/Ms. Symmetry is used extensively to reduce the 
quantity of user input per program run and the program's main 
memory requirement. Data is input for a single "master" E/M, 
which can later be reflected and/or rotated in certain ways to 
generate other E/Ms in that particular "set" (Pig. A2.6). Up to 
five E/M sets may be input sequentially, permitting analysis of 
relatively complex MSBS configurations. Four types of E/M may be 
specified, single wire loops (Fig. A2.7), E/Ms effectively composed 
of straight blocks of conductor of quadrilateral cross section 
(Fig. A2.8) and "lateral" or "axial" pseudo-circular E/Ms (Figs. 

A2.9, A2.10). Pseudo-circular E/Ms are reduced to assemblies of 
circumferential straight line elements, positioned at suitable 
radii so as to create a field at the centre of each loop identical 
to that created by the equivalent circular loop carrying the same 
current. E/M data may be output to user named datafiles or to the 
console device. Data for a complete array (several E/M sets where 
required) may be input from a user named datafile. 

E/M ampere-turn products are input individually at the console 
device in a standard sequence (Fig. A2.11). 

FACE generates the co-ordinates of straight line wire elements 
distributed over the cross section of a block of conductor, per- 
mitting E/Ms to be input as a series of co-ordinates representing 
the interfaces where blocks of conductor join as Fig. A2.8. 

FACE is derived from the TABLE subroutine COOKS (MIT, (65) ) , 
but differs in two important respects. Firstly, since COOKS only 
generated wire element distributions that correctly simulated 
uniform current density over the winding cross section when the 
interfaces were quadrilateral, FACE generates the fourth corner of 
each interface internally. Secondly, COOKS used the input inter- 
face co-ordinates to define directly the location of the outermost 
wire elements. FACE uses the interface co-ordinates to describe 
an outer envelope, inside which all elements are distributed. 

Thus, identical input data to COOKS and FACE will result in different 
wire element distributions. 
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Fig. A2.6 Symmetry options for E/M sets 



As shown NINT=4 
(no. elements) 


Endpoints input 
in sequence thus 

xi ,yi ,Zi ,X2 , 

74 , Z4 


Current flows 
in sense shown 


Fig. A2.7 Single wire loop E/Ms. E/M input option = 1 


Detail of interface 1 
B MM A 



Pig. A2.8 Interface array E/M input. E/M input option = 2 
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AXIAL 




X2 


Z 



NDIVY or 
NDIVX 


NDIVR 


Section on AA showing 
wire loop distribution 

As shown : 

NDIVR=1, NSEG=12 

NDIVY (lateral) =2 

NDIVX (axial) =2 



AXIAL and LATERAL 


Fig. A2.9 and A2.10 Lateral and Axial Pseudo-circular E/Ms. E/M input 

option = 3,4. 
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ISYMM = 4 (4E/MS) 


(14) 



ISYMM = 6 (8 E/Ms) 



Fig. A2.11 Standard sequence of E/M ampere-turn input 
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CALC controls the calculation of E/M fields, via calls to the 
subroutine FIELD, summation of model forces and torques and the 
routing of calculated data. The position and orientation of 
model elements is updated using previously input data (Fig. A2.12). 
FIELD then provides the field and field gradient components at 
each model element in turn, the relevant vector products of field 
and magnetization being calculated and summed. Model forces and 
torques are expressed about the model's own co-ordinate origin, 
in translated but unrotated axes (Pig. A2.13). The calculated 
field at eeich model element may be output to datafile or the 
console device. Force and torque data may be routed likewise. 
Calculations are repeated if stepping of the model position or 
orientation had been previously selected. 

FIEL D calculates the field and field gradient components at the 
single points required, due to the specified E/M array. Individual 
E/M set calculations are made first for the master E/M, then for 
E/Ms generated by the chosen reflection/rotation symmetries. Set 
calculations are repeated for all specified sets with the final 
field data being the sum over all E/Ms of all sets. 
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A2.4 


GE:NEF<AL INSIfi'UCTIONB FOE' USERS OF PROGRAM FORGE ************t*tm**t*t 


THE MAXIMUM MUMBER OF MODEL ELEMENTS CURRENTLY PERMITTED IS 1000 , 

IHE NUMBER OF ELEMENTS REQUIRED FOR A PARTICULAR MODEL MAY BE CALCULATED 

AS follows; 

N=(NOX+D)H<NOY+l))i((NOZ+l) FIELD GRID 

=Nl..)i<NA)KNR CYL. FUSE. 

+NAKNR)l'.NAl)t:2 ELLIP. ENDCAPS 

+NS!«NC)KNT)tt2 SLAB WINGS 


THE MAXIMUM NUMBER OF E/M WIRE ELEMENTS CURRENTLY PERMITTED IS 1000 
THE NUMBER OF ELEMENTS IN A PARTICULAR E/M ARRAY MAY BE CALCULATED 


AS follows; 

M=SUM 


OOER ALL SETS OF! 
NINT 

NINT*MM)KNN 

ND I VRJiiND I V Y)KNSEG 

NDIORKNDIOXJicNSEG 


SINGLE LOOP 

INTERFACE ARRAY 

LATERAL CIRCULAR 

AXIAL CIRCULAR 


ONLY THE MASTER E/M IN EACH BET AFFECTS THE OALUE OF M. E/MS GENERATED 
BY SYMMETRY OPTIONS DO NOT REQUIRE FURTHUR STORAGE ALLOCATION. 


THE NUMBER OF MODEL AND E/M ELEMENTS SHOULD ALWAYS BE CHOSEN TO BE AS 
SMAI I AS WILL YIELD ACCEPTABLE RESULTS. CHOICl: OF UNNECESSARILY LARGE 
NUMBERS OF ELEMENTS WILL RESULT IN INCONUENIENTLY LONG EXECUTION TIMES 
AND PERHAfS IN THE ACCUMULATION OF ROUNDING ERRORS. 

COMPUTAION TIMEi. TAKEN FROM THE COMPLETION OF THE LAST USER 
INPUT TO FINAL PROGRAM COMPLETION^ NEGLECTING HARDCOPY PRINT TIMESr 
HAS BEEN FOUND TO BE APPROXIMATELY; 


T=N)CMM*<STEPS+1 )*0.01 

WHERE MM=-TOTAL NUMBER OF E/M ELEMENTS BUMMED ODER ALL 
E/MS OF ALL SETS. 

PDP 11-34 WITH HARDWARE FLOATING POINT OPTIONf 
RT-ll OPERATING SYSTEMy GIOEN IN SECONDS. 


THE. STANDARD FORMAT FOR THE TITLES OF USER NAMED DATAFILES UNDER RT-ll 

IS AS follows; 

**)!(**)*!. DAT 

WHERE IS AN ARBITRARY 6 CHARACTER ALPHANUMERIC NAME 


ALL INPUT/OUTPUT IS IN SI UNITSy SPECIFICALLY; 

DIMENSIONS IN METRES 
FIELDS IN TESLA 

FIELD GRADIENTS IN TEGLA/METRE 

ANGLES IN DECREES (CONVERTED INTERNALLY TO RADIANS) 
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A2.5 


USER INPUT UARIABLES FOR PROGRAM FORCE ******* HoK *** He 

REFER TO FIGS.A2.1 TO A2.13 FOR DEFINITION OF VARIABLES 
NOT CITED BELOW 

MODEL GE.OMETRY OPTION <MOPT) 

AB where: 

A-WING INPUT OPTION 

0 NO WING 

1 SLAB WING 
B=FUSELAGE INPUT OPTION 

0 FROM DATAFILE (WHOLE MODEL) 

1 CYLINDRICAL FUSELAGE 

2 FIELD GRID 

3 AS liELLIPSOID ENDCAPS 

EG. MOPT"1.1. SPECIFIES CYLINDRICAL FUSE. WITH SLAB WINGS 
MOPT=l SPECIFIES CYLINDRICAL FUSE. ONLY 

MODEL OUTPUT OPTION (lOUTOP) 

0 NONE 

;l TO CONSOLE 

2 TO DATAFILE 

MODEL DISPLACEMENT OPTION (MDISPO) 

0 NONE 

1 SINGLE SPECIFIED DISPLACEMENT 

2 LOOPING OF CHOSEN PARAMETER FROM SPECIFIED INITIAL 

DISPLACEMENT 

~2 AS 2 BUT FROM NULL INITIAL DISPLACEMENT 

WHERE MDISP0--2»-2 

STEPS NO. OF INCREMENTS IN CHOSEN PARAMETER NOT INCLUDING 
START POSITION 

FARM 1-6 SELECTING X»Yr/i> YAWfRITCH/ROLL DISPL. RESPECTIVELY 
SIZE SIZE OF STEPS (METRES OR DEGREES) 

NO. OF E/M SETS (NSET) 

NO. OF INDEPENDANTLY SPECIFIED SETS OF E/MS 

0 INPUT COMPLETE. ARRAY FROM DATAFILE 

(E/M) INPUT OPTION (INPOPT) 

1 SINGLE LOOP 

2 Interface array 

3 LATERAL PSEUDO-CIRCULAR 

4 AXIAL PSEUDO-CIRCULAR 

(E/M) SYMMETRY (ISYMM(IfJ)) 

1 NONE (MASTER E/M ONLY) 

2 REFLECTION IN Y»Z AND X»2 PLANES 

3 1 ROTATION ABOUT X AXIS - BY 180 DEGREES 


4 3 ROTATIONS 90 

5 5 60 

6 7 45 

7 REFLECT IN Y»Z PLANE 

8 REFLECT IN XrY Xi-Z AND Yi-Z PLANES 

9 AS 3 WITH REFLECTION IN Y»Z PLANE 

10 AS 4 YfZ 

11 AS 5 Y»Z 

12 AS 6 YrZ 


E/M OUTPUT (lOUTOP) 

0 NONE 


- 276 



1 interface: array coordinates only 

2 AS 1+WIRE flLEMENT COORDS. 

3 WIRE ELEMENTS TO DATAFILE 

E/M AMPERE TURNS 

TOTAL AMF'ERE- TURM PRODUCT FOR EACH E/M. INPUT IN STANDARD SEQ. 

(WIRE ELEMENT CURRENTS DERIUED APPROPRIATELY) 

THE SIGN FOR CURRENT IN SYMMETRY-DERIMED E/MB IS SUCH THAT 
THE COMPONENT OF FIELD RADIALLY IN TOWARDS THE ORIGIN IS 
MAINTAINED CONSTANT (IN MAGNITUDE. NOT DIRECTION) WITH 
CONSTANT E/M CURRENTS 

<DATA) OUTPUT OPTION <IOUT) 

=AB where: 

A=FIELD LISTING OPTION 

0 NONE 

1 TO CONSOLE 

2 TO DATAFILE (2FIELD.DAT) 

B-FORCE/TORQUE OUTPUT OPTION 

0 TO CONSOLE 

1 TO DATAFILE < 2F0RCE.DAT) 

EG. 21 ROUTES FIELD AND FORCE/TORHUE DATA TO FILES 

0 ROUTES FORCE/TORUUE DATA TO CONSOLE. FIELD DATA LOST 


Refer to subroutine principal variable lists for description of 
detail input variables. 
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C ROOT SEGMENT OF PROGRAM FORCE 

C NO VARIABLES FUNCTIONAL IN THIS SEGMENT 

C 

0001 VIRTUAL CX(1000j 2> » CY < 1000 y 2 ) j CZ ( 1000 r 2 ) > CORE < 1000 ? 6 ) 

0002 D I MENS I ON D ISP < 7 ) ? FORCE ( 9 ) y I SYMM ( 5 r 3 ) ? CURR ( 1 6 s’ 5 ) , I D I SP ( 3 ) 
C 

C MODEL DATA INPUT 
C 

0003 CALL MODINICORErDISP? IDISP) 

C 

C COIL DATA INPUT 
C 

0004 CALL COIL (CX.CYfCZ»ISYMMy CURR) 

C 

C FORCE S’ TORQUE AND FIELD CALCULATION 
C 

0005 CALL CALC<CX.CY7CZs-CORE?DISp5’ISYMMs’CURRs-IDISP) 

0006 STOP ' HAVE A NICE DAY' 

0007 END 


A2a6 Subroutine listings 
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0001 

0002 

0003 

0004 


0005 

0006 

0007 

0008 

0009 

0010 


0011 

0012 

0013 

0014 
0016 

0017 

0018 

0019 

0020 
0021 


0022 

0023 

0024 

0025 

0026 

0027 

0028 
0029 


0030 

0031 

0032 


0033 

0034 

0035 

0036 

0037 

0038 

0039 


SUBROUTINE MODIN(COREi.DISP»;[Dir>P) 

VIRTUAL C0RE(1000»6) 

DIMENSION EiIGP<7) f IDISP(3) 

PI=3. 141592654 
C 

C ZEROISE CORE AND DISFL. ARRAYS 
C 

DO 10 J=l»3 
10 IDISP<J)=0 

DO 20 1=1 » 6 
DO 20 J=l»1000 
20 C0RE<J»I)=0* 

NE=0 

C 

C MODEL INPUT SECTION 
C 

30 TYPE 40 

40 FORMAT!' MODEL GEOMETRY OPTION?') 

READ(5»)KrERR---30) MOPT 
IF(MOPT.EO.O) GO TO 60 
M0PT2=M0PT/10 
M0PT1=M0PT-M0PT2*10 
CALL FUS I N ( MOPT 1 » NE ). CORE) 

CALL W 1 NG < M0PT2 . NE r CORE ) 

IDISP<3)=NE 
GO TO 50 
C 

C MODEL INPUT FROM DATAFILE 
C 

60 TYPE 70 

70 FORMAT!' FILENAME?'/) 

CALL ASSIGN! 1,»-1» 'OLD' »») 

DEFINE FILE 1 ! 1001 » 12»U»NREC) 

DO 80 1=1 » 1000 

80 READ ! 1 ' I ) ! CORE ! I » J ) » J=1 » 6 ) 

READ! 1' 1001) IDISP!3) 

CLOSE ! UNI T=l) 

C 

C MODEL OUTPUT OPTION 
C 

SO TYPE 90 

90 FORMAT!' MODEL OUTPUT?') 

READ ! 5 » * r ERR=50 ) I OUTOP 
C 

C MODEL OUTPUT 
C 

GO T0!100»110rl20) » !I0UT0P+1) 

110 TYPE 130 

TYPE 140 

130 FORMAT!' MODEL COFtE COORDS. AND MAG.') 

1 40 FORMAT ! 4X » 1 1 IX » 7X » 1 NY » 7X » 1 HZ » 8X » 6H ! V ) X » 8X » 6H ! J*V ) Y 

1»8X!-6H! J*V)Z) 

DO 15!) I = lrIDISP!3) 

ISO TYPE 160» !CORE!I, J) » J=lr6) 
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0040 

0041 


0042 

0043 

0044 

0045 
004A 

0047 

0048 


0049 

0050 

0051 

0052 


0054 

0055 

0056 
0058 


0059 

0060 
0061 
0062 

0063 

0064 

0065 


0067 

0063 

0069 

0070 

0072 

0073 

0074 


160 F0RMAT(3F0.4,3E14.4) 

GO TO 100 
C 

C MODEL OUTPUT TO DATAFILE 
C 

120 TYPE 70 

CALL AGS I GN < 2 » , -1 » ' NEW ' , , ) 

DEFINE FILE 2 ( 1001 , 12»U»NREC) 

DO 170 1=1 » 1000 

170 UR I TE ( 2 ' I ) ( CORE (I y J > » J= 1 » 6 ) 

UR I TE < 2 ' 1 00 1 ) i DI SP < 3 ) 

CL0SE<UNIT=2) 

C 

C GET MODEL POSN. AND ATTITUDE 
C 

100 TYPE 180 

180 FORMAT <' MODEL OISP. OPTION?') 

READ<5.:K>ERR=100) MDISPO 
IF<MDISPO,GT.O) GO TO 190 
C 

C MDISP0=0y-2 
C 

DO 200 J=ly7 
200 DISP(J)=0. 

IF ( MDISPO. EQ. -2) GO TO 220 
GO TO 210 
C 

C MDISP0=1 
C 

190 TYPE 230 

230 FORMAK' MODEL DISP .( X » Y » Z ) y ORI ( YAW y PITCH y ROLL )?' ) 
READ ( 5 y * y ERR-^190 ) ( DISP ( J ) y J- 1 y 6 ) 

DISP(7)=0. 

DO 240 J=4y6 

2 40 D I G P < -.1) •= D I SI-' ( J ) * P I / 1 80 . 

IF<MDXGPO.EQ.l) GO TO 210 
C 

C MDISP0=2y-2 
C 

220 TYPE 250 

250 FORMAK' STEPSyPARAMETERySIZE?' > 

READ < 5 y * y ERR =220 ) I DI SP ( 1 ) y I D I SP ( 2 ) y D I SP ( 7 ) 
IF<IDISP<2) .LT.4) GO TO 210 
D I SP ( 7 ) =D I SP < 7 ) *P I / 1 80 . 

210 RETURN 
END 
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0001 

0002 

0003 

0004 


0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 


0016 

0017 

0018 
0019 
0021 


0022 

0023 

0024 

0025 

0026 


0027 

0028 


0029 

0030 

0031 

0032 

0033 

0034 

0035 


0036 


SUEiROUT I NE FUG I N ( MOPT 1 . NE » CORE ) 

OIRTUAL CORE< 1000.6) 

PI=3. 141592654 

GO TO <30.20.10.20)»(M0PT1+1) 

C 

C FIELD GRID SECTION 
C 

10 TYPE 40 

40 FORMAT<' XO.YO.ZO. XINC. YINC.ZINC. NOX . NOY . NOZ? ' ) 

READ (5.*) XO.YO.ZO.XINC.YINC.ZINC.NOX.NOYrNOZ 

DO SO IX=1.N0X 

DO 50 lY-l.NOY 

DO 50 IZ-l.NOZ 

NE:-=NEH 

CORE < NE y 1 ) =X0+ < IX-1 ) sKXINC 
CORE < NE y 2 ) =Y0+ ( I Y--1 ) >KYINC 
50 CORE ( NE y 3 ) ^=Z0+ <IZ-1 ) StZINC 

GO TO 30 
C 

C CYLINDRICAL SECTION 
C 

20 TYPE 60 

60 FORMAK' CYLIND.FUSE, IFLENyFRADyNAyNR?' ) 

READ < 5 y )iO FLEN y FRAD y NL y NA y NR 
IF<NL.EQ.O) GO TO 80 
ANGL=PI>K2./NA 
C 

C LOOP FOR RADIUS 
C 

DO 70 IRADt=:OyNR-l 
R2-= ( IRAD+1 ) JKFRAD/NR 
R1^=^IRAD*FRAD/NR 
XAREA-:--ANGL« ( R2>kR2-Rl>!<Rl )/2. 

D-4.)!c<R2:i<)i<3--Rl!i<>K3))|<SIN(ANGL/2, )/(3.)K<R2)!(R2-Rl»!Rl )>kANGL) 
C 

C LOOP FOR ANGLE 
C 

DO 70 IANG=OyNA-l 
ANG^=IANG>kANGL 
C 

C LOOP FOR AXIAL POS. 

C 

DO 70 ILEN=OyNL-l 
NE=NE+1 

CORE < NE y 1 ) ^=-FLEN/2 . + < ILENiO . 5 ) ikFLEN/NL 
CORE ( NE y 2 ) = D*COS < ANG ) 

C0RE(NEy3)™D:l<SIN(ANG) 

CORE ( NE y 4 ) =--XAREA»!FLEN/NL 
70 CONTINUE 

C 

C END OF CYLINDRICAL SECTION 
C ELLIPSOID ENDCAPS SECTION 
C 

80 IFCMOPTl .NE.3) GO TO 30 
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oo;?8 

0039 

0040 

0041 


0042 

0043 

0044 

0045 

0046 


0047 

0048 


0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 


0061 

0062 


TYPE 90 

90 FORMAT<' ELLIP.ENDCAPSJARrNAl?' ) 

REA0(5»*) ARfNAl 
THETA=PI/<NA1!)(2) 

C 

C LOOP FOR ANGLE FROM AXIS 
C 

DO 100 ITHETA=0,NA1-1 

THETA1==THETA>KITHETA 

THETA2==TH£TA)i:<ITHETA+l) 

D==COS ( THETA 1 ) --COS < THETA2) 

T:=C0S(THETA1*2. )-C0S(THETA2)K2, ) 

C 

C LOOP FOR ANGLE AROUND AXIS 
C 

DO 100 IANG=0»NA-1 
ANG=IANG)KANGL 
C 

G LOOP FOR FRONT/BACK AND PSEUDO RADIUS 
C 

DO 100 IFRBAK=~1»1»2 
DO 100 IRAD^=0»NR-1 

ri=irad;icfrad/nr 

R2= < IRAD+1 ) *FRAD/NR 

R--= < R2:«):<4-Rl!i:*4 ) / ( R2**3-Rl>Kii(3) 

N£==NE+1 

03 . !»:SI N ( ANGL/2 . ) )KR* ( THETA+T ) / ( 4 . )!(ANGUD ) 

CORE ( NE f 1 ) -IFRBAK5K < -FLEN/2 * ~AR*Ri}:T:4!3 . / ( D>K1 6 . ) ) 
C0RE(NE»2)=C!itC0S(ANG) 

C0RE<NEf3):=C)KSIN(ANG) 

CORE ( NE » 4 ) = < R2>J!*3-R1 i»:*3 ) »ANGL>it AR*D/3 . 

100 CONTINUE 
C 

C END OF ELLIP. ENDCAPS 
C 

30 RETURN 
END 
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0001 

0002 

ooo;? 

0004 


0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 


0016 


0017 

0018 


0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 


SUBROUTINE W I NG ( M0PT2 » NE » CORE > 

OIRTUAL C:0RE(1000y6) 

PI=3. 141592654 
IF<M0PT2.E«.0) GO TO 60 
C 

C SLAB WING INPUT 
C 

10 TYPE 20 

20 FORMAT ( ' SLAB WING J 1/2SPAN y 1/2GAP y CHORD y THICK y TAPER y SWEEP? ' ) 

READ ( 5 y )K y ERR-l 0 ) SPAN y FRAD y CHORD y TH I CK y TAPER y SWEEP 
30 TYPE 40 

40 FORMATC' DLONG y DOERT y NSy NCy NT y I WIDIR? ' ) 

READ < 5 y >1! y ERR-30 ) DLONG y OVERT y NS y NC y NT y I WIDIR 
SWEEP:^-=SWEEP*P 1/180. 

SINSWE^'=SIN(SWE£P) 

C0SWE=T:0S< SWEEP) 

SPAN2“SPAN-FRAB 

C 

C LOOP FOR PORT AND STB. WINGS 
C 

DO 50 IWING=-lyly2 
C 

C LOOP FOR SPAN AND THICKNESS 
C 

DO 50 ISPAN=OyNS“l 
DO 50 ITHICK=OyNT-l 
C 

C CALCULATE CORE DIPOLES 
C 

DO 50 ICHORD=OyNC-l 
NE-.^NE-H 

T= ( < 1 . -TAPER ) * ( ISPAN+1 )-NS ) / ( < 1 . -TAPER ) >itISPAN-NS ) 

TYPO-FLOAT ( I SPAN ) /NS+ < 1 . -2 . » ( 1-T ) /3 . ) / < ( 1 . IT ) *NS ) 
CHORDlj:=CHORD* <1 . - < 1 . -TAPER ) «TYPO ) 

OOL-CHORD>KTI 1 1 CK*SPAN2)>: ( 1 . TT ) / ( 2 . >!<NC:i<NT^NS ) 

CORE ( NE y 1 ) =DLONG-CHORDUft < 0 . 5- ( 0 . 5+ ICHORD ) /NC ) -SPAN2:«TYP0 
1-YSINSWE/COSWE 

CORE ( NE y 2 ) == ( SP AN2)i<T YPO+FRAD ) * I W I NG 

CORE < NE y 3 ) ==DOERT-TI LICK* < 0 . 5- < 0 . 5+ 1 TH I CK ) /NT ) 

CORE ( NE y 4 )'=:VOL*SINSWE*IWiriIR 
C0RE<NEy5>=:00L;kC0SWE*IWIDIR*IWING 
50 CONTINUE 

60 RETURN 

END 
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0001 


SUBROUT INE COIL < CX , CY » CZ » ISYMM r CURR ) 

0002 


VIRTUAL CX(1000>2) >CY(1000f2) .CZ(1000»2) 

0003 


niMENSlON XX<4» 10) fYY(4»10) »ZZ<4» 10) > ISYMM(5r3) i.CURR( 16»5) 

0004 


niMENSION IN0C0<12> 

0005 


DATA INOCO/1 »4>2»4»6»8»2f8f4f 8> 12r 16/ 

OOOA 


NSET-0 

0007 


IWIRE=1 

0008 


IOUTOP-0 

0009 


PI=3.141S92654 

0010 


no 10 J=l»5 

0011 


ISYMM<J»1)«0 

0012 


no 10 1 = 1 f 16 

0013 

10 

CURR<I»J)=1. 

0014 

20 

TYPE 30 

0015 

30 

FORMAK' NO. E/M SETS?') 

0016 


REAn<S»*»ERR=20) NSET 

0017 

c 

IF(NSET.GT.O) GO TO 40 


C COMPLETE SYSTEM riATAFILE INPUT 
C 


0019 


TYPE SO 

0020 

50 

FORMAK ' FILENAME?'/) 

0021 


CALL ASSIGN(lr»-l»'OLn'»f) 

0022 


DEFINE FILE 1 < 1003» 12»U»M) 

0023 


no 60 1=1 »3 

0024 

60 

READ ( 1 ' I ) ( ISYMM ( J r I ) » J=1 » 5 ) 

0025 


no 70 K=l»1000 

0026 

70 

REAna'K+3) (CX<Kf J) »CY(K»J) rCZ(Ki- J) r J-1 »2) 

0027 


CL0SE<UNIT=1) 

0028 


DO 80 1=1 »5 

0029 


IFdSYMMdf 1) .EQ.O) GO TO 80 

0031 


NSFIT=NSET+1 

0032 

80 

CONTINUE 

0033 


GO TO 90 


L 

C E/M 
C 

40 

C 

C GET 

SET LOOP START 

0034 

no 100 ISETI=1,NSET 


MASTER E/M GEOMETRY FOR THIS SET 

0035 

L 

TYPE llOvISETl 

0036 

110 

FORMAT<' E/M SET'dS) 

0037 

120 

TYPE 130 

0038 

130 

FORMAK' INPUT OPTION»SYMMETRY OPTION?') 

0039 


READ ( 5 y !« y ERR=120) INPOPT y ISYMMdSETI y 1 ) 

0040 


GO TO (140yl50yl60yl70) INPOPT 


C 

C SINGLE LOOP INPUT SECTION 

c; 


0041 

140 

TYPE 180 

0042 

180 

FORMAT<' NO. ELEMENTS?') 

0043 


READ < 5 y » y ERR=1 40 ) N INT 

0044 

190 

TYPE 200 
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00^5 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 


0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 
0069 


0070 

0071 

0072 

0073 

0074 


0075 

0076 


0077 

0078 

0079 


0080 

0081 

0082 


200 FORMAK' ELEMENT END POINTS?') 

REAri(5r*»l-;RR=190) ( CX (I » 1 ) » CY < I » 1 ) fCZ < I » 1 ) » 1 = 1 i-NINT ) 

NEND=NINT+IWIRE~1 

DO 210 J.---IWIREf (NEND-1) 

CX(Jf2)=CX(J+l»l) 

CYC Jr2)=CY( J+lrl) 

210 CZ( J»2)=CZ(JH»1) 

CX(NEND»2)=CX(IWIREi>l) 

CYC NEND f 2 ) =CY ( lUlRE » 1 ) 

CZ C NEND , 2 ) =CZ ( I W I RE » 1 ) 

ISYMMC ISETI»3)=1 
GO TO 220 
C 

C INTERFACE ARRAY INPUT SECTION 
C 

150 TYPE 230 

230 FORMAT (' NUMDER OF LAYERS, LOOPS/LAYER, INTERFACES?') 

READ ( 5 , * , ERR=150 ) MM , NN , MINT 
240 TYPE 250 

250 FORMATC' INTERFACE COORDINATES?' ) 

DO 260 J=1,NINT 

260 READ ( 5 , * y ERR=240 ) ( XX < I , J ) y YY ( I y J ) , ZZ ( I y J > y 1 = 1 y 3 ) 

NEND=MM*NN*NINTfIUIRE-l 
CALL FACE ( MINT y MM y NN y XX y CX y IWIRE ) 

CALL FACE ( NINT y MM y NN y YY y CY y IWIRE > 

CALL FACE ( NINT y MM y NN y ZZ y CZ y IWIRE ) 

ISYMM<ISETIy3)=MM*NN 
GO TO 220 
C 

C END OF INTERFACE ARRAY 
C LATERAL CIRCULAR E/M SECTION 
C 

160 TYPE 280 

280 FORMATC ' DXyDYyRADl yRAD2yDYl yDY2yNDI0RyNDIVYyNSEGyANGLE?' > 

READ ( 5 y * y ERR=1 60 ) DX y DY y RADI y RAD2 y DYl y DY2 y NDI VR y NDI VY y NSEG y ANGLE 
C0SANG=C0S(ANGLE)|fPI/180. ) 

S I NANG=S I N C ANGLE*P I / 1 80 . ) 

C 

C LOOP FOR AXIAL DEPTH 
C 

DO 290 I==lyNDIVY 
Y=DY1 +( I-O . 5 ) * ( DY2 -DYl ) 7NDI VY 
C 

C LOOP FOR RADIAL THICKNESS 
C 

DO 290 M=lyNDIVR 

R=RAD1 f ( M-0 . 5 ) * < RAD2-RAD1 ) /NDI VR 
R=R*NSEG*SIN C PI /NSEG ) / ( PDKCOS ( PI/NSEG ) ) 

C 

C LOOP AROUND CIRCUMFERENCE 
C 

DO 290 J=lyNSEG 

K=< < I-l )>!<NDIVR+M-1 )*NSEGiJ+IWIRE-l 
ANG=PI)f(2./NSEG 
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0083 ZT1“-SIN(ANG*< J--1) )«R 

0084 ZT2 = SIN(ANG)!(J)*R 

0085 CX ( K y 1 ) =DX+COS < ANG* ( J~1 > ) >KR 

008A CX(Ky2)=DX + C0S(ANG.’)tJ))t!R 

0087 CY(Ky 1 )=Y!XOOSANG--ZT1*SINANG+DY 

0088 OY < K . 2 ) ==Y:i:C03ANG-ZT2>ltSINANG+DY 

0089 CZ<K'y 1 )=Y)t!SINANGTZTl*COSANG 

0090 CZ(K»2)=Y)i:S]:NANG+ZT2)t:C0SANG 

0091 290 CONTINUE 

0092 NENIi==NSEG;KNDIUR:«!NIiIVY+IUIRE-l 

0093 ISYMMaSETIy3:>=NDIOR*NDIOY 

0094 GO TO 220 
C 

C END OF LATERAL CIRCULAR E/H 
C AXIAL CIRCULAR E/M SECTION 
C 

0095 170 TYPE 310 

0096 310 FORMATC' XI y X2 y RADI y RAD2 y NDIUR y NDIOX y NSEG? ' ) 

0097 READ(5y*yERR=170) XI y X2 y RADI y RAD2 y NDI OR y NDIOX y NGEQ 
C 

C LOOP FOR AXIAL DEPTH 
C 

0098 DO 320 I==lyNDIOX 

0099 X=-Xl + < I-0.5)^:<X2~X1)/NDIUX 
C 

C LOOP FOR RADIAL THICKNESS 
C 

0100 DO 320 M=lyNDIOR 

0101 R:=RAD1 f (M-0.5))K(RAn2-RADl>/NDI0R 

0102 R=R;«NSEG*GIN<PI/NSEG)/<Pl!icCOS(PI/NSEG) ) 

C 

C LOOP AROUND CIRCUMFERENCE 
C 

0103 DO 320 J-lyNSEG 

0104 K=< < I-DiKNOIOR+M-DKNSEG+JTIWIRE-l 

0105 ANG=PI?K2./NSEG 

0106 CX<Kyl)=X 

0107 CX<K,2)=X 

0108 CY(Ky 1 )=COS(ANG*( J-U )>i<R 

0 109 C Y < K y 2 ) =COS < ANG>!( J ) *R 

0110 CZ < K y 1 ) =SIN ( ANG* < J -1 ) ) HtR 

0 1 1 1 CZ < K y 2 ) ^=S I N ( ANG* J ) *R 

0112 320 CONTINUE 

0113 N E N D == N S £ Ci * N D I M X ;« N D 1 0 R + 1 U I R E - 1 

0114 ISYMM< ISETIy3)=NDI0X)KNDI0R 
C 

C END OF AXIAL CIRCULAR E/M 
C STORE E/M ELEMENT NO, 

C 

0115 220 ISYMM< ISETI y2)=NEND-IWIRE+l 

C 

C OUTPUT SECTION 
C 

0116 IF(I0UT0P,Ea,3) GO TO 435 
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0118 

Oil'? 

0120 

0121 

0122 

0123 

0124 

0125 
012A 

0127 

0128 
0129 

0131 

0132 

0133 

0134 

0135 

0136 

0137 


0138 

0139 


0140 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 


0150 

0151 

0 152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 


340 TYPE 350 

350 FORMAT!' E/M OUTPUT'?') 

READ ( 5 y )K » ERR=340 ) I OUTOP 

GO TO (435»360»370»435), (IOUTOP+1) 

370 TYPE 380 
TYPE 381 

380 FORMAT!' WIRE COORDINATES') 

38 1 FORMAT ! 5X ? 2HX1 . 8X 1 2HY1 r 8X » 2HZ1 f 8X 1 2HX2 r 8X 1 2HY2 » 8X » 2HZ2 ) 
DO 390 I==IWIRE»NEND 

390 TYPE 400>!CX!I,J)»CY!I»J)»CZ!I»J)»J==1»2) 

400 F0RMAT!6F10.5) 

360 IF!INP0PT.NE.2) GO TO 435 
TYPE 410 
TYPE 411 

410 FORMAT!' INTERFACE NO. AND COORDS') 

4 1 1 FORMAT ! 3X » 4HINT . y 9X » IHX » 14X » IHY » 14X » IHZ ) 

DO 420 J-^--l»NINT 

420 TYPE 430f!JyXX!IyJ)yYY!IyJ)»ZZ!IrJ)fI = ls-4) 

430 FORMAT! I5»3F15. 5) 

C 

C INCREMENT WIRE ELEMENT COUNTER 
C 

435 IUIRE^=NEND+1 
100 CONTINUE 
C 

C END OF E/M SET LOOP 

C COMPLETE SYSTEM DATAFILE OUTPUT 

C 

440 IF!I0UT0P.NE.3) 60 TO 90 
TYPE 50 

CALL ASSIGN!2r y-1 y 'NEW' f y ) 

DEFINE FILE 2 ! 1003 y 12 y U y M ) 

DO 450 1=1 y 3 

450 WR I TE ! 2 ' I ) ! I SYMM ! J y I ) y J=1 y 5 ) 

DO 460 K=lyl000 

46!) WRITE!2'K+3) !CX!Ky J) yCY!KyJ) yCZ! Ky J) y J=1 y2) 

CL0SE!UNIT=2) 

C 

C READ AHEAD E/M CURRENTS 
C 

90 DO 470 ISETI=lyNSET 

NOC--INOCO! ISYMM! ISETI y 1 ) ) 

TYPE 1 10 y ISETI 
490 TYPE SOOyNOC 

500 FORMAT!' E/M AMP . TURNS/1000? !'yI2y')') 

READ ! 5 y )K y ERR=490 ) ! CURR ! I y ISETI ) y 1 = 1 y NOC ) 

DO 470 1=1 y NOC 

CURR! I y ISETI )=CURR! I y ISETI )*1000, 

470 CONTINUE 
RETURN 
END 
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0001 

0002 

0003 

0004 

0005 

0006 
0007 
00 OB 


SUBROUTINE CALC ( CX »CY r CZ f CORE y DISP y ISYMM y CURR y IDISP) 

UIRTUAI... CX( 1000y2) yCY(1000y2) yC7.< 1000y2) yCORE< 1000»6> 

n I MENS I ON POB ( 6 ) y BB ( 3 ) y GB < 9 ) y COSR < 3 ) y G I NR < 3 ) y T ( 3 ) y BPOS < 3 > 

DIMENSION DISP (7) yP0RCE<6) y ISYMM(5y 3) y CURRv 16y5) y IDISP (3) 

PI =3. 141592654 

XMU=0. 00000125664 

LUN0UT==5 

LUMFLD=:=5 


C OUTPUT ROUTING 
C OPEN DATAFILES IF REQUIRED 


0009 

10 

TYPE 20 

0010 

20 

FORMAT (' OUTPUT OPTION?') 

0011 


READ(5yKyERR-10) lOUT 

0012 


IFLIST--I0UT/1O 

0013 


1 0 1 ■ I L E 1 0 1 J T - 1 F L I S T K 1 0 

0014 


IF(IOFILE.NE.l) GO TO 30 

0016 


LUNOUT-1 

0017 


OPEN ( UNIT-1 y NAME- ' 2F0RCE . DAT 



lyRECORDSIZE-78) 

0018 

30 

IFaFLIST.NE.2) GO TO 40 

0020 


LUNFLD-2 

0021 


OPEN ( UN IT=2 y NAME- ' 2FIELD . DAT 



1 yRECORDSIZE-120) 

0022 

40 

I F < I FL. I ST . EQ . 0 . OR . LUNFLD . EQ . 

0024 


URITE(LUNFLDy60) 

0025 

60 

FORMA r < 1 HO y 5X y IHX y 9X y IHY y 9X y 



1 y 3HBXX y 7X y 31 IBXY y 7X y 3HBXZ y 7X y 


C POSITION/ORXENTATION LOOP START 

0026 

50 

DO 70 K-OyIDISP<l) 

0027 


DO 80 N-ly6 

0028 

30 

FORCE <N)-0. 


C CONVERT MODEL ORIENTATION 

0029 


DO 90 J-ly3 

0030 


COSR( J)=COS<DISP< JT3) ) 

0031 

90 

S I NR ( J ) =S I N ( D ISP ( Jf 3 ) ) 


C MODEL ELEMENT LOOP START 


c 


0032 

100 

C 

DO 110 N-lyIDISP(3) 


y TYPE=-- ' NEW ' y FORM== ' FORMATTED' 


■'NEW' yFORM=' FORMATTED' 


2) GO TO 50 

1 HZ y BX y 2HBX y 8X y 21 IBY y 8X y 2HBZ y 8X 


C ROTATE 
C 


ELEMENT 


0033 

0034 

0035 

0036 


DO 120 JJ=ly2 

POS ( J ) -CORE < N y J ) *COSR ( 1 ) 4:C0SR ( 2 ) - CORE < N y J f 1 ) >KSI NR ( 1 ) 

1 >KCOSR < 2 ) ( CORE < N y J+2 ) *SINR < 2 ) 

P 0 S, ( J H ) - C 0 R E ( N y J ) * ( S I N R ( 1 ) >K C 0 S R ( 3 ) -K; 0 S R ( 1 > >!< S I N R < 2 ) Ht 
1 S I NR < 3 ) H CORE < N y J f 1 ) K < COSR (1 ) )«COSR < 3 ) -S I NR (1 ) «S I NR ( 2 ) KB I NR ( 3 ) ) 
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PA(3E 002 


0037 


0038 

0039 


0040 


0041 

0042 

0043 

0044 


0045 

0046 

0047 


0048 

0050 

0051 

0052 


0053 

0055 

0056 

0057 

0058 


0059 

0060 


2TC0RE ( N . JT2 ) ;XCOSR ( 2) >XSINR ( 3 ) 

1 20 F'OG ( J + 2 ) =CORF ( N » J ) * ( SI NR ( 1 ) >XSINR ( 3 ) -COSR ( 1 ) INR ( 2 ) JftCOSR ( 3 ) ) 

l+CORE(N jJf 1 )?KCOSR< 1 )*SINR(3)+SINR( 1 ) *SINR < 2 ) *COSR ( 3 ) ) 

2iC0RF < N , J + 2 ) *COSR < 2 ) JXCOSR ( 3 ) 

C 

C GET MODEL ELEMENT ABSOLUTE POSITION 
C 

DO 130 J-=l»3 

130 DPOS( J)=POS( J)TDISP(J) 

C 

C GET FIELD FOR THIS ELEMENT 
C 

CALL F I ELD ( BB , GB » CX . CY » CZ r I SYMM , CURR , DPOS ) 

C 

C COMPUTE ELEMENT FORCES 
C 

DO 140 J=lf3 
J-1 ))X3+1 

T < J ) ■= ( POS < 4 ) » GB < J J ) EPOS ( 5 ) *GB < J JH ) EPOS ( 6 > >!<GB ( J JE2 ) ) 

1 40 FORCE ( J ) -FORCE ( J ) ET ( J ) /XMU 

C 

C AND ELEMENT TORQUES AND MOMENTS Oh FORCES 
C 

FORCE (. 4 ) -FORCE < 4 ) E ( < POS ( 2 ) *T ( 3 > -POS ( 3 ) »(T ( 2 ) ) E 
1 ( POS ( 5 ) )i;BB ( 3 ) -POS < 6 ) )XBB ( 2 ) ) ) /XMU 
FORCE < 5 ) ^= F ORCE ( 5 ) E < ( POS ( 3 ) >XT ( 1 ) -POS (1 ) >XT < 3 ) ) E 
1 < POS < 6 ) YBB a ) -POS ( 4 ) )XBB ( 3 ) ) ) /XMU 
FORCE ( 6 ) =FORCE ( 6 ) E ( < POS ( 1 ) *T ( 2 ) -POS < 2 ) )XT ( 1 ) ) E 
1 ( POS ( 4 ) EBB ( 2 ) -POS < 5 ) EBB ( 1 ) ) ) /XMU 
C 

C FIELD OUTPUTS? 

C 

IF< IFLIST.EQ.O) C70 TO 110 

WR I TE < LUNFLD y 1 50 ) POS ( 1 ) y POS ( 2 ) y POS < 3 ) y BB ( 1 ) y BB ( 2 > y BB < 3 ) 

1 y GB a ) y GB ( 2 ) y GB < 3 ) y GB < S ) y GB ( 6 ) y GB ( 9 ) 

150 FORMAT < 12F10. 4) 

110 CONTINUE 
C 

C END OF ELEMENT LOOP 
C FINAL OUTPUTS 
C 

IF<K.GT.O.OR.IOFILE.EQ.l) GO TO 160 
URITE(LUNOUTy 170) 

1 70 FORMAT ( 7X y 61 IFORCEX y 7X y 6HF0RCEY y 7X y 6HF0RCEZ y 7X y 6HT0RQUX y 7X 

1 y 6I1T0RQUY y 7X y olFIOROUZ ) 

160 WRITE (LUNOUTy 180) (FORCE ( I ) y 1=1 y6) 

180 F0RMAT(6F13.2) 

C 

C INCREMENT DISPLACEMENT 
C 

D I SF' ( I n I SI"' < 2 ) ) =D ISP ( I DI SP ( 2 ) ) ED ISP ( 7 ) 

70 CONTINUE 
C 

C END OF DISPLACEMENT LOOP 
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C 






0061 


IF(LUNOUT.NE.l) 

GO 

TO 

190 


0063 


CLOSE(UNIT==i) 





0064 

190 

IF<LUNFLD.N£.2) 

GO 

TO 

200 


0066 


CL0SE(UNIT==2) 





0067 

200 

RETURN 





0068 


END 






FORTRAN ly V02.5 PAGE 001 

0001 SUBF?0UT;CNE FACE(NINT,H»Ni.WW>W,IWIRE) 

c 

C THIS SUBROUTINE BERIOED FROM MIT PROGRAM “TABLE" » BUT WITH TWO 
0 SIGNIFICANT ALTERATIONS IN FUNCTION. REFER TO PROGRAM DESCRIPTION 


0002 

0003 

0004 


0005 

0006 


0007 

0008 
0009 
001.0 
0011 


0012 


0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 


C 

OIRTUAL W(1000»2) 

DIMENSION WW(4»10> 

IND"IWIRE"1 

C 

C GENERATE FOURTH CORNERS OF INTERFACES 
C 

DO 10 I==1.NINT 

10 WW<4y I)=WW<3>I)-WW(2»I)+WW(1»I) 

C 

C GENERATE INTERMEDIATE COORDINATES FOR EACH INTERFACE 
C 

DO 20 I^nyM 
DO 20 J=lrN 
DO 20 K==lyNINT 
IND:^=IND+1 

WdNDy 1 ) = (WW(lyK)*(2*<N-J> + l>H((2)K(M~I> + l)FWW(4yK)>!t(2>KJ-l> 
l)!((2)i<<M"-I) i lHWW(2»K))K(2>|cI-l ))K(2*(N--J) + l)iWW(3ylOi|t(2)KI--l> 
2iK<2))tJ-l ) )/<4*M)KN) 

20 CONTINUE 

C 

C GENERATE W2 FROM W1 SO THAT ENDPOINTS OF WIRES ARE TOGETHER 
C 

DO 30 I-^lyM)f<N 

J::=(I-1))KNINT 

W<I^<NINT»2)==W(Jilfl> 

DO 30 K=lyNINT-l 
W(JM<y2)==W( JFK+1»1) 

30 CONTINUE 

RETURN 
END 
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000 .1 SUBROUT I NE F" I rO Ji < BB » GB » CX i- CY i- CZ » I SYHM t CURR > DPOS > 

0002 VIRTUAL CX ( 1000 , 2 ) y CY < 1000 » 2 > » CZ < 1000 » 2 ) 

0003 B 1 ht::NS 1 ON S ( 3 ) f T ( 3 ) y DC ( 3 ) y DD ( 3 > y D6 ( 3 > y DP < 3 ) y DS ( 3 ) 

0004 DIMI- NBION BB< 3) yGB<6) yCURR< 16y5> r ISYMM<Sy 3) y DP0S(3> 

0005 IWIRE^-1 

0006 P I ^=3. 141592654 

0007 DO 10 N=ly9 

0003 10 GB(N)=0. 

0009 DO 20 N^=ly3 

0010 20 BB<N)=0. 

C 

C COIL SET LOOP START 
C 

0011 DO 30 I SET-1 y 5 

0012 INCOIL-0 

0013 IFdSYMMdSETyl) .EQ.O) GO TO 30 

0015 IREFL=1 

0016 IFaSYMMaSETyl) .LT.7) GO TO 40 

0018 IR'EFL—1 
C 

C CALCULATION VECTORING AND CONTROL SECTION 
C 

0019 40 DO 50 IX-lyIREFl.y-2 
C 

C MASTER COIL 
C 

0020 DY==1. 

0021 DZ-1. 

0022 DYZ=0. 

0023 INCOIL==INCOILK 

0024 I CURR- IX 

0025 ASSIGN 60 TO I POINT 

0026 GO TO 1000 

0027 60 IFdSYMMdSETyl) .EO.l.OR.ISYMM(ISETyl) .EQ.7) GO TO 50 

0029 IF<ISYMMdSETyl) .Ea.2.0R.ISYMM<ISETyl> .ECKG) GO TO 70 

C 

C ROTATIONAL SYMMETRIES 
C 

0031 IFdSYMMdSETyl) .GT.6) GO TO 80 

0033 I ROT- ( I SYMM ( I SET y 1 ) -2 ) *2 

0034 GO TO 90 

0035 30 IR0T-(ISYMM<ISETyl)~8)»2 

0036 90 DO 100 N-lvIROT-l 

0037 ANGL-2.)l!PIM<N/IR0T 

0038 DY-COS(ANGL) 

0039 DZ-DY 

0040 DYZ-SIN< ANGL) 

0041 INCOIL-INCOIL+1 

0042 I CURR- IX 

0043 ASSIGN 100 TO IPOINT 

0044 GO TO 1000 

0045 100 CONTINUE 

0046 GO TO 50 
C 
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C REFLECTION SYMMETRIES 
C 


0047 

70 

nz=-i. 

0048 


BY=1 . 

0049 


BYZ=0. 

0050 


INC0IL=INCGIL+1 

0051 


ICURR==-IX 

0052 


ASSIGN no TO IPOINT 

0053 


GO TO 1000 

0054 

no 

DY^=--1. 

0055 


nz=i . 

0056 


INCOIL-INCOIL+1 

0057 


ICURR.--IX 

0058 


ASSIGN 120 TO IPOINT 

0059 


GO TO 1000 

0060 

120 

DY==”1 , 

0061 


DZ=“1. 

0062 


INC0IL=INCD1L+1 

0063 


ICURR:^=IX 

0064 


ASSIGN 105 TO IPOINT 

0065 


GO TO 1000 

0066 

105 

CONTINUE 

0067 

so 

n 

CONTINUE 


C 

C INCREMENT WIRE COUNTER 


C 


0068 

I W I RE::~ I W I RET I SYMM ( I SET 

0069 30 

CONTINUE 


C 

C END OF con. SET LOOP 
C COMPLETE DUMMY VARIABLES 


C 


0070 

GB ( 4 ) =N3B < 2 ) 

0071 

GB(7)“GB<3) 

0072 

GB(8)=GB(6) 

0073 

RETURN 


C * >1< )»:** *)!< )l< >!<;!<>!<* ^ ** * >1< * Jit *** Jt; )|< ;!( )i( ***>!< :lt SK * >k JK* >k >;<«:)<>!< >1< :K >;<* >1< >!( >K * .'f: * ^ !i( )|! >!( >|c 

C 

C FIELD AT POINT SECTION 

i: 

C CORRECT WIRE CURRENTS 
C 

0074 :l, 000 CURREfP^= I CURRJiiCURR ( I NCOIL » ISET ) / 1 SYMM < I SET > S ) 

c; 

C THE REMAINDER OF THIS SUBROUTINE IS TAKEN WITH NO SIGNIFICANT 
C ALTERATION FROM THE PROGRAM " TABLE “ r COURTESY MASSACHUSETTS 
C INSTITUTE OF TECtINOLOGY. 

C 

C WIRE ELEMENT LOOP START 
C 

00/5 no 2000 L"IWIREyIWIREMSYMM(ISET»2)-l 

C 
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0076 

0077 

0078 

0079 

0080 
GOBI 


0082 

0083 

0084 
008 S 
0086 
0087 


0088 

0089 

0090 


0091 

0092 

0093 

0094 

0095 

0096 


0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 


0105 

0106 
0108 


0109 

0110 


C CAL.CULATE ArP»C»D»EfF 
C 

A" ( I X»CX < L f 1 ) -DPOS <1 ) ) 

< I X>KCX < L » 2 ) -DPOS (1 ) ) 

C"(DY1<CY(L» 1 )-DYZ!j<CZ(L»l )“DP0S<2) ) 

( DY)KCY ( L r 2 ) -H YZifcCZ (L » 2 ) -DPOS ( 2 ) ) 

E==(DZ)!<CZ<L» 1 )+DYZ>KOY<L» 1 )-tiP0S<3> ) 

F- ( DZKCZ ( L f 2 ) f-DYZikCY < L » 2 ) -DPOS < 3 ) ) 

C 

C SUBSCRIPT A»B>CrDfE.F FOR LATER USE 
C 

8(1) -^A 
S(2)-C 
8(3) ::=E 
T(1)=B 
T(2)=D 
T ( 3 ) ^F 
C 

C CALCULATE U»0»W 
C 

U:-C)!<F-DKE 

0«E*B-F*A 

W-A)KD-B)KC 

C 

C CALCULATE RHOly RH02 AND THEIR DOT» CROSS PRODUCTS 
C 

R1- ( A*AFC>KC+E)KE ) iloKO .5 
R 2 ■- ( B B I D t D f F * I-' ) ^! >K 0 . 5 
RS-R1+R2 
RM~4-;;l)!<R2 
RDR^=A!}<B(C>KD+E*F 
RXR=aHVFU 
C 

C CALCULATE THE DEKIUATIMES OF THE SUMy ETC. OF RHOl AND RH02 
C 

DO 1010 M^--ly3 

DP < h ) ( 8 ( M ) >KR2/R1 f T ( M ) >KR1/R2 ) 

DS ( M ) '■■■•- ( S ( M ) /R 1 f r ( M ) /R2 ) 

DD ( M ) -- ( S ( M ) FT ( M ) ) 

1010 CONTINUE 

DC( 1 )^=F-EFC-D 
DC(2)-E-FFD-A 
DC;<3)=D-CFA-D 
C 

C CALCULATE AND TEST H TO DETERMINE EON. FOR G TO BE USED 
C 

H== ( RMFRDR ) /RM 
IF(H.LT.O.Ol) GO TO 1020 
G-TiS/(RM)K(RM+RDR) ) 

C 

C CALCULATE G AND ITS XyYyZ DERIOATIVES 
C 

DO 1030 M-ly3 
DGA:=RM* ( RM+RDR ) KDS ( M ) 
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0111 


riGB --RS* < RMi< ( DP < M ) TDD < M ) ) TDP ( M ) * ( RMTFiDR ) ) 

0 1 1 2 

1 030 

DG < M ) ^ 

^ < DGA-DGD ) / ( RM* ( RMTRDR > ) **2 

0113 


GO TO 

1040 

0114 

1020 

( ( RS ) iX < RM-RDR ) ) / < RM*RXR*RXR > 

0115 


DO 1050 M==l>3 

0 1 1 6 


DGA - ( R3!X ( DP ( M > ~DD ( M ) ) TDS < M ) )X ( RM -RDR ) ) >XRM)t!RXRS!)K2 

0117 


DGB^^=R! 

( RM -RDR ) * ( RM*2 . *RXR*DC ( M ) TDP ( M ) *RXR**2 ) 

0 1 1 B 

1050 

DG(M)= 

- ( DGA-DGB ) / < RM*RXR**2 ) **2 


c: calculatih: 

THE FIELD AND GRADIENT CONTRIBUTIONS OF EACH 


C CURRENT ELEMENT 

0119 

1040 

DGX-=D(3(1) 

0120 


DGY=^DG<2) 

0121 


DGZ=DG<3) 

0122 


C U R F' = C U R R' E N / 1 0 0 0 0 0 0 0 . 

0123 


CURM^TJURPiXG 


C 

C SUM 

THE INDIOIDUAL CONTRIBUTIONS TO THE FIELD AND GRADIENT 


C GET 

THE TOTAL FIELD AND GRADIENT 

0124 

C 

DDU)^ 

-BB( 1 )TCURM*U 

0125 


DD<2):- 

--BB<2)TCURM*0 

0126 


BB ( 3 ) : 

-BB < 3 ) TCURM*W 

0127 


GB< 1 

T?(B<1HCURP*U!KDGX 

0128 


GB(2)^ 

=GB < 2 ) TCURP* ( G* ( E-F ) TU>ltDGY ) 

0129 


GB<3): 

=(3B ( 3 ) I CURP* < G* ( D-C ) TU*DGZ ) 

0130 


GB<5)= 

■^GB<S)TCURP*0*D6Y 

0131 


GB(6): 

OD ( 6 ) TCURP* < G* < A -B ) TU*DGZ ) 

0132 


GB<9): 

T3B(9)TCURP*DGZ>KW 

0133 

2000 

1"' 

CONTINUE 


L 

c (-:nd 

OF WIRE ELEMENT LOOP 

0134 

0135 

(2 

GO TO 
END 

IPOINT 
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PRINCIPAL VARIABLE LIST FOR SUBROUTINE MODIN *tttt*t**ttttt**t**t* 


variables passed into MODIN 

NONE 

VARIABLES PASSED OUT OF MODIN 

COKEfDISPfIDISP LOADED WITH MODEL CHARACTERISTICS 


CORE MODEL ELEMENT POSITION AND MAGNETIZATIONS 

<N.*l> X POSITION 

<N*2> Y 

<N*3) Z 

(N-4> X WISE MAGNETX2ATI0N*ELEMENT VOLUME 
(N.5) Y 

<Nf6> Z 

DISP REOUIREti MODEL DISPLACEMENT FROM ZERO POSITION 

<1> X DISPLACEMENT 

Y 

<3) r 

iA> YAW DISPLACEMENT 

(5> PITCH 

<A> ROLL 

<7) STEP SIZE FOR DISPLACEMENT LOOP WHERE REQUIRED 

IDISP LOOP CONTROL PARAMETERS 

<1) NO. OF REQUIRED STEPS IN POSITION/ORIENTATION LOOP 
(2> parameter (IN ARRAY DISP) TO BE INCREMENTED 

<3> TOTAL NO. OF MODEL ELEMENTS 

MOPT MODEL INPUT OPTION 

DECODED INTO 

MOPTl fUSELAGE INPUT OPTION 

M0PT2 WlNii INPUT OPTION 

NE NPnDFR OF MODEL ELEMENTS (PASSED TO IDISP<3>) 

ICUTCP MUDLL OUTPUT UPTION 

0 NONE 

1 LISTING ON CONSOLE DEVICE 

2 to DATAFILE 

MDISPO MODEL DISPLACEMENT OPTION 

0 NONE 

1 SINGLE SELECTED POSITION/ORIENTATION 

2 SU.ELIED INITIAL DISPLACEMENT WITH STEPPING OF ONE 
SCI crrru, pARAMrirR 

-2 AS 2 PUT NULL INITIAL DISPLACEMENT 



Principal variable lists 



principal variable list for subroutine wing **************************** 
variables passed into wing 

hOPT2 FROM MODIN 

NE SET POINTING AT LAST ENTRY IN CORE ARRAY 

CORE MAY BE PARTIALLY LOADED 

variables passed out of wind 

NE UPDATED WHERE APPROPRIATE 

CORE DITTO 


I 

to 

VO 

a\ 

I 


SPAN HALF WING SPAN 

FRAD HALF WING CENTRE GAP (FUSELAGE RADIUS?) 

CHORD WING ROOT CHORD 

THICK WING THICKNESS 

lA^ER TIP CHOKD/KOOT CHORD 

SWtEP WING SUECr-BACK 

DlONG X WISE rUSPLACEMENT OF WING FROM MODEL DATUM 

rVEKT 2 WISE DISPLACEMENT OF WING FROM MODEL DATUM 

(DLONGfDVERT ONLY USEFUL WHERE FUSELAGE PRESENT) 
NS NO. OF ELEMENTS ALONG SINGLE WING PANEL 

NC NO. OF LlEMENTS ALONG WING CHORD 

NT NO. OF elements THROUGH WIND THICKNESS 

IWXDIR WING MAGNETIZATION DIRECTION 

SINSWE SINE OF SWEEP ANGLE 

CCSWE COSINE OF SWEEP ANGLE 

SPAN2 WING PANEL LENGTH 

CHOKDL LOCAL CMUKD 

distance OF centroid of ELEMENT FROM WING ROOT 
VOL Ft EMFNT VOLUME 


EQUIVALENT VALUE OK TAPER FOR SINGLE ELEMENT 


PRINCIPAL VARIABLE LIST FOR SUBROUTINE FACE 

VARIABLES PASSED INTO FACE 
NINTrM»NfWUf IWIRE 


VARIABLES PASSED OUT OF FACE 

W FILLED WITH WIRE COORDINATES 


NINT NO. INTERFACES 

M NO. LAYERS OF WIRE 

N NO. UIRES/l.AYER 

WW INTERFACE COORDINATES <XrY OR 2) 

(IfN) CORNER 1 

<2»N) 2 

<3»N) 3 

<4fN) A (PROGRAM GENERATED) 

IWIRE POINTER TO FIRST FREE ENTRY IN W 

U GENERATED WIRE COORDINATES <XfY OR 2) 

<N*1) WIRE START 

<Nt2> WIRE END 



PRINCIP«5iL VARIABLE LIST FOR SUBROUTINE COIL *******!*! !*:**)*!******i»!*i|c*)lt**** 


VARIABLES PASSES INTO COIL 
NONE 

VARIABLES PASSED OUT OF COIL 

CXyCYrCZi'ISYMMjCURR E/M ARRAY DATA 


X»Y»Z COORDINATES OF WIRE ELEMENTS 
WIRE ELEMENT START 
WIRE ELEMENT END 
E/M SYMMETRY CODES AND SET DATA 
SYMMETRY OPTION FOR THIS SET 
NO, OF WIRE ELEMENTS PER E/M IN THIS SET 
NO. OF LOOPS OF ELEMENTS IN THIS SET 
E/M TOTAL AMPERE-TURNS 
IN FIRST E/M OF SET 
SECOND 
ETC. 

COORDINATES OF E/M INTERFACES (WHERE USED) 

FIRST CORNER 
SECOND 
THIRD 

FOURTH (PROGRAM GENERATED) 

NO. OF E/MS PER SET 

NO. IN SET WITH SYMMETRY OPTION N 

NO. OF E/M SETS 

POINTER TO E/M SET DATA 
POINTER TO FIRST WIRE ELEMENT OF E/M 
POINTER TO LAST WIRE ELEMENT OF E/M 
E/M OUTPUT OPTION 
E/M INPUT OPTION 
NO. OF E/MS IN THIS SET 

SOME VARIABLES ACTIVE IN THE SINGLE LOOP INPUT SECTION 
NINT NO. OF ELEMENTS 

SOME VARIABLES ACTIVE IN THE INTERFACE ARRAY INPUT SECTION 

MM»NN no. OF LAYERS » LOOPS/LAYER OF WIRE ELEMENTS 

NINT NO. OF INTERFACES 

SOME VARIABLES ACTIVE IN THE LATERAL CIRCULAR E/M INPUT SECTION 

BX.DY X»Y wise DISPLACEMENT OF ORIGIN OF E/M AXIS 

RAD1.RAD2 INNER, OUTER RADII OF E/M 

DY1,DY2 NORMAL DIST. FROM X AXIS TO E/M FRONT, REAR FACE 

NDIVR,NDIVY NO. OF RADIAL, AXIAL LOOPS 

NSEG NO. OF WIRE ELEMENTS/LOOP (CIRCUMFERENTIALLY) 

ANGLE ANGLE FROM Y AXIS (IN YZ PLANE) TO E/M AXIS 

K COUNTER OF WIRE ELEMENTS IN THIS E/M 

R RADIUS FROM E/M AXIS TO END OF WIRE ELEMENT 

SOME VARIABLES ACTIVE IN THE AXIAL CIRCULAR E/M INPUT SECTION 
XI, X2 AXIAL DISTANCE TO E/M INNER, OUTER FACE 

NDIVR,NDIVX NO. OF RADIAL, AXIAL LOOPS 

NSEG, RADI, 

PAD2,R,K AS LATFZRAL CIRCULAR E/M SECTION 


CX,CY,CZ 

(N,l) 

(N,2) 

ISYMM 

(N,l ) 
(N,2) 
(N,3) 

CURR 

(1,N) 

(2,N) 

ETC. 

XX,YY,ZZ 

(1,N) 

(2,N) 

(3,N) 

(4,N) 

INOCO 

(N) 

NSET 
I SE T I 
I WIRE 
NEND 

lourop 

INPOPT 

NOC 
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PRINCIPAL VARIABLE LIST FOR SUBROUTINE CALC *******1^******t**t***t*■¥***** 

VARIABLES PASSED INTO CALC 

CX,CY,CZ,ISYMM»CURR FROM COIL 

DISP»IDISP»CORE FROM MODIN 

VARIABLES PASSED OUT OF CALC 
NONE 


XMU PERMEABILITY OF FREE SPACE(SI) 

LUNOUT FORTRAN LOGICAL UNIT NUMBER FOR MAIN OUTPUTS 

LUNFLD FORTRAN LOGICAL UNIT NUMBER FOR FIELD LISTINGS 

lOUT OUTPUT OPTION 

DECODED INTO 

IFLIST =2 FIELD LISTING TO DATAFILE 

=1 FIELD LISTING ON CONSOLE DEVICE 

=0 NONE 

lOFILE =1 MAIN OUTPUTS TO DATAFILE 


=0 MAIN OUTPUTS TO CONSOLE DEVICE 


K 


DUMMY VARIABLE CONTROLLING 

POSITION/ORIENTATION LOOP 

FORCE 


MAIN RESULTS ARRAY 




(1) 

FORCE IN X DIRECTION 




<2) 

Y 




(3) 

2 




<4) 

TOROUE ABOUT X AXIS 




<5) 

Y 




<6> 

Z 



COSR 


COSINES OF MODEL ROTATIONS 




<1) 

YAW 




<2) 

PITCH 




(3) 

ROLL 



SI NR 


SINES OF MODEL ROTATIONS 




<1) 

YAW 




(2) 

PITCH 




<3) 

ROLL 



N 


MODEL ELEMENT NUMBER 



POS 


POSITION AND MAGNETIZATION 

COMPONENTS 

OF MODEL ELEMENT 


(1) 

X POSITION 




<2) 

Y 




(3) 

2 




<4) 

X WISE MAGNETIZATION 




<5) 

Y 




<6) 

Z 



DPOS 


ABSOLUTE POSITION OF TRANSLATED MODEL 

ELEMENT 


<1) 

X 




<2). 

Y 




<3) 

Z 






PRINCIPAL MARIABLE LIST FOR SUEiROUTlNE FIELD 

VARIABLES PAGBED INTO FIELD 

CXyCYrCZ»ISYMM»ClJRR FROM COIL 

BPOS FROM CALC 


VARIABLES PASSED OUT OF FIELD 

BByOB FIELD AT REOUIRED POINT 


BB 

GB 


ru 

•: 2 ) 

< 3 ) 


y X / 
(2) 
< 3 ) 
< 4 ) 
<S) 
( 6 ) 
(7) 
< 8 ) 
(9) 


I SET 
INCQIL 
I WIRE 
IREFL 

1 

.... I 

DXyDYyDYZr IX 
I ROT 
ICtJRR 
CURREN 


TOTAL FIELD AT POINT (TESLA) 

X COMPONENT 

Y 

I 

FIELD GRADIENTS AT POINT (TESLA/METRE) 

GRADIENT OF X FIELD IN X DIRECTION 
X Y 

X Z 

Y X 

Y Y 

Y Z 

Z X 

Z Y 

Z Z 

BY FREE SPACE FIELD EQUATIONS t 

(4)“<2) y (7)“(3) y (8)«(6) 
VECTOR FOR E/M SET DATA IN ISYMMy ETC. 

VECTOR FOR E/M DATA IN CURR 

VECTOR FOR WIRE ELEMENT DATA IN CXy ETC. 

AUXILIARY SYMMETRY CODE 

NO REFLECTION IN YZ PLANE 

REFLECTION IN YZ PLANE 

CONTROL PARAMETERS FOR REFLECTIONS AND ROTATIONS 
NO. OF E/M ROTATIONS ABOUT X AXIS 
CURRENT DIRECTION CORRECTION 
CURRENT IN WIRE ELEMENT 


IPOINT RETURN VECTOR FROM MIT DERIVED SUBPROGRAM 


VARIABLES SPECIFIC TO MIT DERIVED SUBPROGRAM 
NOT CLASSIFIED 
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> 

NJ 

CO 


.RUN 2FORCE 

fiODEU GEOMETRY OPTION? 

II 

CVLINO.FUSE. :FLENfFRAO*NA»NR? 

1 .0*0. 1 1 

SLAB UING: l/2SFAN#l/2GAP*CH0RPfTHICK*TAPERfSUEEP'? 

0. A. 0. : *0. 1*0.01*1 .0.0. 

PLONG » UVERT. NS * NC* NT* I UI PI R? 

0 . 0 . 3 . 2 * 1*1 

hOPfL OUTPUT? 

1 

hODEL CORE COORDS. AND MAG. 



X 

Y 

7 

( JTV)X 

< J*U)Y 

< J*V>Z 


-0.3333 

0.0600 

0.0000 

0.2618E-Q2 

O.OOOOE+00 

O.OOOOE+OO 


0.0000 

0.0300 ■ 

0.0000 

0.261SE-02 

O.OOOOE+00 

O.OOOOE+OO 


0.3333 

0.0600 

0.0000 

0.2618E-02 

O.OOOOE-fOO 

O.OOOOE+OO 

I 

-0.3333 

0.0000 

0.0600 

0. 2618E-02 

O.OOOOE+00 

O.OOOOE+OO 

LO 

0.0000 

0.0000 

0.0600 

0.2618E-02 

O.OOOOE+OO. 

O.OOOOE+OO 

0.3333 

0.0000 

0.0600 

0.2618E-02 

O.OOOOEIOO 

O.OOOOE+OO 

o 

-0.3333 

-0.0600 

0.0000 

0.261SE-02 

O.OOOOEIOO 

O.OOOOE+OO 


0.0000 

-0.0600 

0.0000 

0.261 BE -02 

O.OOOOEIOO 

O.OOOOE+OO 

1 

0.3333 

-0.0600 

0.0000 

0.261BL-02 

O.OOOOEl 00 

O.OOOOE+OO 


-0.3333 

0.0000 

-0.0600 

0.261SE-02 

O.OOOOEIOO 

O.OOOOE+OO 


0 . 0000 

0.0000 

-0.0600 

0.2v61SE-02 

O.OOOOEIOO 

O.OOOOE+OO 


0,3333 

0.0000 

-0.v)600 

0.26 IDE -02 

O.OOOOi; ♦ 00 

O.OOOOE+OO 


-0.0250 

-0. 1;^33 

0.0000 

O.OOOOi 100 

-0.8333E-04 

O.OOOOE+OO 


0.0250 

-0.1833 

0,0000 

O.OOOOEfOO 

-0. U333E- 04 

O.OOOOE+OO 


-0.02b0 

-0.3500 

0.0000 

O.OOOOE 100 

-0.83331-:- 04 

O.OOOOE+OO 


0.0250 

-0.3500 

0.0000 

O.OOOOEIQO 

-0.S333E-O4 

O.OOOOEIOO 


-0.0250 

-0.5167 

0.0000 

O.OOOOEIOO 

-0.C333E-04 

O.OOOOEIOO 


0.0250 

-0.5167 

0.0000 

O.OOOOEIOO 

-o.8;r3;5E-04 

O.OOOOE+OO 


-0.0250 

0. 1833 

0.0000 

O.OOOOt 100 

0.B333t-04 

O.OOOOE+OO 


0 . 0 250 

0. 1833 

0 . 0000 

O.OOOOi: 100 

O.U;33;‘lE- 04 

O.OOOOE+OO 


-0 . 0250 

0.3500 

0,0000 

0.00001 100 

0.8;?3:3E-<>4 

O.OOiOOE+OO 


0. 0250 

0 • 3500 

0.0000 

O.OOOOL 100 

0.U333E -04 

O.OOOOEIOO 


-0 . 0250 

0.51*7 

0.0000 

O.OOOOt 100 

0.8333E-04 

O.OOOOEIOO 


0.0250 0.5167 0.0000 

hOOEL DISP. OPTION? 

O.OOOOEIOO 

0.8333E-04 

O.OOOOE+OO 


1 


hODEL MSP. ( J(*Y.2) * ORH YAU *P I TCM. ROLL ) ? 
0*0*0*0*-?5. *0 
NO. E/rt SETS? 

I 

E//i SET 1 

ivnn OPTION*SYMhtlKY OPTION? 

3*10 

DX.DYfRADI *RAU2»DY1 *DY2*NDIVR*NItIVY*NSEG*ANGLE? 
2.0*0.* I .0*1 .S*I .5*2.0*2*2*8*4ti. 

E/M OlUPUT-? 


WIRE COORDINATES 


XI 

3.18663 
2.83908 
2.00000 
1 .16092 
0.81337 
1.16092 
2.00000 
2.83908 
3.45033 
3.02554 
2.00000 
0.97446 
0.54967 
0.97446 
2.00000 
3.02554 
3.18663 
2.83908 
2.00000 
1 .16092 
0.31337 
1 .16092 
2.00000 
2.83908 
3.45033 
3.02554 
2.00000 
0.97446 
0.54967 
0.97446 
2.00000 
3.02554 
E/M SET 
£/M ArtP.TI 


Y1 

1.14905 
0 . 55573 
0.30997 

0 • 55573 
1.14905 
1.74237 
1.98813 
1.74237 
1.14905 
0.42388 
0.12351 
0.42388 
1.14905 
1 .87421 
2.17459 
1.87421 
1.32583 
0.73251 
0.48675 
0.73251 
1.32583 
1.91914 
2. 16490 

1 .91914 
1.32583 
0.60066 
0.30029 
0.60066 
1 . 32583 
2.05099 
2.35136 
2.05099 

1 

NS/1000? 


Z1 

X2 

Y2 

22 

1.14905 

2.83908 

0 • 55573 

1.74237 

1.74237 

2.00000 

0.30997 

X . 9831 3 

1.98813 

1.16092 

0.55573 

1.74237 

1.74237 

0.81337 

1.14905 

1 .14905 

1 . 14905 

1 . 16092 

1.74237 

0.55573 

0.55573 

2.00000 

1.98313 

0.30997 

0,30997 

2.83908 

1.74237 

0.55573 

0.55573 

3.18663 

1.14905 

1 .14905 

1.14705 

3.02554 

0.42393 

1.37421 

1.87421 

2.00000 

0.12351 

2.17459 

2.17459 

0.97446 

0.42383 

1 .S/421 

1.87421 

0.54967 

1 .14905 

1 .14905 

1.14705 

0.97446 

1 .87421 

0.42383 

0. 42385 

2.00000 

2.17459 

0.12351 

0.12351 

3.02554 

1.87421 

0.42333 

0.42388 

3.45033 

1.14905 

1.14905 

1.32583 

2.83903 

0.73251 

1.91914 

1.91914 

2.00000 

0.4B675 

2.16490 

2.16490 

1 . 16092 

0.73251 

1 .91914 

1.91914 

0.81337 

1.32583 

1 .325B3 

1.32583 

1.16092 

1.91914 

0. 7:::25i 

0.73251 

2.00000 

2.16490 

0.48675 

0.48675 

2.83908 

1.91914 

0.73251 

0.73251 

3.18663 

1 .32583 

1.32583 

1 .32583 

3.02554 

0.60066 

2.05099 

2.05099 

2.00000 

0.30029 

2.35136 

2.35136 

0.97446 

0 . 60066 

■ 2.05099 

2.05099 

0.54967 

1.32583 

1 .32S8 3 

1 .32583 

0.77446 

2.05099 

0.60066 

0.60066 

2.00000 

2.35136 

0.30029 

0.30027 

3.02554 

2.05077 

0.60066 

0.60066 

3.45033 

1.32583 

1 .32583 


( 8 ) 


1000 . * 1000 . * - 1000 . *- 1000 . *- 1000 . *- 1000 . * 1000 . * 1000 . 

OUTPUT OPTION? 

10 


Example program run 
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B 

B etc 

• 

B 




X 

y 


XX 

X 

Y 

Z 

BX 

BY 

BZ 

BXX 

- 0.2357 

0.0600 

0.2357 

- 0.0443 

0.0004 

0,0429 

0.0050 

0.0000 

0.0600 

0.0000 

- 0.0000 

0.0000 

- 0.0000 

- 0.0000 

0.2357 

0.0600 

- 0.2357 

0.0443 

0.0004 

- 0.0429 

0.0050 

- 0.1933 

0.0000 

0.2731 

- 0.0518 

0.0000 

0.0346 

0.0053 

0.0424 

0.0000 

0.0424 

- 0.0079 

0.0000 

- 0.0079 

- 0.0002 

0.2781 

0.0000 

- 0.1933 

0.0366 

0.0000 

- 0.0513 

0.0042 

- 0.2357 

- 0.0600 

0.2357 

- 0.0443 

- 0.0004 

0.0429 

0.0050 

0.0000 

- 0.0600 

0.0000 

- 0.0000 

- 0.0000 

- 0.0000 

- 0.0000 

0.2357 

- 0.0600 

- 0.2357 

0.0443 

- 0.0004 

- 0.0429 

0.0050 

- 0.2781 

0.0000 

0.1933 

- 0.0366 

- 0.0000 

0.0513 

0.0042 

- 0.0424 

0.0000 

- 0.0424 

0.0079 

- 0.0000 

0.0079 

- 0.0002 

0.1933 

0.0000 

- 0.2781 

0.0518 

- 0.0000 

- 0.0346 

0.0053 

- 0.0177 

- 0.1833 

0.0177 

- 0.0033 

- 0.0000 

0.0033 

0.0000 

0.0177 

- 0.1833 

- 0.0177 

0.0033 

- 0.0000 

- 0.0033 

0.0000 

- 0.0177 

- 0.3500 

0.0177 

- 0.0034 

- 0.0000 

0.0034 

0.0001 

0.0177 

- 0.3500 

- 0.0177 

0.0034 

- 0.0000 

- 0.0034 

0.0001 

- 0.0177 

- 0.5167 

0.0177 

- 0.0035 

- 0.0000 

0.0035 

0.0001 

0.0177 

- 0.5167 

- 0.0177 

0.0035 

- 0.0000 

- 0.0035 

0.0001 

- 0.0177 

0. 1833 

0.0177 

- 0.0033 

0.0000 

0.0033 

0.0000 

0.0177 

0.1833 

- 0.0177 

0.0033 

0.0000 

- 0.0033 

0.0000 

- 0.0177 

0.3500 

0.0177 

- 0.0034 

0.0000 

0.0034 

0.0001 

0.0177 

0.3500 

- 0.0177 

0.0034 

0.0000 

- 0.0034 

0.0001 

- 0.0177 

0.5167 

0.0177 

- 0.0035 

0.0000 

0.0035 

0.0001 

0.0177 

0.5167 

- 0.0177 

0.0035 

0.0000 . 

- 0.0035 

0.0001 

FORCEX FORCEY 

FORCEZ 

TORQUX 

TORQUY 

3329 . 

85 

- 0.00 

- 3131.94 

0.00 


0.00 


STOP — HAVE A NICE HAY 


B etc. 

xy 


BXY 

BXZ 

BYY 

BYZ 

BZZ 

- 0.0021 

- 0.1844 

0.0072 

O.OOIS 

- 0.0122 

- 0.0000 

- 0.1873 

- 0.0000 

0.0000 

- 0.0000 

0.0021 

- 0.1844 

0.0072 

- 0.0018 

- 0.0122 

- 0.0000 

- 0.1808 

0.0069 

0.0000 

- 0.0122 

- 0,0000 

- 0.1370 

- 0.0002 

- 0.0000 

0.0004 

- 0.0000 

- 0.1369 

0.0072 

- 0.0000 

- 0.0115 

0.0021 

- 0. 1844 

0.0072 

- 0.0013 

- 0.0122 

- 0.0000 

- 0.1873 

0.0000 

- 0.0000 

- 0.0000 

- 0.0021 

- 0.1844 

0 . C 072 

0.0013 

- 0.0122 

- 0.0000 

- 0. 1869 

0.0072 

0.0000 

- 0.0115 

- 0.0000 

- 0.1870 

- 0,0002 

0.0000 

0.0004 

0.0000 

- 0.1808 

0.0069 

- 0.0000 

- 0.0122 

0 . 0004 

- 0.1891 

0.0000 

- 0.0004 

- 0.0001 

- 0.0004 

- 0.1891 

0,0000 

0.0004 

- 0.0001 

0.0006 

- 0.1939 

0.0000 

- 0.0006 

- 0.0001 

- 0.0006 

- 0.1939 

0.0000 

0.0006 

- 0.0001 

0.0007 

- 0.2003 

- 0.0000 

- 0.0007 

- 0.0001 

- 0.0007 

- 0.2003 

- 0.0000 

0.0007 

- 0.0001 

- 0.0004 

- 0.1891 

0.0000 

0.0004 

- 0.0001 

0.0004 

- 0.1891 

0.0000 

- 0.0004 

- 0.0001 

- 0.0006 

- 0.1939 

0.0000 

0.0006 

- 0.0001 

0.0006 

- 0.1939 

0.0000 

- 0.0006 

- 0.0001 

- 0.0007 

- 0.2003 

- 0.0000 

0.0007 

- 0.0001 

0.0007 

- 0.2003 

- 0.0000 

- 0.0007 

- 0.0001 


TORQUZ 


0.00 



A2.9 Diagnostics 


RT-11 FORTRAN generates WARNING diagnostics for the program 
segments MODIN, COIL and FIELD. These may be ignored 

Diagnostics for program vmit MODIN 

Warning in line 0026 Possible modification of index "I" 

Warning in line 0045 Possible modification of index "I" 

Diagnostics for program xinit COIL 

Warning in line 0023 Possible modification of index "I" 

Warning in line 0144 Possible modification of index "I" 


Diagnostics for program xinit FIELD 


Warning in line 0012 
Warning in line 0019 
Warning in line 0036 


DO loop entry 
DO loop entry 
DO loop entry 


at 

label 

"60" 

at 

label 

"60" 

at 

label 

"100" 
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APPENDIX 3 


Scaling of results to other physical sizes 

The majority of the theoretical analysis presented herein is 
unspecific as to the choice of physical scale. However, numerical 
computations have generally been made at particular fixed scales 
for convenience. Scaling of computed performance is thus of interest 
and is relatively straightforward provided some care is taken. 

Maximum force and torque capabilities 

Exact magnetostatic equivalence of the model (any core material) 
will only be obtained if the model's magnetic characteristics and the 
applied field remain equal at all corresponding points during scaling. 
The former condition may be taken for granted, the latter generally 
requires adjustment of E/M characteristics. With consistent E/M 
geometry and characteristics (notably current density) the field 
strength at geometrically similar points varies as; 

H a Scale 

Thus to preserve field equivalence, the E/M current density, for 
instance, could be varied as: 

J a Scale 

With such a condition it is clear that for fixed magnetic 
configurations: 

2 

Model forces a Scale (Vol.* Field gradient) 

3 

Model torques a Scale (Vol.* Field) 

This conflicts with the commonly quoted "Square-Cube" law (66 

and numerous others) , which implied that magnetic force capability 

varied as core volume (Scale^) whereas aerodynamic force varied as wing 
2 

area (Scale ). With regard to the maximum force and torque capability 
of MSBSs, the relations presented above (the "Square-Square" law?) 
are undoubtedly valid where permanent magnet cores are used, since here 
the performance of the core is limited by demagnetization, which in 
turn defends only on the total local field strengths and is not other- 
wise scale dependent. The relevance for soft iron cores is less 
obvious, since there is no particular restriction to operation of the 
core beyond saturation. However, when the core is far into saturation 
there is some evidence to support the validity of the Square-Square law 
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as a measure of maximum capability, for instance Section 5.5.1, 

3 

where torque capability was found to vary as Scale , implying of 

2 

course that force capability may vary as Scale . 

The absolute maximum capability for iron cores may prove to be 
unrealisable with exist ing technology E/Ms, as appears to be the 
case for the SIM computations in Section 5, thence the practical 
scaling law may lie between the Square-Cube and the Square-Square 
laws, inclusive. 

General 

Where maximum capabilities, core demagnetization etc. are not 
critical factors, the following conditions are easily identified 
(15, 49). 

For consistent balance geometry, model magnetizations, E/M 
current densities, etc.; 

Model forces n Scale^ 

4 

Model torques a Scale 

technological limitations on E/M performance are functions of 
scale, in particular the usable current density falls with increasing 
scale, but the trends are not continuous, apparently involving 
enforced abandonment of particular conductor technologies at specific 
limiting scales. Further treatment of these effects cannot be 
attempted here. 
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APPENDIX 4 


M agnetic units and definitions 

All equations are given in SI units, whence: 

B = y IH in free space (y = 4 x lO ^Hm B in Tesla; 

° ° H in A/m) 

However, two alternative subsystems exist, the Kennelley and 
Sommerfeld systems. The Kennelley system is used in this report, 
in the belief that it is somewhat more convenient where permanent 
magnetic material is present. The key definitions of this system are: 


^ = Hq ii + J. (flux through permanent magnet material) 

5f = J . V H^6v (force on a dipole) 

6t = X fl 6V (torque on a dipole) 

Both formulations lead acceptably to the classical: 

B = yyH j = yXH y = l+ x 

In this report J is replaced by M where the equations are relatively 
unspecific regarding choice of material (permanent magnet or iron core) . 
This is a non-standard symbology since the Sommerfeld system leads to 
(introducing primes to distinguish M' (traditional, in A/m) from M 
(Tesla) ) : 

B = y^ (H + M* ) 

It is clear that M in fact equals ■ The use of M partially 

avoids the unfortunate clash between the standard symbology for 
polarization (J) and current density (J) . There should be no con- 
textual difficulty in distinguishing the two otherwise. 

The first subscript to any field property describes the component 
under consideration, the second (where present) describes the gradient 
direction, e.g; 

H = field strength in direction a 

a 

= gradient of in the b direction 

In free space H , = H, 
ab ba 

Occasionally a gradient of a gradient is required, this being 
given thus: 

(H , ) s the rate of change (gradient) of the H , 
ab c ab 

component taken in the c direction. 
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APPENDIX 5 



It happens that classical applied fields in MSBSs with high levels 
of E/M symmetry may be identified and classified by the character- 
istic symmetry of those fields taken about the origin of axes and, 
strictly, in balance axes. This method of representing fields has 
been used extensively in Section 4 to derive relatively generalised 
expressions for the magnetic couplings occuring with simple SPM wing 
cores. 

When a model is displaced or rotated away from the usual datum 
position and orientation, the strict conditions of applied field 
symmetry break down. However, the adjustments to E/M currents that 
would be required to maintain the principal components of applied 
field, hence the generated force and torque vectors, as fixed in 
model axes, will tend to approximately preserve the symmetry 
conditions cited herein, also in model axes. Analysis based on 
applied field symmetry thus at least holds good for small displacements 
and rotations of the model. 

The E/M configuration employed in this Appendix is identical to 
that used in Section 4 for the SPM roll coupling study, and readers 
should refer to this Section (Figs. 4.8 - 4.10) for precise 
definition of this geometry. However, the fundamental results remain 
valid for any configuration of vaguely conventional form and all 
geometries thereof, provided the E/M array incorporates the required 
E/M symmetry. Fig. A5.1 thus illustrates the general E/M configura- 
tion with Table A5.1 showing the definition of the applied fields for 
this configuration. 

Using FORCE, the field and field gradient components around 
the origin of axes (at Hh 0,1m in x,y and z directions in fact) are 
easily found (Fig. AS. 2). The respective signs of the field and 
field gradient components are of principal interest, the magnitudes 
being included here for information only. 

The applied field symmetry matrices may be deduced directly from 
Fig. A5.2 and represented as matrices of Odd (O) or Even (E) 
functions of x,y and z, following Section 4, and Equations 4.8 and 
4.9. The full derived matrices are presented in Fig. A5.3. 
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Fig. A5. 1 Schematic diagram showing E/M configuration for field 

symmetry analysis. Refer to Figs, 4.8 - 4.10 for precise 
array geometry. 


Table A5.1 Definition of E/M current directi ons for classical 
applied fields. After Fig. A5.1 


Classical 

field 

Primary 

component 

I, I., I, I. I, 1^ 
123456789 10 

Lift 

Drag 

Sideforce 

Pitch 

Yaw 

Magnetizing 

(axial) 

B 

xz 3 
B 

B j 

B 

y B 

X 

+ 1 

+ + 

I + 

+ 1 

+ 1 

1 + 

+ + 

1 1 

I 1 

+ + 



I I I I 

11 12 13 14 

Roll 

B 

yz 

- + - + 


Current directions are chosen to make the classical primary field 
component for each applied field positive. 
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LIFT 
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(B ) 
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DRAG 

FIELD 

(B ) 

XX 
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(B ) 
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Position (m) 


X 

Y 

- 0.1000 

- 0.1000 

- 0.1000 

- 0.1000 

- 0.1000 

0.1000 

- 0.1000 

0.1000 

0.1000 

- 0.1000 
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- 0.1000 
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0.1000 

0.1000 


X 

Y 
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- 0.1000 
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- 0.1000 
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0.1000 

0.1000 


X 

Y 

- 0.1000 

- 0.1000 

- 0. 1000 

- 0.1000 

- 0.1000 

0.1000 

- 0.1000 
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0.1000 
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0.1000 

- 0.1000 

0.1000 

0.1000 

0.1000 

0.1000 


z 

- 0.1000 
0.1000 
- 0.1 000 
0.1000 
- 0.1000 
0.1000 
- 0.1000 
0.1000 


z 

- 0.1000 

0.1000 

- 0.1000 

0.1000 

- 0.1000 

0.1000 

- 0.1000 

0.1000 


z 

- 0.1000 

0.1000 

- 0.1000 

0.1000 

- 0.1000 

0.1000 

- 0.1000 

0.1000 


Field component (T) Field gradient component (T/m) 


PX 

tiY 

BZ 

PXX 

BXY 

- 0 . 02 A 2 

0.0002 

- 0.0244 

- 0.0022 

- 0.0016 

0.0242 

- 0.0002 

- 0.0244 

0.0022 

0.0016 

- 0.0242 

- 0. 0002 

- 0.0244 

- 0.0022 

0.001 6 

0.0242 

0.0002 

- 0.0244 

0.0022 

- 0.001 6 

- 0.0242 

- 0.0002 

0.0244 

0.0022 

- 0.0016 

0.0242 
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0.0244 

- 0,0022 

0.001 6 
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0.0016 
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0.0244 
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- 0.0016 


BX 

BY 

BZ 
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- 0.0226 
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- 0.0113 
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0 . 0226 
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0.0007 
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- 0.0113 
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- 0. 0113 
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- 0.0007 
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BYY 
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-0.0016 
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O 
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-0.0016 
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-0.0037 
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0.0016 

0 , 00 1 6 
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Ui 
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Ui 
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0 
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Cb 
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o 
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3 
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BX 

BY 

BZ 
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PITCH 

X 

- 0.1000 

Y 

- 0. 1000 

z 

- 0.1000 

BX 

0.0014 

BY 

- 0.0014 

FIELD 

- 0.1000 

- 0.1000 

0.1000 

- 0.0014 

0.0014 

- 0.1000 

0.1000 

- 0.1000 

0.0014 

0.0014 


- 0.1000 

0.1000 

0.1000 

- 0.0014 

- 0.0014 

0.1000 

- 0.1000 

- 0.1000 

- 0.0014 

- 0.0014 


0.1000 

- 0.1000 

0.1000 

0.0014 
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0.1000 

0.1000 

- 0.1000 

- 0.0014 

0.0014 


0.1000 

0.1000 

0.1000 

0.0014 

- 0.0014 


YAW 

X 

V 

z 

BX 

BY 

- 0.1000 

- 0.1000 

- 0.1000 

0.0014 

0.1865 

FIELD 
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0.1000 
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0 . 1 865 
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0.1000 
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0.1 365 
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Y 
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0.1000 
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0.1000 
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X 
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Y 
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z 
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BX 
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BY 
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0.0001 
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BXZ 
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BYZ 
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0.0137 
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Fig. A5.2 continued 
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Fig. A5.2 continued 



Fig. A5.3 Derived symmetry matrices for classical applied fields 
Following Figs. A5.1, A5.2 

The field symmetry is presented as in Eqns. 4.8, 4.9 
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